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INTRODUCTION 


The  Laser-Plasma  Interaction  group  at  NRL  is  evaluating  for  the 
Department  of  Energy  the  feasibility  of  using  direct  laser  drive  to 
implode  fusion  pellets. *  Mission  Research  Corporation  (W1C)  has 
contracted  to  support  this  experiment  by  using  its  best  effort  to  perform 
the  tasks  summarized  below: 

•  A  parametric  study  shall  be  performed  using  the  double  foil 
technique  on  actual  laser-foil  experiments.  Design  parameters,  to 
be  varied  as  needed,  include  laser  intensity,  material  and 
thickness  of  the  first  foil,  material  and  thickness  of  the  second 
foil,  and  the  spacing  between  foils.  Notched-edge  configurations 
will  also  be  investigated. 

•  In  support  of  the  above  task,  other  diagnostic  measurements  will  be 
performed.  Such  diagnostics  include  calorimeters,  Faraday  cups, 
and  ballistic  pendula. 

•  Analytical  techniques  and  methods  analyzing  the  experimental  data 
will  be  developed  and  Implemented  on  the  existing  PDP-11 /10  system. 

•  We  will  investigate  whether  the  x-ray  tracer  technique  ...  is 
suitable  to  track  plasma  flow  from  targets  accelerated  to  high 
velocities.  We  will  also  investigate  whether  this  technique  is 
suitable  for  the  study  of  plasma  hydrodynamic  instabilities.  For 
this  purpose  we  will  experiment  with  target  and  tracer  materials 
that  are  density  matched  and  at  the  same  time  maintain  a  sufficient 
tracer/target  x-ray  emission  contrast.  We  will  also  use  plasma 


*  B.  H.  Ripin  et  al.,  Phys.  Fluids  23(5),  1012  (1980);  24,  990  (E) 
(1981). 


tasks. 


calorimeters  and  Ion  collectors  to  investigate  laser  absorption  in 
long-density-scalelength  plasmas.  These  studies  will  include,  as 
needed,  parametric  variations  of  laser  Intensity,  laser  energy, 
laser  profile,  target  thickness,  and  target  density.  We  will 
operate  and  develop  software  for  the  4052  Tektronix  data 
acquisition  system  if  this  becomes  necessary  for  the  above 
experiments. 

Mission  Research  Corporation  (MRC)  will  interface  existing  NRL 
x-ray  detectors  and  particle  detectors  with  digital  data 
acquisition  systems  which  are  owned  by  htflL .  MRC  will  operate, 
maintain,  and  upgrade  these  systems  for  the  purpose  of 
investigating  the  properties  of  laser  interaction  with  matter. 
Software  will  be  developed  for  these  systems  if  necessary  for  the 
above  experiment. 

The  data  obtained  In  the  above  tasks  will  be  reduced  and  analyzed. 

The  following  sections  summarize  our  fulfillment  of  these 
Details  are  provided  in  the  appendices. 
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DOUBLE-FOIL  AND  NOTCHED-EDGE  EXPERIMENTS 


Experiments  using  double-foils  and  notched-edge  configurations 
have  been  performed.  In  particular: 

•  The  validity  of  the  double-foil  method  was  verified  using  clear 
double-foils  and  x-ray  backlighting  methods. 

•  A  low  density  plasma  traveling  ahead  of  the  accelerated  target  has 
been  identified  and  diagnosed. 

•  A  triple-foil  configuration  was  used  to  diagnose  the  impact  foil. 

•  Target  speeds  up  to  160  km/sec  have  been  measured. 

•  Target  pressure  has  been  measured  using  the  double- foil  technique. 
These  measurements  verify  earlier  results  obtained  with  ballistic 
pendula,  plasma  calorimeters,  and  ion  collectors. 

•  The  relationship  between  irradiance  uniformity  and  target  velocity 
uniformity  was  measured. 

•  A  notched  double-foil  configuration  was  used  to  diagnose  shocks  in 
the  impact  foil. 

Appendix  A  provides  more  details. 

ABLATION-PLASMA  ION  MEASUREMENTS 


Earlier  experiments  which  inferred  target  pressure  from 
measurements  on  the  ablation  plasma  were  compared  to  more  recent 
experiments  which  Inferred  target  pressure  from  its  acceleration  as 
measured  with  x-ray  backlighting.  Good  agreement  was  found. 

Appendices  B  an  C  provide  more  details. 
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ANALYSIS  OF  ION  MEASUREMENTS 


A  program  called  MINCAL  was  written  to  analyze  plasma 
calorimeter  data  aquired  by  NRL's  PDP-11 /10-based  data  acquisition 
system.  MINCAL  is  a  flexible,  user-oriented  program  with  the  following 
features: 

•  Parameters  are  stored  permanently  on  disk. 

•  The  program  provides  options  to  list  and/or  change  parameters. 

•  Energy  and  velocity  data  reduction  is  possible. 

•  Plots  are  done  on  a  line  printer  using  a  line  printer  graphics 
(LPG)  routine. 

•  Program  stores  all  system  constants  (amplifier  gain,  etc.)  and 
automatically  reduces  data  to  absolute  values. 

•  Program  reduces  plasma  calorimeter  data,  interpolates  between  data 
points,  and  plots  energy  and  velocity  versus  angle. 

•  Program  outputs  useful  calculations  for  user:  9 1/2  (angle 

containing  half  the  mass,  momentum  versus  angle,  others  as 
specified  by  the  user). 

•  Program  is  designed  to  easily  incorporate  new  data  reduction 
routines. 

•  Program  is  implemented  on  the  DEC-10  and  the  PDP-11/10. 

Program  Listing  is  provided  in  Appendix  D. 
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TRACER  METHOD 


The  tracer  method,  coinvented  by  an  MIC  employee  (3.  Grun),  has 
been  used  to  study  the  ablation  plasma  flow  from  planar  targets  irradiated 
by  uniform  and  perturbed  laser  beams.  The  experiment: 

•  Discovered  that  the  ablating  plasma  flows  away  from  the  target  as 
if  it  were  an  incompressible  fluid. 

•  Explained  the  angular  distribution  of  ablation  plasma  momentum, 
velocity,  and  mass  as  measured  by  diagnostics  placed  far  from  the 
target. 

•  Provided  qualitative  evidence  for  velocity  smoothing  caused  by 
thermal  conduction  in  the  overdense  ablation  plasma. 

Initial  experiments  to  spectrally  and  spacially  resolve  the 
emission  from  a  planar  target  irradiated  by  a  nonuniform  laser  beam  have 
been  performed.  From  such  measurements,  the  temperature  and  density  of 
the  plasma  along  the  target  surface  may  be  resolved  and,  thus,  the  "cloudy 
day  effect"  may  be  quantitatively  studied. 

MRC  has  applied  the  tracer  method  to  the  study  of 
hydrodynamically  unstable  targets.  This  aspect  of  our  work  will  be 
discussed  below. 

MRC  has  also  helped  to  use  tracers  to  make  local  spectroscopic 
measurements  within  laser-produced  plasmas.  This  application  of  the 
tracer  method  constitutes  a  major  advance  in  plasma  spectroscopy,  since  it 
potentially  eliminates  problems  of  interpreting  spectra  originating  from 
regions  of  varying  density  and  temperature  in  an  inhomogeneous  plasma  and 
reduces  complications  due  to  plasma  opacity. 

Appendices  A,  E,  and  F  provide  more  detail. 


5 


TEKTRONIX  DATA  ACQUISITION  SYSTEM  AND  SOFTWARE 


MRC  recommended  the  purchase  of  a  data  acquisition  system 
consisting  of  a  Tektronix  4052  computer,  4907  disk  drive,  7612D 
digitizers,  4631  hard  copy  unit,  and  an  HP- IB  37203A  extender.  The 
installation  and  operation  of  this  system  was  overseen  by  MRC.  This 
system  is  being  used  to  acquire  ion  time-of- flight  data  and  plasma 
calorimeter  data. 

The  data  acquisition  software  for  this  system  was  custom  built 
by  Tektronix  to  MRC  specifications.  Subsequently,  the  Tektronix  program 
was  modified  by  MRC  to  help  on-line  analysis  of  the  data.  A  listing  of 
this  software  is  in  Appendix  G. 

Other  software  written  to  aid  analysis  of  the  data  includes  two 
programs  to  reduce  focal  spot  size  data  (Appendix  H)  and  a  program  used  to 
calculate  the  speed  and  the  number  of  classical  Rayleigh-Taylor  e- foldings 
of  an  ablatively-accelerated  target  (Appendix  I). 

TARGETS  AND  TRACERS  FOR  RAYLEIGH-TAYLOR  EXPERIMENTS 

MRC  suggested  that  tracers  embedded  in  perturbed  targets  be  used 
to  study  the  ablation  plasma  flow  characteristic  of  hydrodynamically 
unstable  targets.  In  our  opinion  these  targets  should  have  the  following 
characteristics  -  some  of  which  are  necessary,  others  desirable: 

•  Targets  should  have  uniform  density. 

•  Periodic  structure  should  have  a  wavelength  10 nn  and  perturbation 
magnitude  of  0-2. 

•  Tracer  emission  should  be  significantly  larger  than  target  emission 
in  the  energy  band  being  measured. 
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•  Target  perturbation  period  should  be  larger  than  the  tracer  period 
for  acceptable  resolution. 

•  Tracers  should  not  perturb  target  areal  density. 

•  Tracers  should  not  perturb  the  dynamics  of  the  ablation  plasma  (for 

example  pv2  should  be  the  same  for  the  target  and  the  tracer). 

•  It  is  desirable  that  the  tracer  period  match  the  perturbation 
period  for  some  experiments. 

•  It  is  desirable  that  the  target  manufacturing  methods  match  the 

capabilities  of  the  target  fabrication  group  at  NRL.  In  particular 
polystyrene  targets  are  desirable  since  this  group  has  had 

extensive  experience  making  them. 

MRC  has  evaluated  over  40  materials  for  their  suitability  and 
has  examined  various  manufacturing  methods.  We  found  that  carbon  targets 
with  silicon  or  aluminum  tracers  met  all  our  criteria  except  the  last.  To 
extend  the  capabilities  of  the  NRL  target  fabrication  group  MRC  located  a 
used,  ion-mill  machine  which  was  purchased  by  NRL  permitting  the 

manufacture  of  the  above  targets. 

EXPERIMENTAL  METHODS  FOR  RAYLEIGH- TAYLOR  MEASUREMENTS 

In  addition  to  the  tracer  method  as  applied  to  the 

Rayleigh-Taylor  problem,  MRC  invented  a  magnesium-layer  method  that  may 
provide  a  measurement  of  the  Rayleigh-Taylor  growth  rate. 

The  tracer  and  magnesium-layer  methods  and  initial  experiments 
are  described  in  Appendix  3. 
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LONG  SCALELENGTH  EXPERIMENTS 


MRC  has  participated  in  long  scalelength  experiments  in  support 
of  the  NRL  effort.  Details  are  provided  in  Appendix  K. 

X-RAY  DETECTORS 

The  existing  NRL  x-ray  data  acguisition  system  has  been  upgraded 
and  operated  by  MRC.  In  particular: 

•  The  detectors  and  their  associated  electronics  have  been 

reconfigured  (with  the  addition  of  several  new  components  and 
electrical  shielding)  in  order  to  eliminate  a  pre-existing 

electrical  noise  problem. 

•  The  x-ray  detectors  have  been  operated  in  conjunction  with  the  NRL 
long-scalelength  experiments. 

•  MRC  has  maintained  and  upgraded  the  NRL's  PDP-11/10  computer  system 
and  software.  This  has  included  the  purchase  and  installation  of  a 
new  RT-11  operating  system,  as  well  as  streamlining  the  computer 
interface  for  more  reliable  operation. 

•  A  program  titled  PADCAL  has  been  written  to  facilitate  calibration 

of  the  x-ray  detector  system.  This  program  has  significantly 

reduced  both  the  time  and  effort  previously  required  in  the 
performance  of  this  task  and  allows  calibrations  to  be  performed 
much  more  frequently  in  the  course  of  an  experiment. 

•  MRC  has  investigated  available  options  towards  the  possibility  of 
integrating  the  Tektronix  4052  and  PDP-11/10  data  acquisition 
systems  into  a  single,  coherent  system.  These  options  will  be 
considered  by  NRL  personnel  for  future  implementation. 
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Ablative  acceleration  of  planar  targets  to  high  velocities 

J.  Grun,*1  S.P.  Obenschain,  B.  H.  Ripin.R.R.  Whitlock,  E.  A.  McLean,  J.  Gardner,  M.J.  Herbst, 
andJ.A.  Stamper 

Naval  Research  Laboratory.  Washington.  D.  C.  20375 
(Received  8  January  1982;  accepted  20  October  1982) 

Laser  irradiated  targets  are  ablatively  accelerated  to  velocities  near  those  required  for  fusion  pellet 
implosions  while  remaining  relatively  cool  and  uniform.  The  target  velocities  and  velocity  profiles  are 
measured  using  a  double-foil  method,  which  is  described  in  detail.  Also,  the  ablation  plasma  flow  from  the 
target  surface  is  spatially  resolved,  and  the  scalings  with  absorbed  irradiance  of  the  ablation  pressure, 
ablation  velocity,  and  mass  ablation  rate  are  determined.  Results  are  compared  with  hydrodynamic  code 
calculations. 


I.  INTRODUCTION 

One  approach  to  inertia!  confinement  fusion  involves 
the  use  of  a  multi-nanosecond,  moderate-irradiance  laser 
pulse  to  implode  a  hollow  pellet  containing  DT  fuel.1  If  this 
concept  is  to  succeed,  the  imploding  pellet  shell  must  be  effi¬ 
ciently  and  uniformly  accelerated  to  a  high  velocity,  and  the 
DT  fuel  must  be  kept  cool.  Only  then  will  the  pellet  compress 
to  densities  sufficient  for  a  high  gain  thermonuclear  bum. 

In  this  paper,  we  describe  experiments  which  utilize  ab¬ 
latively  accelerated  planar  targets  to  model  large  pellet  shells 
early  in  the  implosion  phase.  We  have  made  detailed  mea¬ 
surements  of  the  motion  of  such  targets2,3  4  and  of  the  abla¬ 
tion  plasma  that  accelerates  them. 2,5,6  These  measurements 
include:  emitted  light  studies;  visible  light  shadowgraphy 
and  x-ray  shadowgraphy  on  a  novel  target-impact  foil  (dou¬ 
ble-foil)  configuration  which  is  used  to  determine  the  mean 
velocities  and  velocity  profiles  of  the  dense  part  of  the  accel¬ 
erated  target,  to  study  irradiance  symmetrization,  and  to 
study  double-foil  interactions;  diagnosis  of  the  impact  foil 
using  a  triple-foil  configuration;  spatially  resolved  studies  of 
the  ablation  plasma  flow  from  uniformly  and  nonuniformly 
irradiated  targets;  as  well  as  measurements  of  the  ablation 
plasma  pressure,  velocity,  and  mass  ablation  rate.  The  target 
velocity,  velocity  profiles,  and  the  temperature  of  the  target 
rear  have  been  studied  as  a  function  of  target  thickness,  irra¬ 
diance,  and  irradiance  uniformity.  We  have  accelerated  rela¬ 
tively  cool  and  uniform  targets  to  velocities  of  160  km/sec — 
a  velocity  which  is  above  the  minimum  thought  to  be  re¬ 
quired  for  laser  fusion.'  Our  experimental  results  are  com¬ 
pared  to  hydrodynamic  code  calculations. 

The  experimental  arrangement  is  described  in  Sec.  II. 
Measurements  on  the  accelerated  target  are  described  in  Sec. 
Ill;  and  measurements  of  the  ablation  parameters  are  de¬ 
scribed  in  Sec.  IV.  Conclusions  are  offered  in  Sec.  V. 

II.  EXPERIMENTAL  ARRANGEMENT 

These  experiments  are  performed  using  the  NRL 
Pharos  II  Ndrglass  laser  U  =  1.05  (tm\  which  has  been  de¬ 
scribed  elsewhere.7  The  laser  beam  is  typically  focused  to  a 
large  diameter  |~  1  mmi  spot  on  the  target  surface  by  a  1.2 


*  Prevent  address:  Mission  Research  Corporation.  Alexandria.  Virginia. 


m,  f/6  aspheric  lens.  Such  large  diameter,  quasinear  field 
laser  spots  are  used  whenever  possible  to  maximize  irra¬ 
diance  uniformity  and  to  minimize  focal-spot-edge  ef¬ 
fects.2,5,8  Irradiance  profiles,  such  as  those  shown  in  Fig.  1, 
are  measured  on  each  discharge  in  an  equivalent  target 
plane.9  These  profiles  are  fairly  uniform  with  small  and  large 
scale  peak-to-valley  variations  of  about  +  35%.  Unless 
otherwise  stated,  irradiance  is  varied  by  adjusting  the  laser 
energy  and  not  by  changing  the  focal  spot  size.  The  laser 
pulse  duration  is  3-5  nsec  FWHM. 


FIG.  I  Azimuthally  averaged  radial  intensity  profiles  at  the  target  plane. 
Focal-spot  radius  is  varied  by  moving  the  focusing  lens  Arrows  indicate 
radii  within  which  50%  or  90%  of  the  energy  is  contained  Inserts  show  ihe 
temporal  behavior  of  the  laser  pulse  and  constant  irradiance  contours  of  a  I  • 
mm  diameter  laser  spot. 
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Figure  2  shows  the  experimental  arrangement.  The  tar¬ 
get  normal  is  tilted  6°  from  the  laser  axis  to  make  it  experi¬ 
mentally  accessible.  Surrounding  the  target  are  arrays  of 
plasma  calorimeters,  in-situ  calibrated  ballistic  pendula,2 
and  time-of-flight  ion  collectors  which  monitor  the  ablation 
plasma  and  target  energies,  momenta,  and  velocities,  respec¬ 
tively.  (In  this  long  pulse  regime  ion  collector  f  aces  are  nar¬ 
row  and  single  peaked,  making  them  ideal  for  time-of-flight 
measurements.)  From  the  angular  distributions  of  these 
quantities  we  infer  the  ablation  pressures,  mass  ablation 
rates,  and  ablation  velocities.  Although  any  two  of  the  mea¬ 
sured  quantities  are  sufficient  to  calculate  the  third,  we  mea¬ 
sure  all  three  with  independent  diagnostics  to  verify  the  con¬ 
sistency  and  validity  of  the  results.  Spatially  resolved 
pictures  of  the  ablation  plasma  flow  are  also  obtained  using 
tracer  elements  deposited  in  the  target  surface. 

Temporally  and  spatially  resolved  measurements  of  the 
target  motion  are  made  using  a  double-foil  technique.  This 
method  consists  of  placing  a  thin  foil  (impact  foil]  behind  the 
laser  accelerated  target  and  observing  the  reaction  of  the 
impact  foil  to  its  collision  with  the  target  |Fig.  3).  We  observe 
the  two-foil  interaction  in  three  ways:  ( 1 1  from  the  side  with 
optical  shadowgraphy  or  optical  streaking  (Fig.  2),  (2]  from 


FIG  J  Schematic  of  the  double-foil  concept  Behavior  of  (he  laser  acceler¬ 
ated  target  is  inferred  from  the  response  of  the  impact  foil  to  the  collison 


the  side  using  an  x-ray  backlighter,  and  (3)  from  the  rear 
recording  light  emitted  by  the  impact  foil  when  it  collides 
with  the  target  (Fig.  2).  The  details  of  particular  double-foil 
configurations  depend  on  the  intent  of  the  specific  experi¬ 
ment  and  will  be  discussed  in  later  sections. 

Temporally  resolved  temperature  at  the  rear  surface  of 
*he  accelerated  target  is  determined  by  comparing  the  abso¬ 
lute  visible  continuum  emission  from  the  target’s  rear  to 
black  body  emission. 10  The  continuum  intensity  is  measured 
at  two  different  wavelengths  using  3/4-m  spectrographs 
coupled  to  1-nsec  risetime  photomultipliers. 

III.  TARGET  DYNAMICS 

In  earlier  studies  of  ablative  acceleration  we  inferred  the 
target  velocity  from  the  velocity  of  its  debris  measured  with 
time-of-flight  ion  collectors,  and  from  the  velocity  of  its  opti¬ 
cal  shadow.8  Each  of  these  methods  however,  has  limita¬ 
tions.  Ion  collectors,  for  example,  measure  properties  of  the 
target  very  far  from  its  acceleration  region,  while  optical 
shadowgraphs  can  not  distinguish  the  high  density  part  of 
the  target  from  any  low  density  plasma  that  may  precede  it. 
In  the  present  work  we  use  a  double-foil  technique"  that 
discriminates  against  such  a  low  density  plasma,  and  permits 
temporally  and  spatially  resolv  ed  diagnosis  of  the  dense  part 
of  the  target.  Double  foils  may  also  be  used  to  investigate 
hypervelocity  impact  phenomena  such  as  the  generation  of 
ultra  high  pressures'2  or  to  model  double-shell  fusion  pellets. 

With  the  double-foil  method,  the  target  velocity  is  de¬ 
termined  from  its  flight  time  to  the  impact  foil,  placed  a 
known  distance  away.  The  flight  time,  in  turn,  is  determined 
by  noting  the  moment  at  which  the  impact  foil  reacts  to  its 
collision  with  the  target.  Such  reaction  includes  any  of  the 
observable  phenomena  normally  associated  with  colli¬ 
sions — for  example,  spall  or  fluff  from  the  rear  of  the  impact 
foil,  movement  or  deformation  of  the  impact  foil,  or  visible 
light  emission  from  its  rear  surface.  The  target  position  as  a 


FIG  4  Optical  shadowgraphs  of  a  single  la.b.c)  and  a  double-foil  Id.e.fi 
target  irradiated  at  5  >- 10' ;  W/cm;  The  limes  notes  are  relative  to  the  peak 
of  the  laser  pulse  Photographs  b.  c,  and  e,  f  are  taken  during  the  same  dis¬ 
charge. 
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function  of  time  is  measured  by  varying  the  initial  target- 
impact  foil  spacing,  thereby  controlling  the  instant  of  colli¬ 
sion.  In  this  manner  we  can  diagnose  the  target  during  as 
well  as  after  its  acceleration  phase.  Also,  the  visible  light 
emission  from  the  rear  of  the  impact  foil  gives  spatial  resolu¬ 
tion  along  the  target  surface  so  that  the  relationship  between 
velocity  and  irradiance  uniformity  may  be  studied. 

The  choice  of  the  target  and  impact-foil  material,  their 
dimensions,  and  the  spacings  between  the  two  foils  depend 
on  the  purpose  of  the  experiment.  To  infer  target  velocities, 
for  example,  it  is  desirable  that  the  impact  foil  be  thin  enough 
so  that  the  sound-transit-time  through  it  is  very  short  com¬ 
pared  to  the  flight  time  of  the  target.  In  this  case,  the  sound- 
transit-time  may  be  ignored,  and  the  collision  time  of  the 
target  with  the  front  surface  of  the  impact  foil  is  directly 
related  to  the  reaction  of  the  impact  foil’s  rear  surface. 

The  ability  of  the  double-foil  method  to  discriminate 
against  a  low  density  plasma  preceding  the  dense  part  of  the 
target  is  demonstrated  in  Fig.  4.  Here,  we  compare  shadows 
cast  by  a  single  and  a  double-foil  target  that  are  temporally 
resolved  by  multiple-frame  shadowgraphy  using  a  0.5  nsec, 
5270-A  laser  probe  flash.  Both  types  of  targets  are  photo¬ 
graphed  before,  as  well  as  1.3  nsec,  and  3.8  nsec  after  the 
peak  of  the  laser  pulse.  Figures  4(a),  (b),  and  (c)  show  shadow¬ 
graphs  of  the  single  target  which  is  a  7-^im  thick  CH  foil.  The 
shadow  of  this  target  moves  with  a  velocity  of  1.6  X  107  cm/ 
sec  and  traverses  a  distance  of  about  400  /i m  in  2.5  nsec. 
Figures  4(d),  (e),  and  (f)  show  a  corresponding  sequence  for 
the  double  foil  which  is  composed  of  a  10-/rm  thick  CH  tar¬ 
get  and  a  7-/rm  thick  aluminum  impact  foil.  Note  that  al¬ 
though  the  target  shadow  has  already  reached  the  impact 
foil  in  Fig.  4(e),  the  impact  foil  itself  does  not  react  at  this 
time.  In  fact,  the  first  indication  of  the  impact  foil’s  reaction 
occurs  in  Fig.  4(f) — more  than  2.5  nsec  after  the  optical  shad¬ 
ow  of  the  accelerated  target  reached  the  impact  foil.  Since 
the  shock-transit-iime  through  the  impact  foil  is  less  than  1 .4 
nsec  (sound  speed  in  cold  Al  is  5  X  105  cm/sec),  much  of  this 
delay  is  caused  by  the  material  in  the  leading  edge  of  the 
target  shadow  exerting  insufficient  pressure  for  the  impact 
foil  to  react.  Because  the  pressure  exerted  by  a  moving  mass 
is  proportional  to  its  density,  we  conclude  that  the  leading 
edge  of  the  accelerated  target  shadow  is  composed  of  low 
density  material.  This  low  density  material  obscures  the 
dense  part  of  the  accelerated  target  from  optical  shadow- 
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graphy  and  interferometry  diagnostics,  but  is  discriminated 
against  with  the  double-foil  method. 

Another  way  to  view  the  delayed  reaction  of  the  impact 
foil  is  presented  in  Fig.  5,  which  shows  the  streak  record  of  a 
double-foil  collision  backlit  by  a  10-nsec,  5270-A  laser 
probe.  This  one-dimensional,  but  temporally  continuous, 
shadowgraph  shows  the  low  density  plasma  advancing  to¬ 
ward  the  impact  foil,  the  low  density  plasma-impact  foil  col¬ 
lision,  and  the  reaction  of  the  impact  foil  when  the  dense  part 
of  the  accelerated  target  finally  strikes.  The  sudden  reaction 
time  (~0.5  nsec)  of  the  impact  foil  suggests  a  localized,  dense 
accelerated  target.  The  time  of  the  dense  target-impact  foil 
collision  is,  therefore,  a  good  marker  for  the  target  time-of- 
flight. 

We  determine  a  target  velocity  by  varying  the  initial 
target-impact  foil  separation  and  accumulating  distance 
traveled  versus  time-of-travel  data  for  many  discharges. 
Such  data  for  1-fxm  thick  and  10-//m  thick  CH  targets  irra¬ 
diated  at  6x  1012  W/cm2  are  shown  in  Fig.  6.  For  most  of 
these  cases,  the  target-impact  foil  separation  is  large  enough 
so  that  the  collision  occurs  after  the  ablative  acceleration 
phase  is  nearly  over.  Therefore,  the  data  for  each  target 
thickness  lie  on  a  straight  line  whose  slope  is  the  final  target 
velocity.  Velocities  obtained  in  this  way  agree  with  predic¬ 
tions  based  upon  measurements  of  the  ablation  momentum 
(pressure),  mass  ablation  rate  (Sec.  IV)  and  momentum  con¬ 
servation,  and  with  the  predictions  of  a  hydrodynamic  code 
described  below.  However,  they  are  significantly  lower  than 
the  low  density  plasma  velocities  inferred  from  the  motion  of 
the  target’s  optical  shadow  in  the  double-foil  streak  pictures. 
Averaged  over  many  discharges,  the  ratio  of  the  optical 
shadow  velocity  to  dense  target  velocity  is  2  ±  0.8,  with  the 
spread  mostly  due  to  discharge-to-discharge  variations  in 
the  low  density  plasma  velocity.  This  ratio  is  consistent  with 
the  leading  edge  of  the  low  density  plasma  being  caused  by 
spall  or  fluff  from  the  unloading  of  a  shock  at  the  target’s 
rear  surface.” 

Like  the  target  shadow,  the  impact-foil  shadow  is  also 
cast  by  low  density  material  and  not  by  the  dense  part  of  the 
impact  foil.  Thus,  the  acceleration  of  the  impact  foil  or  the 
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FIG  5.  Streaked  optical  .hadows  of  a  target  colliding  with  an  impact  foil: 
lai-before  discharge.  and  ibv-alter  discharge.  The  target  is  7-/im  thick  CH 
and  the  impact  foil  is  7-ftm  thick  Al. 


FIG.  6.  A  graph  of  target-impact  foil  separation  versus  the  time  Iwith  re¬ 
spect  to  an  arbitrary  origin)  al  which  the  impact  foil  reacts  to  its  collision 
with  the  target.  The  slopes  of  lines  drawn  through  points  accumulated  over 
many  discharges  give  the  dense  target  velocities.  The  laser  pulse  duration  is 
4- nsec  FWHM 
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FIG.  7.  Triple-foil  experiment-  A  third  foil  diagnoses  the  second  foil.  The 
User  accelerated  target  is  J-fim  thick  C  and  the  two  impact  foils  are  7-/im 
thick  C.  respectively. 


pressure  history  in  the  collision  cannot  be  gotten  from  mea¬ 
surements  using  optical  probe  light  (Fig.  5).  This  is  demon¬ 
strated  in  Fig.  7  which  shows  that  the  leading  edge  of  the 
impact  foil  shadow  exerts  insufficient  pressure  to  cause  a 
third  foil  to  react. 

Backlighting  the  double  foils  with  x  rays  provides  a  way 
to  directly  measure  the  motion  of  the  dense  portion  of  both 
foils.  An  x-ray  shadowgraph  of  a  double  foil  photographed 
using  a  J-nsec  duration  A1  x-ray  flash  (K  shell  emission, 
mostly  at  1 .6  keV;  ( mil  Be  filter  on  camera)  as  a  backlighting 
source,  is  shown  in  Fig.  8.4  This  figure  is  the  x-ray  shadow¬ 
graph  of  a  7-/rm  thick  carbon  target  irradiated  with 
3.5x  1013  W/cm3,  and  movingat3-4x  1 06  cm/sec  before  its 
collision  with  a  l-/tm  thick  carbon  impact  foil.  The  shadows 
visible  in  Fig.  8  are  the  unaccelerated  foil,  the  accelerated 
target,  and  the  impact  foil.  Bright  regions  are  due  to  the  x- 
ray  flash  and  x-ray  emission  from  the  plasma  corona  on  the 
target’s  front  surface.  We  estimate,  from  x-ray  absorption 
data  in  cold  carbon  and  the  observed  path  lengths,  that  the 
shadows  are  caused  by  material  that  is  at  least  3%  of  solid 
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FIG  8.  The  x-ray  shadow  graph  of  a  double  foil  before  the  target-impact  foil 
collision  The  x-rjv  flash  is  I  nsec  long  and  occurs  5  nsec  after  the  peak  of  the 
laser  pulse.  Pinhole  smeanng  and  time  smearing  are  about  20 /im  each 


FIG.  9  Density  profiles  of  a  thin  and  a  thick  target  accelerated  to  high 
velocities.  The  laser  is  incident  from  the  right;  is  the  critical  electron 
density.  The  simulation  conditions  are:  (a)  1.6  X  1011  W/cm2  absorbed.  2.6- 
nsec  pulse  duration,  6.5-fjm  thick  CH  target  and  (b)  1 X  10”  W/cm2  ab¬ 
sorbed,  a- nsec  pulse  duration  and  20-/jm  thick  CH  target. 

density.  Thus,  the  figure  shows  that  the  dense,  accelerated 
target  traversed  most  of  the  distance  to  the  impact  foil — a 
distance  30  times  its  thickness — without  noticeable  breakup. 
It  is  evident  too  that  the  impact  foil  does  not  react  apprecia¬ 
bly  before  it  is  hit  by  the  high  density  part  of  the  target.  These 
observations  support  our  interpretation  of  results  from  the 
double-foil  method.  Also,  velocities  obtained  by  x-ray  back¬ 
lighting  agree  with  those  obtained  by  the  double-foil  meth¬ 
od:  velocities  obtained  with  the  two  methods,  on  different 
shots,  are  the  same  to  within  35%. 

It  is  also  interesting  to  note  that  the  accelerated  target  is 
not  “cut  out”  from  the  rest  of  the  target  foil,  but  is  connected 
to  it  with  a  near  solid  density  bridge.  The  gap  between  the 
accelerated  target  and  the  target  foil  is  thus  effectively  closed 
so  that  front  surface  plasma  cannot  flow  to  the  impact  foil 
and  disturb  our  measurements.  This  bridge,  in  addition,  pre¬ 
vents  front  surface  plasma  from  heating  the  rear  of  the  tar¬ 
get,  and  thereby  spuriously  increasing  the  rear  surface  tem¬ 
peratures  measured  in  Ref.  10  and  later  in  this  paper. 

We  have  computed  target  velocities  and  density  profiles 
using  a  one-dimensional  fluid  code  that  has  been  described 
elsewhere. 14  This  code  utilizes  a  sliding  zone  Eulerian  grid 
with  flux  corrected  transport  (FCT).  The  physics  included  in 
the  code  are:  a  single  temperature  fluid,  classical  transport 
physics  with  Spitzer  heat  conduction,  and  the  same  temporal 
irradiance  history  as  the  experiment.  Absorption  in  our 
<  1014  W/cm3  long  pulse  regime  is  assumed  to  be  by  inverse 
bremsstrahlung  in  the  underdense  region  with  any  remain¬ 
ing  energy  dumped  smoothly  into  two  zones  surrounding 
the  critical  surface.  The  only  target  preheating  considered  is 
by  shocks.  Examples  of  computed  density  profiles  for  thin, 
fast  targets  and  thick,  slow  targets  are  shown  in  Fig.  9.  Both 
the  thick  and  thin  target  in  the  figure  are  compressed  during 
the  acceleration  phase.  Afterwards  they  expand  somewhat 
and  develop  a  low  density  foot  that  precedes  the  dense  target. 
Such  a  low  density  foot  in  front  of  the  dense  target  material  is 
inferred  in  the  experiment.  The  expansion  rate  of  the  dense 
target  is  much  smaller  than  its  directed  velocity.  However, 
the  expansion  rate  and  the  length  of  the  low  density  foot  are 
probably  underestimated  in  the  calculation  because  preheat 
mechanisms  other  than  shock  heating  are  not  included. 
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In  addition  to  measuring  target  velocity  magnitude, 
measuring  the  velocity  profile  across  the  target  is  also  impor¬ 
tant.  This  is  because  nonuniformly  accelerated  pellet  shells 
will  not  achieve  high  compression  ratios  and  will  thereby 
degrade  the  pellet  gain.  Nonuniform  target  velocities  could 
arise  from  any  number  of  sources  such  as  hydrodynamic 
instabilities,  nonuniform  targets  or  nonuniform  irradiation. 
We  now  address  the  third  issue,  namely:  Are  target  velocity 
profiles  and  irradiance  profiles  directly  related,  or  will  a 
mechanism  such  as  thermal  conduction  in  the  ablation  plas¬ 
ma  smooth  irradiance  nonuniformities?  Answers  to  this 
question  will  influence  the  choice  and  design  of  lasers  for 
direct  illumination  pellet  implosions. 

To  measure  target  velocity  profiles  we  utilize  the  visible 
light  which  is  emitted  from  the  rear  surface  of  the  impact  foil 
during  its  collision  with  the  target.  The  time  of  this  emission 
corresponds  closely  to  the  time  of  the  dense  target-impact 
foil  collision  as  determined  from  double-foil  streak  shadow- 
graphy  (Fig.  10),  making  the  emission  a  good  time  marker 
from  which  the  target  velocity  may  be  derived.  Now,  if  a 
nonuniformly  accelerated  target  collides  with  an  impact  foil 
the  regions  of  the  impact  foil  struck  by  fast  portions  of  the 
target  will  emit  first.  Conversely,  the  regions  of  the  impact 
foil  struck  by  slower  portions  of  the  target  will  emit  later. 
Consequently,  target  velocity  profiles  may  be  determined 
from  differences  in  emission  time  across  the  impact  foil’s 
surface. 

The  relationship  between  a  dense  target-impact  foil  col¬ 
lision  and  the  visible  emission  from  the  impact  foil’s  rear 
surface  is  graphically  illustrated  by  Figs.  11(a)  and  11(b).  Fig¬ 
ure  1 1  (a)  shows  the  x-ray  shadowgraph  of  an  impact  foil  after 
it  has  collided  with  a  nonuniformy  accelerated  target.  Non- 
uniform  acceleration  was  purposefully  induced  by  a  laser 
pulse  whose  spatial  intensity  profile  was  distorted  with  a 
rectangular  mask  that  blocked  a  portion  of  the  laser  beam 
thus  casting  a  low-intensity  shadow  across  the  center  of  the 
focal  spot.  Under  these  conditions,  the  directly  irradiated 


FIG.  10  Fcil  separation  versus  collision  time  determined  in  two  ways:  111 
from  double-foil  streak  shadow  graphy  dashed  linei.  and  |2|  from  light  emis¬ 
sion  at  the  rear  of  the  impact  foil  Pi.  Light  emission  is  measured  at  10%  of 
the  full  scale  reading  of  the  spectrograph  photomultiplier. 
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FIG.  1 1.  Two  views  of  the  collision  between  a  nonuniformly  accelerated 
target  and  an  impact  foil:  (a)  An  x-ray  shadowgraph  taken  after  the  target 
and  the  impact  foil  collided.  The  x-ray  flash  and  pinhole  smearing  condi¬ 
tions  are  the  same  as  in  Fig.  8.  (b)  Visible  emission  from  the  rear  of  a  nonuni¬ 
formly  accelerated  target  and  from  the  rear  of  an  impact  foil  struck  by  that 
target. 


regions  of  the  target  accelerate  ahead  of  the  shadowed  region 
and  subsequently  impress  their  shape  on  the  impact  foil. 
Streaked  emission  from  the  rear  of  an  impact  foil  struck  by  a 
nonuniformly  accelerated  target  is  shown  in  Fig.  1 1(b).  In 
this  figure,  we  used  a  double-foil  configuration  in  which  the 
impact  foil  was  transparent  so  that  emission  from  the  accel¬ 
erating  target  rear  and  the  impact  foil  rear  could  be  recorded 
on  the  same  discharge.  The  temporal  history  of  this  emission 
is  as  follows:  Before  the  collision,  and  while  the  impact  foil  is 
transparent,  only  the  emission  from  the  target  is  seen.  This 
emission  emanates  initially  only  from  those  target  regions 
that  are  directly  irradiated — as  seen  on  the  left  of  Fig.  1 1(b). 
As  time  goes  on,  however,  the  emitting  region  spreads  in 
space  causing  both  the  directly  irradiated  and  shaded  target 
regions  to  appear  luminous.  Such  spread  may  be  caused  in 
part  by  the  expansion  of  the  low  density,  luminous  plasma 
that  precedes  the  dense  part  of  the  target  and,  perhaps,  in 
part  by  shock  heating  of  the  shaded  region  of  the  target. 
Whatever  the  mechanism,  when  the  low  density  portion  of 
the  target  reaches  the  impact  foil,  the  foil's  optical  opacity  is 
increased  so  that  the  target  emission  is  no  longer  visible,  and 
the  middle  of  the  streak  record  appears  dark.”  Afterwards, 
when  the  high  density  portion  of  the  target  finally  collides 
with  the  impact  foil,  the  impact  foil’s  rear  begins  to  emit 
from  two  well  defined  lobes.  These  luminous  lobes  corre- 
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spond  to  the  two  impact  foil  bulges  in  Fig.  1 1(a).  Thus,  we 
have  further  evidence  that  emission  from  the  impact-foil  rear 
is  sensitive  to  velocity  nonuniformities  in  the  dense  portion 
of  the  target.  Target  structure  which  may  be  obscured  from 
direct  view  by  a  low  density  plasma  is  unmasked  by  the  dou¬ 
ble-foil  method. 

The  speed  of  the  light  emission  front  propagating 
through  aluminum  impact  foils  has  been  directly  measured 
by  using  impact  foils  with  a  thickness  step  on  their  rear  sur¬ 
face  and  noting  the  time  difference  between  emission  from 
the  two  sides  of  the  step.  We  found  that  the  emission  front 
was  supersonic  (10*’  cm/secj  when  the  foils  collided  during 
the  laser  pulse,  i.e.,  when  the  target  was  still  accelerating;  but 
it  was  close  to  sonic  if  the  collision  occurred  after  the  accel¬ 
eration.  The  supersonic  velocities  are  consistent  with  shock 
heating  causing  the  rear-surface  light  emission.  The  lower, 
sonic  velocities  can  be  explained  by  assuming  that  de¬ 
compression  of  the  target  after  acceleration  cushions  the  col¬ 
lision  impact  preventing  the  formation  of  a  strong  shock. 

Figure  12(a)  shows  the  velocity  profiles  of  dense  targets 
accelerated  to  160  km/sec — a  velocity  above  the  lower  limit 
thought  to  be  required  for  laser  fusion.1  The  rear-surface 
brightness  temperature  of  these  targets  is  measured  to  be 
below  10  eV  throughout  the  acceleration  [Fig.  12(b)].  Their 


speed  is  consistent  with  the  predictions  of  our  hydrodynamic 
code.  The  code,  for  example,  predicts  the  measured  target 
velocities  at  an  absorbed  irradiance  of  1.6x  1011  W/cm2 
(Fig.  9a),  while  the  experiment  utilizes  an  incident  irradiance 
of  2.9  XlO13  W/cm2  which — with  80%  absorption* — is 
equivalent  to  an  absorbed  irradiance  of  2.3  X 1013  W/cm2. 
Focal-spot-edge  effects2,5’*  at  the  relatively  high  irradiance 
and  small  laser  spot  size  of  this  measurement  (600/rm  diame¬ 
ter)  probably  reduce  the  effective  irradiance  below  2.3  x  10' 3 
W/cm2  and  bring  it  closer  to  the  predicted  value. 

The  fast  moving  targets  in  Fig.  12(a)  exhibit  a  velocity 
nonuniformity  -  Vmm)/Vmm  of  only  15%,  even 
though  the  peak-to- valley  nonuniformities  of  our  laser  beam 
are  about  2  to  1  (Fig.  1).  This  smoothing  of  laser  nonunifor¬ 
mities  may  be  due  to  the  so-called  “cloudy  day  effect’”6  in 
which  laser  radiation  absorbed  near  the  critical  surface  is 
thermally  smoothed  by  the  time  it  is  transported  to  the  abla¬ 
tion  surface,  where  most  of  the  pressure  is  applied  to  the 
target.  We  did,  in  fact,  make  measurements  on  the  ablation 
plasma  which  qualitatively  support  such  a  mechanism. 
These  measurements  will  be  shown  at  the  end  of  Sec.  IV. 

W  e  have  also  conducted  quantitative  studies  of  the  rela¬ 
tionship  between  laser  nonuniformities  and  target  velocity 
nonuniformities.3  To  simplify  the  analysis  we  used  slowly 


TIME  (nsac) 


FIG.  12.  Velocity  profiles  and  positions  lal,  and  rear  surface  temperatures 
Ibi  of  targets  accelerated  to  IftO  km/sec.  The  dashed  line  ts  a  theoretical  fit  to 
the  data  l°l  based  on  the  ablation  pressure  and  mass  ablation  rate  scaling 
la»s  described  in  Sec.  IV  The  peak  of  the  laser  pulse  in  ibi  is  at  0  nsec.  The 
targets  are  6.3/im  thick  CH.  the  laser  pulse  length  is  2  6  nsec,  the  laser-spot 
diameter  is  60num.  and  the  incident  irradiance  is  2.9  <  10' '  W/cmv  Veloc¬ 
ity  profiles  are  measured  with  the  double-foil  method 
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FIG.  13.  Velocity  nonuniformities  decrease  with  irradiance:  lal.  Ibi  Streaked 
emission  from  the  rear  surface  of  an  impact  foil  struck  by  a  target  irradiated 
with  a  structured  laser  profile.  The  mean  irradiance*  are  (ai  4  x  10"  W/cm" 
andibl9x  10"  W/cmv  Spatial  resolution  is  20 /im.  (c)  Target  velocity  non- 
uniformities  graphed  as  functions  of  average  irradiance  for  different  dips  in 
the  incident  beam:  1 1 1 ^  10. 220->im  FWHM  dip;  l2l 6, 
1 40-/i m  FWHM  dip;  I  J|  /I^  =*2,  !40*m  FWHM  dip  Case  121  has  the 
laser  focal-spot  contours  shown  in  the  figure.  Finite  nsetime  of  the  impact 
foil  emission  limits  velocity  nonuniformity  resolution  to  approsimately  . 
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moving  Toils  accelerated  to  no  more  than  SO  km/sec.  Under 
these  conditions,  the  fraction  of  mass  ablated  from  these  foils 
is  small  so  the  velocity  nonuniformities  are  expected  to  scale 
almost  linearly  with  irradiance  nonuniformities  in  the  ab¬ 
sence  of  smoothing,  i.e.,  (^/^Joc 

where  the  scaling  of  pressure  with  irra¬ 
diance  derived  in  Sec.  IV  is  used.  Velocity  nonuniformities 
below  this  scaling  indicate  the  presence  of  a  smoothing 
mechanism. 

Some  velocity  nonuniformity  results  are  shown  in  Fig. 
1 3.  Figures  1 3(a)  and  (b)  show  the  emission  from  impact  foils 
struck  by  targets  irradiated  with  the  spatially  structured  la¬ 
ser  beam  shown  in  the  inset  at  two  different  average  irra- 
diances.  Note  how  the  velocity  nonuniformities  are  marked¬ 
ly  reduced  at  the  higher  irradiance.  In  Fig.  13(c)  velocity 
nonuniformities  for  several  laser  irradiance  profiles  are  plot¬ 
ted.  In  all  cases  there  is  an  increase  of  irradiance  profile 
smoothing  with  higher  average  irradiance.  At  about  10 13  W/ 
cm2, 140-//m  FWHM  dips  in  intensity  are  almost  completely 
smoothed  out  but  the  200-//m  FWHM  dips  in  intensity  are 
imprinted  on  the  target.  Calculations14  indicate  that  in  this 
irradiance  regime  the  ablation  to  critical  surface  separation 
increases  with  irradiance  as  1° 7  and  is  about  100  fim  at  1013 
W/cm2.  Therefore,  these  results  are  qualitatively  consistent 
with  a  thermal  conduction  mechanism  which  smoothes  laser 
nonuniformities  when  the  ablation  surface  to  critical  surface 
distance  is  comparable  to  the  nonuniformity  scalelength.  It 
remains  to  be  seen  whether  larger  wavelength  perturbations 
at  somewhat  higher  irradiances  are  similarly  smoothed  out. 
If  they  are,  then  the  l%-3%  ablation  pressure  uniformities 
required  for  a  high-gain  pellet 17  may  be  achievable  with  real¬ 
istic  laser  systems. 

IV.  ABLATION  PLASMA  STUDIES 

In  this  section  we  present  measurements  of  ablation 
plasma  properties  relevant  for  the  target  velocity  measure¬ 
ments  described  above.  Also,  using  a  tracer  technique  to  spa¬ 
tially  resolve  the  ablation  plasma  flow,  we  provide  evidence 


FIG  14  Distribution  ofblowolT  velocity,  energy,  and  momentum  mea¬ 
sured  with  ttme-of-llight  ion  collecton,  plasma  calorimeters,  and  ballistic 
penduia  respectively  The  target  isa  1.2-mm  diameter  CH  disk  irradiated  at 

Jy  10,:  W7cm;. 
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that  the  observed  smoothing  of  target  velocity  profiles  is  in¬ 
deed  caused  by  smoothing  in  the  ablating  plasma. 

Both  foils  and  disk  targets  are  chosen  for  these  studies. 
Each  type  of  target  has  its  own  advantage:  Wide  foils,  for 
example,  do  not  permit  ablation  plasma  to  flow  around  the 
foil  edge  and  interfere  with  rear  surface  measurements.  Disk 
targets,  on  the  other  hand,  can  be  irradiated  relatively  uni¬ 
formly  by  using  a  focal-spot  diameter  larger  than  the  disk 
diameter,  and,  since  the  laser-target  interaction  area  A  is 
fixed  by  the  finite  geometry,  average  irradiances  as  well  as 
ablation  pressures  («  momentum/area)  and  mass  ablation 
rates  (amass/area)  may  be  reliably  calculated.  Also,  any 
energy  that  escapes  the  interaction  area  is  easily  detected  and 
its  effect  on  the  results  easily  estimated.  A  small  amount  of 
ablation  plasma  energy  has  indeed  been  measured  at  the  rear 
of  disk  targets  with  plasma  calorimeters,  but  its  effect  is  min¬ 
imal  in  our  experiments.  For  example,  for  300-/im  and  600- 
fim  diameter  disks  less  than  20%  of  the  total  absorbed  (plas¬ 
ma  and  target)  energy  is  detected  at  the  rear  of  the  disk,  while 
for  1000-/4R1  and  1200 -fim  diameter  disks  less  than  5%  of 
the  absorbed  energy  flows  rearward.  At  the  irradiances  em¬ 
ployed  in  these  experiments,  this  excess  energy  is  probably 
due  to  the  flow  of  hot  thermal  plasma  around  the  disk  edges. 
Very  energetic  suprathermal  electrons  observed  in  C02  laser 
experiments'*  and  at  higher  irradiance  Nd  laser  experi¬ 
ments19  are  not  the  dominant  mechanism  here,  as  evidenced 
by  the  low  level  of  x-ray  emission  observed  above  20  keV.20 

When  a  ballistic  pendulum  measurement2  is  used,  as  it 
is  here,  special  care  must  be  taken  that  events  beyond  the 
focal-spot  periphery  do  not  influence  the  experimental  re¬ 
sults.  Such  events,  if  significant,  could  not  only  complicate 
interpretation  of  the  results  but  could  make  them  irrelevant 
to  laser  fusion,  since  a  fusion  pellet  surface  has  no  edges.  For 
many  diagnostics  the  finiteness  of  the  laser  spot  is  not  a  ma¬ 
jor  concern  since  the  phenomenon  investigated  requires  high 
energies  that  exist  only  within  the  focal  spot  or  the  diagnostic 
has  time  or  space  resolution.  But  a  pendulum  is  time  and 
space  integrating,  and  a  small  amount  of  energy  can  easily 
heat  large  amounts  of  mass  to  produce  a  large  momentum. 
Small  amounts  of  such  extraneous  energy  may  be  provided 
by  thermal  conduction  through  the  focal-spot  periphery,  by 
plasma  flow  along  the  target  surface,  or  by  radiation  from 
the  plasma  plume.  The  situation  is  further  complicated  since 
these  heating  mechanisms  may  vary  with  laser-spot  size,  en¬ 
ergy,  or  irradiance.  We  have  observed,  in  fact,  that  under 
some  circumstances  the  momentum  of  material  from  be¬ 
yond  the  focal-spot  region  is  seven  times  as  large  as  the  mo¬ 
mentum  from  within  the  focal-spot  region.2  For  these  rea¬ 
sons  we  ordinarily  use  only  disk  targets  with  the  pendulum 
measurement. 

Angular  distributions  of  the  ablation  plasma  energy 
E  (0 ),  momentum  p(0 ),  and  velocity  u{& )  for  an  isolated  disk 
target  are  shown  in  Fig.  14;  u  is  a  mean  velocity  unfolded 
from  ion-collector  traces.2  Angular  distributions  from  var¬ 
ious  disk  (0.3  to  1.2  mm  diameter)  and  wide  foil  targets  have 
similar  shapes.  All  detectors  are  in  the  plane  of  incidence  and 
cylindrical  symmetry  about  the  target-normal  is  assumed. 
From  such  angular  distributions  we  determine  the  ablation 
pressure  Pt  -  PJrA,  mass  ablation  rate  m  -  m./Ar,  and 
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ablation  velocity  ut  =  Pl/ma  where  m„  is  the  ablated  mass, 
r  is  the  FWHM  pulse  duration,  and  PL  is  the  normal  compo¬ 
nent  of  the  total  momentum  P  obtained  by  integrating  p{8 ) 
over  all  solid  angles.  The  letter  A  in  the  definition  of  9 1  and 
m  represents  either  the  surface  area  of  a  disk  target  or,  for  a 
foil  target,  the  focal-spot  area  containing  90%  of  the  laser 
energy.  In  the  latter  case,  the  focal-spot  diameters  are  chosen 
to  be  large  (~  1  mm)  so  that  extraneous  focal-spot-edge  ef¬ 
fects  are  minimized. 

Scaling  of  ablation  pressure  with  absorbed  irradiance 
9 l  <*/“'  is  shown  in  Fig.  15(a). 2,5  Momentum  (pressure)  in 
the  case  of  disk  targets  is  determined  in  two  independent 
ways:  first  directly,  using  the  pendulum  array  and  second, 
indirectly,  using  the  energy  from  plasma  calorimeters  and 
velocity  data  from  time-of-flight  detectors;  wide  foil  results 
use  the  latter  method  only.  Agreement  between  all  these  in¬ 
dependent  measurements  increases  our  confidence  in  the  re¬ 
sults.  The  scalings  in  the  absorbed  irradiance  of  the  ablation 
velocity  ul  a.  I°2,  and  the  mass  ablation  rate  ma/J‘  are 
shown  in  Fig.  1 5(b). 

The  momenta  of  the  ablation  plasma  (Fig.  1 5(a)]  and  the 
momenta  of  the  accelerated  targets,  inferred  from  the  veloc¬ 
ity  measurements  in  Sec.  Ill,  balance.  For  example,  the  tar¬ 
get  momentum  for  the  cases  in  Fig.  6  is  about  0.7  of  the 
ablation  plasma  momentum.  Similar  agreement  was  ob¬ 
tained  on  a  few  irradiations  where  the  target  velocities  and 
the  ablation  plasma  momenta  were  simultaneously  mea¬ 
sured. 


ABSORBED  IRRADIANCE  I,  (W/cmJ) 


FIG  15  Ablation  pressure  la),  mass  ablition  rate  and  ablation  velocity  lb) 
versus  absorbed  laser  irradiance  for  disks  and  wide  foil  targets.  The  data  •. 
a.  A  are  inferred  using  momenta  measured  with  ballistic  pcndula.  The  data 
C.  Z.  A  are  obtained  using  calorimeter  and  time-of-flight  measurements. 
The  laser  pulse  duration  is  4-nsec  FWHM. 
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Figure  1 5  also  shows  that  the  measured  ablation  param¬ 
eters  agree  well  with  like  quantities  computed  with  our  hy¬ 
drodynamic  code.  In  this  code,  the  time  averaged  ablation 
parameters  are  calculated  by  computing  the  total  mass  and 
momentum  moving  away  from  the  target  late  enough  in  the 
run  (~14  nsec)  so  that  little  further  mass  or  momentum  is 
added  to  the  plasma.  This  is  done  so  that  the  computed  quan¬ 
tities  are  nearly  the  same  as  the  quantities  measured  experi¬ 
mentally.  The  calculated  scalings  are:  9 1  a/°”,  oc/°24, 

and  m  <x/°”. 

Table  I  compares  our  ablation  results  to  other  pub¬ 
lished  theories.  All  of  these  theories  assume  a  steady  or  qua¬ 
sisteady  ablation  plasma  flow  which  is  described  by  conser¬ 
vation  of  mass,  energy,  and  momentum:  but  they  differ  in 
other  details.  Some  theories  (Kidder,  Caruso,  and  Grat- 
ton),21-22  which  are  in  planar  geometry,  neglect  details  of  the 
absorption  and  energy  transport  to  the  ablation  layer.  They 
also  assume  that  all  of  the  absorbed  energy  supports  the  ex¬ 
pansion  of  the  ablation  plasma,  and  that  the  underdense 
plasma  expands  with  sonic  velocity.  Some  theories  assume 
sonic  flow  at  the  critical  surface  with  (Manheimer  et  al ,)2J 
and  without  (Jarboe  et  al.)2*  thermal  conduction,  and  some 
(Ahlbom  et  al.)2s  eliminate  the  sonic  flow  hypothesis  alto¬ 
gether.  Divergent  plasma  flow  from  planar  targets  was  in¬ 
cluded  by  Puell26  who  treated  plasma  hydrodynamics  in 
planar  geometry  using  a  model  similar  to  Refs.  21  and  22, 
but  considered  plasma  expansion  when  calculating  absorp¬ 
tion  of  laser  light  and  the  temperature  at  the  sonic  radius. 
Scalings  of  ablation  parameters  with  irradiance  in  spherical 
geometry  were  determined  by  Nemchinov27  who  neglected 
thermal  transport.  Scaling  of  the  mass  ablation  rate  was  de¬ 
termined  by  Gitomer  et  al.2*  who  included  thermal  trans¬ 
port  but  assumed  that  the  sonic  radius  does  not  vary  with 
irradiance.  The  effects  of  inhibited  thermal  conductivity  in 
spherical  geometry  (at  irradiances  higher  than  ours)  upon 
ablation  was  considered  by  Max  et  al 29 

Our  results  agree  with  all  the  theories  except  that  of 
Max  et  al.  which  assumes  strongly  inhibited  thermal  trans- 


TABLE I.  Companion  of  experiment  to  theoretical  scaling  laws.  Table  en¬ 
tries  are  the  exponents  b  of y*  /)  where  y  is  an  ablation  parameter,  and  /, 
the  absorbed  irradiance.  The  theories  may  be  found  in  Refs.  21-29. 


Theory 

•1 

m 

Kidder 

0.75 

0.25 

0.50 

Caruso 

0.75 

0.25 

0.50 

Ahlbom 

0.78* 

0.22* 

0.56* 

Jarboe 

0.66* 

0.33* 

0.33* 

Manheimer 

0.66 

0.33 

0.33 

Puell 

0.78 

0.22 

0  56 

Nemchinov 

0.78 

0.22 

0.56 

Gitomer 

0.56 

Max 

0.57 

0.09 

0.48 

Experiment 

0.8 

0.2 

0.6 

•7  * 7:'4  in  Ahlbom'i  expressions  have  been  treated  as  constants.  7  is 
the  absorption  fraction. 


‘Jarboe  e<  at.  explicitly  derive  a  relation  for  u,  only.  To  get  the  other  two 
quantict  we  used  the  relations  J*,  an,  y\  and  m  a  n,  a,  implied  by  his 
theory;  n,  is  the  critical  density. 
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FIG.  16.  Experimental  setup  (a)  to  spatially  resolve  plasma  flow  (b)  near  the 
target  surface.  The  five  images  are  from  five  pinholes  in  the  camera. 


port,  and  that  of  Jarboe  et  al.  and  Manheimer  el  al.  who  fix 
the  sonic  flow  point  at  critical  density  and  do  not  let  the 
exhaust  density  vary  with  irradiance. 

Until  now,  we  have  discussed  measurements  of  ablation 
plasma  parameters  far  from  the  target  surface,  i.e.,  after  the 
parameters  reached  their  asymptotic  values.  It  is,  of  course, 
desirable  to  also  know  the  hydrodynamic  development  of 
these  quantities  as  the  plasma  flows  away  from  the  target.  In 
particular,  if  the  results  of  our  experiments  are  to  be  extrapo¬ 
lated  to  spherical  geometry,  then  the  difference  between 
plasma  flow  from  our  planar  targets  and  spherical  pellets 
must  be  taken  into  account.  Toward  this  end,  we  have  devel¬ 
oped  a  tracer  technique  that  lets  us  spatially  resolve  the  plas¬ 
ma  flow  from  a  target  surface.6  Using  this  method  we  can 
also  track  plasma  flow  from  nonuniformly  irradiated  targets 
so  that  smoothing  effects  may  be  studied. 

The  arrangement  for  this  experiment  is  shown  in  Fig. 
16(a).  An  x-ray  pinhole  camera  is  placed  to  view  edge-on  a 
polystyrene  (CH)  target  with  a  row  of  small  Al  dots  (0.5  pm 
thick,  25-/rm  diameter,  spaced  50 pm  apart)  imbedded  in  it. 
When  the  target  is  irradiated  both  the  aluminum  and  the 
polystyrene  ablate  from  its  surface.  However,  since  the  alu¬ 
minum  is  a  much  stronger  emitter  of  x  rays  in  the  spectral 
band  of  the  pinhole  camera  (>  1  keV)  than  polystyrene,  the 
presence  of  aluminum  ions  in  the  plasma  flow  is  identified  by 
the  presence  of  strong  x-ray  emission.  We  assume  that  the 
perturbation  of  the  flow  pattern  by  the  aluminum  tracer  is 
small.  This  assumption  is  reasonable  since  for  a  given  elec¬ 
tron  density  n,  in  the  ablation  plasma  the  fully  ionized  ion 
mass  density  p,  for  Al  and  CH  targets  differ  only  by  12% 
Ip,  —  n,M  /Q,  where  Q  and  M  are  the  ion  charge  and  mass, 
respectively).  Also,  the  speed  of  Al  and  CH  target  ions,  mea- 
st  red  with  time-of-flight  ion  collectors,  is  the  same. 

Figure  161b)  is  a  photograph  of  the  plasma  flow  from  a 
target  irradiated  by  a  nonperturbed  laser  beam.  The  five  im¬ 
ages  in  the  figure  result  from  a  five  pinhole  array  (5-55 /im) 
through  which  the  target  was  photographed.  It  is  evident 
from  Fig.  16!b|  that  the  ablation  plasma  flow  is  reminiscent 
of  a  steady  state  fluid  flow  from  a  circular  orifice.  Note  that 
the  plasma  flow  is  planar  near  the  target  surface,  and  that 
distinct  regions  of  the  target  map  into  distinct  solid  angles. 

In  Sec.  Ill  it  was  shown  that  the  velocity  profiles  of 
accelerated  targets  were,  under  certain  conditions,  smoother 
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FIG.  17.  Ablation  plasma  flow  from  thick  CH  targets  irradiated  by  a  struc¬ 
tured  laser  profile.  The  irradiance  conditions  are:  (a|  Mean  irradiance  of 
3x  I0,J  W/cm:,  c=r8,  220-^m  FWHMdip.  (b)  Mean  irradiance  of 

6xl0'!  W/ cm;,  ~C  140-//m  FWHM  dip.  Mean  irradiance  is 

changed  here  by  reducing  the  laser  spot  diameter. 


than  expected  given  the  irradiance  nonuniformity.  We  sur¬ 
mised  that  these  smoother  velocity  profiles  resulted  from 
ablation  pressures  which  were  themselves  smoothed  by  la¬ 
teral  heat  flow  in  the  ablation  plasma.  Such  smoothing  in  the 
ablation  plasma  is  also  observed  using  the  tracer  technique. 
Figure  17(a)  shows  spatially  resolved  plasma  flow  from  a 
thick  CH  target  irradiated  by  a  laser  beam  distorted  as  in 
Fig.  1 1.  The  average  irradiance  is  3  X  'O'2  W/cm2,  the  non¬ 
uniformity  depth  is  about  8,  and  the  nonunifor¬ 

mity  width  is  220-pm  FWHM.  Under  these  conditions  no 
significant  ablation  smoothing  is  expected  and  none  is  seen 
here.  However,  smoothing  is  expected  and  seen  in  Fig.  17(b) 
where  the  average  irradiance  is  raised  by  a  factor  of  2,  the 
nonuniformity  depth  is  reduced  to  about  6,  and  the  nonuni¬ 
formity  width  is  reduced  to  140-^m  FWHM.  Under  these 
conditions,  lateral  energy  flow  from  the  directly  irradiated 
regions  into  the  central  shaded  region  ablates  the  target  sur¬ 
face  there  and  causes  distinct  streamlines  of  plasma  to  eman¬ 
ate  from  the  shaded  region.  These  streamlines  are  com¬ 
pressed  only  slightly  by  the  surrounding  plasma  which 
indicates  that  the  pressures  in  the  shaded  and  unshaded  re¬ 
gions  are  comparable,  and  qualitatively  confirms  the  pre¬ 
vious  results.  Quantitative  measurements  of  plasma  condi¬ 
tions  in  the  shaded  region  are  the  subject  of  a  future  study. 

V.  CONCLUSION 

In  this  work  we  used  ablatively  accelerated  planar  tar¬ 
gets  to  model  the  dynamics  of  large,  hollow,  ablatively  driv¬ 
en  fusion  pellets.  The  velocities  and  velocity  profiles  of  the 
dense  regions  of  highly  accelerated  planar  targets  were  mea¬ 
sured.  Despite  two-to-one  laser  nonuniformities,  we  did  ac¬ 
celerate  thin,  relatively  uniform,  and  cool  targets  to  160  km/ 
sec.  We  have  also  measured  an  increase  in  irradiance  profile 
smoothing  with  increasing  laser  irradiance.  This  is  encour¬ 
aging  for  laser  fusion.  However,  the  behavior  of  targets  with 
multimillimeter  density  and  nonuniformity  scalelengths 
that  are  typical  of  fusion  scenarios  is  not  yet  known. 

We  also  made  simultaneous  measurements  of  the  abla¬ 
tion  pressure  „  ablation  rate  m,  and  ablation  velocity  u, . 
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In  these  measurements  we  utilized  large  laser  spots  and  disk 
targets  so  that  our  results  are  not  sensitive  to  focal-spot-edge 
effects.  We  found  that  & A  and  ul  «/°J 

where  is  the  absorbed  irradiance.  Also,  the  ablation  plas¬ 
ma  momenta  and  the  momenta  of  ablatively  accelerated  tar¬ 
gets  balance.  These  results  agree  with  a  planar  geometry  flu¬ 
id  code  that  uses  classical  transport  physics. 
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Characteristics  of  ablation  plasma  from  planar,  laser-driven  targets 

J.  Gain.  *’  R.  Decoste, w  B.  H.  Ripin,  and  J.  Gardner 

A 'aval  Research  Laboratory,  Washington  D.C.  20375 

(Received  20  March  1981;  accepted  for  publication  21  July  1981) 

The  momentum,  energy,  and  velocity  characteristics  of  plasma  ablating  from  planar  targets 
irradiated  by  long  Nd-laser  pulses  (4  ns,  <  10'4  W/cnr)  are  measured  and  the  dependence  of 
ablation  parameters  upon  absorbed  irradiance  is  determined.  Large  laser  spots  are  used  in  these 
experiments  so  that  the  results  are  not  sensitive  to  boundary  effects. 

PACS  numbers:  52.50.Jm,  52.70.Nc,  52.30.  +  r 


Considerable  effort  is  underway  to  study  the  feasibility 
of  imploding  moderate  aspect  ratio,  high-gain  fusion  pellets. 
In  this  concept,  long  duration  laser  or  ion  beam  pulses  ablate 
pellet  surface  material  whose  rapid  evaporation  creates  a 
pressure  that  implodes  the  pellet  walls.  This  pressure,  and 
the  speed  of  the  ablated  material  or  the  mass  ablation  rate, 
help  determine  the  velocity  of  the  imploding  walls,  their  rate 
of  acceleration,  and  the  hydrodynamic  efficiency  of  the  im- 
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plosion.  Accurate  measurements  of  these  ablation  param¬ 
eters  are,  therefore,  basic  to  the  calibration  of  computational 
codes  used  to  design  fusion  pellets,  and  are  basic  to  the 
search  for  an  irradiance  regime  consistent  with  a  stable  and 
efficient  implosion. 1  Measurements  of  ablation  pressures  on 
planar  targets  using  small-diameter  laser  spots  and  theoreti¬ 
cal  estimates  of  the  effective  laser-target  interaction  area 
have  been  previously  reported.2  In  this  letter  we  present  si¬ 
multaneous  measurements  of  the  ablation  pressure,  the  mass 
ablation  depth  ( <x  mass  ablation  rate),  and  the  ablation  veloc¬ 
ity.  We  also  use  large-diameter  laser  spots  and  uniformly 
irradiated  disks,  so  that  the  results  are  not  sensitive  to  focal- 
spot  edge  effects  associated  with  energy  flow  through  the 
focal-spot  periphery.  The  laser-target  interaction  area  is  ex- 
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perimentally  determined.  We  utilize  redundant  measure¬ 
ments  and  diagnostics  to  verify  the  validity  and  self-consis¬ 
tency  of  these  results. 

Configuration  of  the  experiment  is  similar  to  one  pre¬ 
viously  described.5  Polystyrene  (CH)  targets  are  irradiated 
by  a  1.05 -/jm  wavelength,  4-ns  FWHM  duration  laser  pulse 
focused  through  an//6  aspheric  lens  to  an  average  irradi- 
ance  <  10'4  \V/cmJ.  This  pulse  length  is  long  enough  so  that 
stationary  ablation  physics  dominate  the  laser  plasma  inter¬ 
action.4  Ablation  pressure  &  i=Pl  / tA ,  mass  ablation 
depth  d  /pA,  and  ablation  velocity  u1  =P1  /ma  are  de¬ 
termined  from  angular  distributions  of  the  ablation  plasma 
energy  E  (6 ),  momentum  P(6 ),  and  velocity  u(6 )  measured 
with  arrays  of  plasma  calorimeters,  ballistic  pendula,  and 
time-of-flight  ion  collectors  arranged  around  the  target.  (In 
the  long-pulse  regime  ion  collector  traces  are  narrow  and 
single  peaked,  making  them  ideal  for  time-of-flight  measure¬ 
ments.)  The  quantities  r,  A,  andp  are  the  pulse  duration,  the 
laser-target  interaction  area,  and  the  target  density,  respec¬ 
tively;  PL  is  the  normal  component  of  the  total  momentum  P 
calculated  by  integrating  the  momentum  distribution  P(0)  or 
the  distribution  IE  /u  over  all  solid  angles;  ma  is  the  ablated 
mass.  Although  any  two  of  the  measured  quantities  would 
have  been  sufficient,  we  measured  all  three  with  independent 
diagnostics  to  verify  the  self-consistency  and  validity  of  our 
results. 

The  pendula  employed  here  were  tested  and  calibrated 
in  situ.*  In  particular,  we  verified  that  pendulum  recoil  is 
caused  by  plasma  impact  only  and  that  the  pendulum  array 
balances  bluwoff  plasma  and  recoiling  target  momenta  (to 
within  30%).  We  found,  however,  that  material  reflected  or 
sputtered  from  a  pendulum  surface  affects  its  absolute  cali¬ 
bration.  To  determine  the  contribution  of  material  reflection 
or  sputter  to  the  pendulum  response  and  to  thereby  calibrate 
it  absolutely,  we  built  a  double-pendulum  device.5  It  consists 
of  a  primary  pendulum  directly  facing  the  blowoff  plasma 
and  a  secondary  pendulum  oriented  to  detect  most  of  the 
material  reflected  or  sputtered  from  the  primary  pendulum 
but  shielded  from  the  blowoff  plasma.  This  in-situ  calibra¬ 
tion  showed  that  the  pendula  in  these  experiments  are 
2.4  ±  0.6  more  sensitive  then  their  bench  calibration  would 
indicate. 

To  model  the  surface  of  large  spherical  pellets  we  utilize 
wide  foil  targets  irradiated  by  large,  1-mm  diam  laser  spots. 
The  targets  are  sufficiently  thick,  so  that  their  motion  has 
little  effect  on  the  measurements.  In  some  cases  we  also  use 
uniformly  irradiated  disks.'1  These  disks  have  two  advan¬ 
tages  over  foil  targets:  first,  the  laser-target  interaction  area 
A  that  enters  into  calculations  of  absorbed  irradiance,  abla¬ 
tion  pressure,  and  mass  ablation  depth  is  experimentally  de¬ 
termined;  second,  any  energy  that  escapes  the  interaction 
area  is  easily  measured  with  plasma  calorimeters  so  that  its 
effect  on  the  results  can  be  estimated.  A  small  loss  of  ablation 
plasma  energy  to  the  rear  of  disk  targets  does  in  fact  occur, 
but  this  has  minimal  effect  on  the  ablation  parameters.  For 
example,  for  300-  and  bUQ-^im  diam  disks  less  than  20%  of 
the  total  absorbed  energy  reached  the  rear  disk  surface, 
while  for  1000-  and  1200-^m  diam  disks  less  than  5%  of  the 
absorbed  energy  reached  the  rear  disk  surface.  At  our  irra- 


diances,  the  energy  observed  behind  the  target  is  probably 
due  to  the  flow  of  hot  thermal  plasma  around  the  disk  edges: 
Very  energetic  suprathermal  electrons1  are  not  the  dominant 
mechanism  as  evidenced  by  the  low  level  x-ray  emission  ob¬ 
served  above  20  keV." 

We  note  that  the  angular  distributions  of  energy,  veloc¬ 
ity,  and  momentum  from  various  diameter  10.3-1.2  mm| 
disks  and  wide  foils  irradiated  by  various  diameter  laser 
spots  (0.3-2  mm)  are  similar.  The  momentum  angular  distri¬ 
bution,  for  example,  is  characterized  by  a  half-cone  angle 
cos~'[Pl  /P )  =  40*  for  both  the  disk  and  wide  foil  cases. 

Figure  1(a)  shows  the  scaling  of  ablation  pressure  with 
absorbed  irradiance  to  be  d? l  /“*  for  disks  and  wide  foil 
targets.  Momentum  in  the  case  of  disk  targets  is  determined 
in  two  independent  ways:  first  directly,  using  the  pendulum 
array  (dark  markers)  and  second  indirectly,  using  energy  and 
velocity  data  (open  markers);  wide  foil  results  (X )  use  the 
latter  method  only.  Agreement  between  the  two  methods 
increases  our  confidence  in  the  results.  The  scalings  with 
absorbed  irradiance  of  the  ablation  velocity  it,  or  1° 2  and  the 
mass  ablation  depth  d  *  I  aa6  are  shown  in  Fig.  1(b). 

We  also  verified  that  these  pressures,  determined  from 
asymptotic  measurements,  are  the  same  (within  30%)  as  the 
ablation  pressures  deduced  from  measurements  of  target 
acceleration.5 

The  measured  ablation  parameters  agree  well  with  like 
quantities  calculated  using  a  one-dimensional  (planar  geom¬ 
etry),  single-temperature  hydrodynamic  code  (Fig.  1).  This 
code  is  a  sliding-zone  Eulerian  FCT  code  that  requires  no 
artificial  viscosities,  utilizes  classical  thermal  transport 


ABSORBED  IRRADIANCE  I,  IW'cmJl 


FIG  I  lai  Ablation  pressure  vs  absorbed  laser  irradiance  for  disks  and  wide 
foil  targets  The  data  denoted  by  ■  •,  and  A  are  inferred  using  momenta 
measured  with  ballistic  pendula.  The  data  Ct.  O.  and  A  are  obtained  using 
calorimeter  and  ttme-of-fltght  measurements.  The  symbol  v  denotes  theo¬ 
retical  results  ibl  Mass  ablation  depth  d  and  normal  ablation  velocity  u  vs 
absorbeJ  irradiance  /„ 
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physics  with  no  ad-hoc  inhibition,  and  uses  the  same  irradi- 
ance  profile  as  the  experiment.  Absorption  is  by  inverse 
bremsstrahlung  in  the  underdense  region  with  any  remain¬ 
ing  laser  energy  dumped  smoothly  into  two  zones  surround¬ 
ing  the  critical  surface.  To  make  a  direct  comparison  with 
our  data,  we  calculate  the  momentum  and  mass  of  all  plasma 
moving  toward  the  laser  after  these  quantities  reach  their 
asymptotic  values  (14  ns  after  the  laser  pulse  peak),  and  de¬ 
termine  the  temporally  averaged  ablation  parameters 
p  x  4  m  t  in  the  same  manner  as  the  experiment.  (Peak  calcu¬ 
lated  pressures  which  occur  near  the  ablation  surface  are 
only  20%  higher  than  the  temporally  averaged  calculated 
pressures.)  The  good  agreement  between  measured  and  cal¬ 
culated  ablation  parameters  supports  the  view  that  basic 
conservation  laws  without  any  transport  inhibition  govern 
laser-target  interaction  physics  in  this  irradiance  and  laser- 
pulse  length  regime. 

Since  the  flow  characteristics  of  ablation  plasma  from 
planar  targets  are  set  by  its  nozzlelike  expansion  from  the 
target  surface,  plasmas  from  planar  targets  are  not  truly  one 
dimensional.  The  good  agreement  between  our  experiment 
and  a  one-dimensional,  planar  geometry  code  may  be  ex¬ 
plained,  however,  if  the  final  nature  of  the  ablation  param¬ 
eters  is  determined  near  the  target  surface,  when  the  plasma 
flow  is  still  planar.  This  picture  is  born  out  by  the  code  which 
shows  that  most  of  the  ablation  plasma  acceleration  occurs 
between  the  ablation  surface  and  a  few  times  the  distance 
from  the  ablation  surface  to  the  critical  surface.  This  abla¬ 
tion  surface  to  critical  surface  separation  is  about  100 gm  at 
1  X  10' 3  W/cm:  (and  less  at  lower  irradiance),  which  is 
smaller  than  the  laser-spot  diameters  or  the  disk,  diameters  in 
the  experiment.  One-dimensional,  spherical  calculations 
with  a  large  ratio  of  pellet  radius  (0.5  cm)  to  ablation  to  criti¬ 
cal-surface  distance  gave  irradiance  scaling  laws  similar  to 
those  in  planar  geometry. 

Studies  using  planar  targets  are  subject  to  criticism  that 
phenomena  at  the  focal-spot  periphery  distort  effects  as¬ 
cribed  to  the  spot  as  a  whole,  so  that  the  results  do  not  ex¬ 
trapolate  to  spherical  geometry  pellets.  Such  phenomena 
may  include  energy  leakage  laterally  across  the  focal-spot 
edge,  which  alters  the  irradiance  distribution,  and  extrane¬ 
ous  plasma  from  target  areas  outside  the  focal  spot,  which 
contributes  to  the  measured  results.  To  check  for  these  ef¬ 
fects  we  placed  planar  targets  in  the  near  field  of  the  focusing 
lens  and  varied  the  spot  size  by  aperturing  the  laser  beam 
while  monitoring  the  ablation  velocity.  Because  of  laser  en¬ 
ergy  limitations,  this  experiment  used  an  average  irradiance 
of  1  x  10'*'  W/om\  But  even  at  this  low  irradiance,  the  abla¬ 
tion  velocity  does  vary  w  ith  spot  size  (Fig.  2).  The  variation  is 
small,  however,  if  the  laser  spot  diameter  is  sufficiently  large 
I  *  1  mm  I  That  is  why  experiments  on  foils  reported  here 
utilize  large  laser  spots,  even  though  this  lowers  the  maxi¬ 
mum  irradiance  achievable  w  ith  our  laser.  It  is  also  clear  that 
controlling  irradiance  by  changing  the  focal-spot  size,  as  is 
often  done,  may  lead  to  misleading  results. 

Two  parameters  which  may  change  the  effective  spot 
size  are  the  spatial  profile  of  the  incident  irradiance  and  the 
lateral  heat  flow  through  the  focai  spot  periphery.  Since  the 
number  of.  ,r.s  at  a  particular  velocity  depends  on  absorbed 
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FIG.  2.  Ablation  velocity  u  vs  interaction-area  diameter  for  wide  foil  tar¬ 
gets.  'Interaction  area  is  defined  as  the  area  that  contains  90^  of  the  inci¬ 
dent  laser  energy.} 


irradiance,  both  these  phenomena  may  contribute  to  the  ve¬ 
locity  variation  in  Fig.  2.  Lateral  heat  flow  especially  could 
shift  energy  from  a  higher  to  a  lower  portion  of  the  effective 
irradiance  profile  and  increase  the  fraction  of  low  velocity 
ions.  Large  spot  sizes,  with  a  relatively  uniform  and  flat- 
topped  illumination,  have  a  smaller  fraction  of  absorbed  en¬ 
ergy  escaping  through  their  edges  and,  therefore,  behave 
somewhat  independently  of  spot-size  diameter. 

In  summary,  we  studied  the  ablation  characteristics  of 
plasma  from  planar,  laser-driven  targets.  We  found  that  ab¬ 
lation  pressure,  ablation  velocity,  and  mass  ablation  depth 
scale  as  the  0.8,  0.2,  and  0.6  power  of  the  absorbed  irradi¬ 
ance.  Our  results  agree  well  with  a  one-dimensional  hydro- 
dynamic  code  that  uses  classical  transport  physics.  Pres¬ 
sures  of  10  Mbar  exist  at  about  5x  10' '  W/cm:. 
Hydrodynamically  efficient  comm  ession  of  moderate  aspect 
ratio  pellets  at  this  or  slightly  higher  irradiance  may  thus  be 
possible.  We  also  found  that  due  to  edge-effects  ablation  ve¬ 
locity  varies  with  laser-spot  size,  but  the  variations  are  small 
for  sufficiently  large  spots. 
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tance  of  N.  Nocerino,  E.  Turbyfill,  M.  Fink,  L.  Seymour, 
and  R.  McGill  is  greatly  appreciated.  This  work  was  sup¬ 
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Flash  x  radiography  of  laser-accelerated  targets 
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Flash  x  radiography  of  ablatively  accelerated  planar  foils  has  provided  quantitative 
measurements  and  qualitative  observations  regarding  several  parameters  of  critical  interest  to 
direct  illumination  laser  fusion.  A  1.05-//,  3  3-ns  driver  beam  was  focused  onto  carbon  foils  in  a 
large  (0 .7- 1  -mm  diameter!  spot  to  reduce  edge  effects.  From  images  produced  by  a  backlighting  x- 
ray  flash,  we  hav  e  measured  overall  coupling  efficiency,  smoothing  of  laser  nonuniformities, 
target  velocity,  and  ablation  pressure.  The  high  velocity  targets  maintain  a  localized,  high  density 
(  >  3%  of  solid).  In  contrast  toother  workers'  recent  measurement  of  pressure  from  x-ray  imaging, 
our  x-radiographic  results,  including  pressure,  are  in  general  agreement  with  earlier  NRL  studies. 
Our  results  have  also  provided  further  insights  into  double  foil  interactions  and  planar  target 
preheat  measurements. 

PACS  numbers:  52.50.Jm,  52.25.Ps,  52.25. Fi 


Direct  illumination  of  spherical  pellet  shells  by  laser 
light  in  a  moderate  intensity  regime  somewhat  above  10u 
W/cm3  may  permit  inertial  confinement  fusion  to  be 
achieved.1  To  that  end,  studies  at  NRL3'"  have  been  direct¬ 
ed  to  the  physics  of  the  laser-pellet  interaction  and  of  the 
efficient  acceleration  of  pellet  shells  to  fusionlike  velocities.1 
Toexperimentally  model  spherical  pelletsearly  in  the  implo¬ 
sion  phase,  we  use  planar  targets  and  large  focal  spots.  These 
experiments  have  shown  that  the  moderate  intensity,  direct 
illumination  approach  scales  favorably  in  several  areas  of 
physics  which  have  been  identified  as  critical  elements  for 
laser  fusion,3  and  that,  therefore,  further  extension  of  these 
experimental  studies  to  a  regime  closer  to  reactor  conditions 
is  warranted.  In  the  present  work,  we  used  a  flash  x-ray 
probe  to  temporally  resolve  the  position  of  the  dense  region 
of  accelerated  foils;  from  this,  we  quantitatively  measured 
the  ablation  pressure,  the  target  velocity,  and  the  uniformity 
of  target  velocity.  We  also  gained  further  insights  into  dou¬ 
ble  foil  interactions  and  target  preheat  measurements. 

In  earlier  work,7' v  visible  probes,  operating  within  limi¬ 
tations  imposed  by  absorption  and  refraction,  have  been 
used  to  image  ablatively  accelerated  targets.  Flash  x  radiog¬ 
raphy13- Misapplied  here  to  planar  target  geometry  to  extend 
the  probing  photon  energy  by  about  10'  beyond  the  visible, 
thereby  gaining  simultaneous  access  to  the  front  and  rear,  as 
well  as  the  dense  middle,  of  accelerated,  high  velocity  tar¬ 
gets.  The  experimental  setup  is  diagrammed  in  Fig.  1.  A 
pinhole  camera  using  photographic  film  recorded  the  image 
of  an  x-ray  source,  an  aluminum  plasma  produced  by  a  20-J, 
0.6-iis  t.'iser  pulse."  The  object  to  be  radiographed  was 
placed  between  the  source  and  the  pinhole;  source  x  rays 
absorbed  by  the  object  result  in  a  shadow  image  of  the  object 
on  the  film.  The  objects  radiographed  were  planar  foil  tar¬ 
gets  of  pyrolytic  graphite  having  areal  densities  in  the  range 
of  1.0-1. 5  mg/cm-.  The  backlighting  x-ray  flash  was  short 
enough  that,  for  these  targets,  the  travel  distance  during  the 
flash  li  e.,  the  velocity  smcaringl  was  comparable  to  the  pin¬ 
hole  camera  resolution  labout  20//  for  a  10-//  pinhole,  mag¬ 
nification  1 .4),  thus  enabling  good  image  contrast  to  be 
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achieved.  The  carbon  |C)  foils  were  accelerated  by  NRL's 
Pharos  II  laser  operating  at  1.05//  with  3. 3-ns|  ±  1  )  driv¬ 

er  pulses  focused  toO  7-1.0-mm  diameter.  The  x-ray  source 
was  approximately  half  that  diameter.  For  some  shots,  a 
double  foil  target4■'''t•",  was  used,  in  which  an  impact  foil, 
located  behind  the  irradiated  first  foil  (see  Fig.  1),  serves  as  a 
diagnostic  of  the  motion  of  the  first  foil. 

The  source  must  be  properly  placed  and  timed  to  flash 
when  the  accelerated  foil  target  moves  into  the  view  axis.  A 
complication  arises  since  the  foil  is  also  an  x-ray  emitter; 
therefore,  we  chose  to  radiograph  the  accelerated  foil  after  it 
had  ceased  to  emit,  i.e.,  at  a  time  4-5  ns  after  the  peak  of  the 
main  (driven  pulse.  An  x  radiograph  of  a  C  foil  irradiated  at 
6.5  X  10'3  W/cm3  is  shown  in  Fig.  2(a).  The  backlighting 
source's  emission  is  imaged  on  the  left  and  the  foil’s  emission 
is  imaged  on  the  right.  Furthest  to  the  right,  x  rays  emitted 
by  the  C  blowoff  plasma  are  imaged.  Since  the  C  foil  is  wider 
than  the  laser's  focus,  only  the  central,  irradiated  portion  of 
the  foil  is  accelerated,  leaving  the  rest  of  the  foil  behind.  The 
shadow  of  this  stationary  margin  may  be  seen  in  the  image  of 
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FIO  !  Setup  for  flash  x  radiography  mot  to  scale,  angles  arc  correctly 
drawn  i  Two  I  05  ft  User  beams  were  focused  through  an // f>  lens  onto 
separate  targets  An  absorption  radiograph  of  the  backlit  object  was  record¬ 
ed  on  x-ray  film  by  a  pinhole  camera.  Single  or  double  carbon  foils  were 
employed  as  objects  Inset  scope  trace  shows  temporal  shape  of  the  mam 
laser  beam  peak  Jt  t  --  0i  and  delay  of  the  backlighter's  laser  pulse  5  n\>. 
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FIG  2  (a)  Flash  x  radiograph  of  a  single,  1.3-mg/cm2,  carbon  foil  moving  at 
3  X  10*  cm/s:  I  C  plasma  emisston.  2  Al  source  emission,  3.  shadow  of 
unaceelcrated  C  foil.  4.  shadow  of  accelerated  C  foil  lb|  Intensity  contour 
profile  of  the  main  laser  beam’s  focal  spot.  The  middle  contour  level  is 
drawn  in  thicker  lines  for  comparison  with  (a).  Incident  irradiance  was 
6.5  X  lO^W/cm*.  with  90%  of  the  laser  energy  in  a0.7-mm-diam  spot.  The 
peaks  of  the  two  laser  pulses  were  4  9  ns  apart.  Shot  81-10636. 

the  foil's  own  emission,  i.e.,  as  a  self-radiograph  [see  Fig. 
2(a)].  Since  the  accelerated  portion  of  the  C  foil  moves  se¬ 
veral  foil  thicknesses  while  radiating,  some  of  the  self-emis¬ 
sion  occurs  to  the  left  of  the  stationary  foil's  shadow.  After 
the  driver  laser  pulse  shuts  ofT,  the  accelerated  portion  ceases 
to  radiate  but  continues  to  fly  to  the  rear  (left).  The  flash  x- 
ray  source  was  positioned  just  to  the  left  of  the  stationary 
foil's  shadow,  but  extended  several  hundred  microns  further 
to  the  left  than  the  self-emission.  Within  this  several  hundred 
microns,  the  shadow  of  the  moving  foil  is  clearly  seen  and 
appears  connected  to  the  stationary  margin.  Several  interest¬ 
ing  measurements  can  be  made  from  this  representative  im¬ 
age. 

Using  the  stationary  foil's  shadow  as  a  fiducial  of  the 
original  position  and  taking  the  accelerated  foil's  shadow  as 
the  final  position,  an  average  velocity  v  =  j/rrf  may  be  com¬ 
puted.  with  Td  the  delay  of  the  x-ray  flash.  The  result  is 
3X  106  cm/s  for  this  shot,  for  an  energy  density  of  5  X  10 5  J/ 
gon  the  axis  of  the  driver  laser.  The  energy  in  the  accelerated 
foil  was  £,  =  (1/2 ]Apa  v 2,  where  A  is  the  area  and  pa  is  the 
areal  density.  The  initial  1.3  mg/cm 2  may  be  taken  as pQ, 
since  other  measurements0  indicate  that  less  than  20%  of 
the  thickness  would  ablate  by  the  time  of  observation.  This 
leads  to  an  estimate  of  total  kinetic  energy  in  the  target  of 
E,  =  (l/2)mu2  =  1.5  J.  The  incident  laser  energy  FL  was  76. 
J,  from  which  the  coupling  efficiency,  one  of  the  critical  ele¬ 
ments,2  was  7)  =  £,/£/.  x  100  =  2%;  higher  values  are 
achieved  when  a  larger  fraction  of  the  target  is  ablated. 1  This 
efficiency  combines  laser  absorption  and  hydrodynamic  ac¬ 
celeration  processes. 

Ablation  pressure  may  also  be  obtained  from  these 
data,  using  P  =  p< i/r, ,  where  rL  is  the  full  width  at  half- 
maximum  of  themain  laser.  This  leads  to  a  pressure  estimate 
of  1.4  Vbar,  at  an  incident  irradiance  of  6.5  X  10u  W/cm2. 
This  is  in  good  agreement  with  pressures  inferred  from  mea¬ 
surements  on  ablation  plasma,'*  which  showed  a  value  of 
1.7  ^  0  3  Mbar  at  an  absorbed  irradiance  of  6.5  X  1012.  The 
pressure  and  velocities  obtained  using  the  double  foil  tech¬ 


nique'  are  in  similar  agreement  with  the  value  obtained  here 

However,  our  pressure  results  do  not  agree  with  the 
recently  reported"  1.0  Mbar  (at  increased  irradiance, 

3.6  x  10,J  W/cmJ)  which  was  found  to  be  5-10  x  lower  than 
cited  measurements  on  ablation  plasma. 14  In  that  work,  a 
smaller  spot  size  (200  X  300 /rm2),  a  shorter  driver  pulse  (0.8 
ns),  and  a  longer  probe  pulse  (0.8  ns)  were  employed.  At  our 
3-ns  pulse  length,  lateral  energy  flow  within  the  blowoff  * 
extends  as  far  as  140 p  at  about  I  X  10"  W/cm2;  therefore, 
we  have  used  focal  spot  diameters  several  times  this  length  to 
avoid  a  reduction  in  pressure  due  to  lateral  flow  of  energy  to 
regions  outside  the  focal  spot.  Also  the  longer  pulse  lengths 
achieve  a  more  predominately  steady  state  ablative  accelera¬ 
tion.  However,  the  actual  source  of  the  discrepancy  noted  in 
Ref.  14  remains  unresolved. 

A  further  measurement  of  the  image  in  Fig.  2  shows 
that  the  on-axis  thickness  of  the  shadow  of  the  moving  foil  is 
about  35-70,m.  Pinhole  resolution  (20 p),  velocity  smearing 
(about  20-30  p),  and  geometric  factors  including  parallax, 
foil  rotation,  and  foil  curvature,  make  the  70  p  an  upper 
bound  for  the  axial  extent  of  this  highly  absorbing  foil.  It  is 
estimated  from  the  x-ray  attenuation  of  un-ionized  C  that 
the  observed  absorption  images  (measured  to  be  mostly  1.6- 
keV  photons,  taking  a  1/4-mm  photon  path)  represent,  as  a 
lower  bound,  a  density  of  about  3%  of  solid.  Shorter  path 
lengths,  higher  photon  energies,  and  increased  ionization 
would  raise  this  estimate.  Since  the  x-ray  range  at  1.6  keV  is 
only  8  p,  there  were  at  most  a  few  solid  bodies,  if  any,  of  this 
size  or  larger  outside  the  shadow  region  of  the  accelerated 
foil;  the  high-density  region  was  therefore  localized  to  within 
the  shadow. 

The  observation  of  a  connecting  structure  of  the  shad¬ 
ow  [Fig.  2(a)]  strengthens  the  validity  of  both  the  double  foil 
technique  and  the  measurements  of  low  preheat  obtained  for 
thin,  planar  targets.6  One  might  intuitively  expect  that  an 
accelerated  disc  would  part  from  the  unilluminated  foil, 
leaving  a  break,  through  which  plasma  could  flow  and  heat 
the  rear  of  the  foil.  Double  foil  collisions  would  then  be  in¬ 
fluenced  by  such  spurious  preheat.  Previous  investigation 
using  optical6  and  double  foil5  techniques  showed  no  evi¬ 
dence  of  such  a  plasma  flow  (see  Fig.  1  in  Ref.  6).  The  con¬ 
necting  structure  observed  in  Fig.  2(a)  evidently  isolates  the 
cool,  rear  surface  from  a  hot  plasma  flow  from  the  laser  side. 
Thus,  the  x  radiographs  explain  how  low  rear  surface  tem¬ 
peratures6  of  below  8  eV  are  maintained. 

Looking  at  the  finer  details  in  Fig.  2(a),  there  is  an  evi¬ 
dent  waviness  to  both  front  and  rear  sides  of  the  accelerated 
foil’s  shadow.  However,  neither  the  source  emission  nor  C 
foils,  of  the  type  used,  have  nonuniformities  of  the  amplitude 
and  wavelength  needed  to  explain  the  waviness.11  On  the 
other  hand,  irradiance  nonuniformities  are  known  to  exist 
and  to  influence  target  motion.  The  laser’s  focal  spot  distri¬ 
bution  was  taken  on  the  same  shot  and  shows  illumination 
nonuniformities  of  about  3:1,  with  a  spatial  wavelength  of 
about  75-100 p.  The  waviness  in  Fig.  2(a),  which  may  be 
quantified  as  about  25%  of  the  distance  travelled,  represents 
a  lateral  profile  of  nonuniformities  in  the  target's  velocity. 
The  relationship  between  velocity  nonuniformities  and  irra¬ 
diance  nonuniformities  has  been  investigated  previously  and 
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APPENDIX  D 


PROGRAM  LASER 


C 

C  THIS  PROGRArt  WRITTEN  BY  U.  MICHAEL  POLLEN,  MISSION 

L  RESEARCH  CORPORA] ION  FOR  THE  LASER  PHYSICS  BRANCH  AT 

C  I  HE  NAVAL  RESEARCH  LABORATORY. 

C 

C  =========s======cOHMON  BLOCK ========================== 

C 

INTEGER  CHAR , CMND 

COMMON/CL ASER  / 1 ADE  ( 72 ) ,  NCIIAR ,  Cl  lAK ,  CMND ,Ni ,  N2 ,  N'3 ,  N4 , 1  GOT , 
4  VAR (30,7), VAKR ( i 0 ) 

4  /CENG/EBL<9i>,EBPC?i),EFL(?l),EFP(?i),IL<20> 

4  /CVLL/V8(?i (9i),RttF(9i),RMB(91) 

OPEN  (UNIT=1,  l-TLE= 'LASERS .  DAT ' ) 

Rl  AD(1)VAR,VARR,IL 
CLOSE ( UNIT=i ) 

VARR(3)=-i 
100  WRITE (S ,110) 

110  FORMAT ( '  >'$) 

C  READS  A  STRING 

CALL  ACEIN(IGOT,IADE> 

NCHAR=i 
CALL  NXTCHR 
ISO  FORMAT ( X,3I20 ) 

If  I  CHAU. GE. 65  AND.  CHAR  .LE. 90)00  TO  480 
313  rOKMAK'  LIN  I  NIELLI  GABLE  INPUT  -  TRY  AGAIN.') 

GO  10  130 
C 

C  DECODE  COMMAND 

C 

40C  CMND=CHAR-&4 

CALL  NXTCHR 
Ni-  (I 
N2=0 
N3=3 
N4=0 

IF(NCKAR.GT.  IGOTH)GO  TO  500 
call  numbplnd 

II  (CHAR  .NE.  A 4 ) GO  TO  S00 
CALL  NXTCHR 
CALL  NUMBR ( N2 ) 

1F(CKAR.NE.44)G0  TO  500 
CALL  NXTCHR 
CALL  NUMBR (H3) 
lF(CHAR.NE.44)G0  TU  500 
CALL  NXTCHR 
CALL  NUMBR (N4) 


C 


o  o 


SI!  CALL  EXCL'TE 

IMIUOT.EQ. -?9)GL»  TO  600 
GO  TO  100 
00  CONI INUL 

UPEN(UNIT=1, FILE3 'LASERR.DAT ' ) 

UNITE  (i)yAR,VAf?R,IL 
UR  I TE  <  5 , 6i  0 ) 

.10  FORMAKX,  'PROGRAM  LASER  TERMlNAltD' ) 

CLOSE(UNIT=l) 

SI  01' 

END 

SUBROUTINE  NXTCHR 

C  GET  NEXT  CHARACTER  FROM  INPUT  STRING  (IGNORING  SPACES) 

C 

C  ===========r=====-=COMMON  EtOCK======*================ 

C 

INTEGER  CHAR , CMND 

CCMM0N/CLA3ER  / 1  A1)L  ( 72 ) ,  NCHAR ,  CHAR ,  CMND ,  Ni ,  N2 ,  Hi , N4 , 1 GOT 
4  VAR  <30,7), V ASR ( i 0 ) 

C 

L  =======«============================*===============" 

id  CHAR='IADE(  NCHAR) 

NCKAR=NC!!AR*i 
IF(CMAR.EQ.32)  GO  TO  10 
RF.  TURN 
END 

SUBROUTINE  NUMBK(NUM) 

C 

C  GET  A  NUMBER  FROM  THE  INPUT  STRING  IF  THERE  IS  QNL 
C  IF  THERE  IS  NO  NUMBER  RETURN  WITH  NUM'O 

C 

===s=s==============COMMON  BtOCK==="-=============s=== 

INTEGER  CHAR, CMND 

COMMUN/CLASCR  / 1 ADL ( 72  >  ,NCVIAR  ,CKAR  ,CMND  ,Ni  ,N2,N'i,N4,lG0T 
4  VAR(iO,7),VARR(iO) 

C 

NUM=0 
IGiGN-i 
C 

I;  CHECK  FOR  NUMBER  SIGN 

C 

IF( CHAR. N£. TOGO  TO  10 
ISIGNs- 1 

CKAR=IADE( NCHAR) 

NCHAR =NCHAR+i 


o  o  o 


LX  IT  IF  NEXT  CHAk  N01  NUMBER  Ok  I!  0L5T  OF  CHARS 

10  IF < CHAR. LT. 43. OR  CHAR.GT.57)GU  10  20 

1HNCHAR .  GT .  1GCT  i  1  )GU  TO  20 
C 

C  ADD  A  NEW  DIGIT 

C 

NUM=NUM*10+(CHAR-48) 

CI!AR=IADE(NC!IAR) 

NCIIAR=NCHARU 
GO  TO  10 

20  NUM=ISIGN*NUM 

IF ( CHAR. EG. 32 > CALL  NXTCHR 

RETURN 

END 

SUBROUTINE  ADEINiNCHAR ,1 AH  RAY ) 

DIMENSION  IARRAf (72) 

RC  AD  <S, t Q ) I ARRAY 
10  FORMA! (72R1) 

C 

C  THIS  READ  AND  DECODE  IS  PARTICULAR  TO  THE  BEC-iO 
C  IN  ORDER  TO  USE  ON  THE  11/10  OR  THE  11/03  IT  MUST 
C  BE  MODIFIED.  THE  INPUT  ARRAY  ELEMENT  MUST  BE  MASKED 

C  SUCH  THAT  THE  RIGHT  MOST  BUS  CONTAIN  THE  ASCII  CODE 

C  AND  IKE  REST  OF  THE  WGKO  IS  c ILLED  «I(H  ZEROS . 

C 

C 

C  riND  THE  END  OF  THE  LINE 

C 

DO  230  J=i,72 

i!-UAkRAY(73-J>.Nt."40)GO  TO  40 
200  CONTINUE 

J=7  3 

40  NCHAR*73-J 

RETURN 
END 

SUBROUTINE  EXCUIt 
C 

L  EXECUTE  A  COMMAND 

C 

C  ^ssssrsssssssrMsCCMMON  BLOCK=====«========  =========== 

C 

INTEGER  CHAR ,CMND 

COMMON  /CLASER/ 1 ADL  (72) ,  NCKAU ,  C'.'.AH ,  CMND ,  Ni ,  N2 ,  N3 ,  N4 ,  IGOT , 
1  yAR(20,7),yARS(i0) 

i  /CENG/EBL(?i),EDP(71),EFLi?l),EFP(?i>,lL(20) 

4  /C‘.*EL/VB(?1) ,VF(°1 ) ,RMF(?i),RMB(?i ) 


c 

L  S:==:=«S==:r=t===S!z::==S=:=r===t«r===:=:==zr====:=:= 

1  GO  TO  (100,200, 300,400,500, 600,700, 800,700, 1000,1100,1280  , 

41200, 1400, 1500,1600, 1700, 1800,1*00, 2000, 2100, 2200, 2300, 
42400, 2000, 2500), CHND 
C 

l  A 

C 

100  IOf!jTE(S,iOi)C«NMl,N2,N3,N4 

101  F0RKAT(X,315> 

GO  10  4000 

C 

C  B 

C 

200  GO  TO  5000 
C 

C  CHANGE 

C 

300  CONTINUE 

II  (Ni.EE.O)GO  TO  5000 
IF(N1.GT.30)GU  TO  5000 
IF(N2.GT.O)GO  TO  305 
GO  TO  30? 

305  IT£MP=N2 

GO  TO  335 
30?  CONTINUE 

WRITE(5,310) 

310  FORMAT!'  THE  PRESENT  VALUES  ARE:') 

URirE(5,1210) 

WRITER,  320 

320  FORMAT (2X, ,5X, ' (1 )' ,SX, ' (2) ' ,SX,' (3) ' ,?X,' (4) ' ,3X, 

4  ' (5) ' ,SX, ' (6) ' ) 

WHITE (5, 1220)  (VAK(N1,J),  J=i,7> 

WHITE (5, 330) 

330  FORMAT ( '  WHAT  VALUE  00  YOU  WISH  10  CHANGE?',*) 

REAIKS,*>I7EMP 
335  N2=llErtP*l 
URITE<5,340) 

340  FORMAT ( '  WHAT  IS  THE  NEW  VALUE?', 1) 

IRN2.LQ.2)G0  TO  3 ?7 
Rl  AD(5,*)VAR(N1,N2) 

345  WRI 1E< 5,350) VAR (Nl,i) 

350  fCRMAMX/TO  CHANGE  MOkE  VALUES  FOR  DET#',F3.0, 

4  'ENTER  VAR  #(i:-  NUT  ENTER  C):'» 

KEAD(5,#)ITEMP 
I T ( I  TEMP.EQ.O)GO  TO  3?0 
N2=ITEMP+i 
GO  TO  333 


39&  UR ITE ( 5, 395 ) 

395  I'ORMAT ( '  THE  HEW  VALUES  ARE 
WRITE(S,12iO) 

WHETS, 1220)  (VAK(Hi,J),J=i 
GO  TO  6000 

357  READ (5, 398)  VAR  (Hi  ,2) 

398  FURHAT(RS) 

GO  TO  345 
C 

C  0 

c 

400  CONTINUE 
GO  TU  5000 
C 

C  E 

C 

500  CONTINUE 

CALL  ENERGY 
GO  TU  6000 
C 

C  [- 

C 

600  CONTINUE 

GO  TO  5000 
C 

C  G 

C 

V00  CONTINUE 

GO  TO  5000 
C 

C  H 

C 

800  CONTINUE 
GO  TU  5000 
C 

C  I 
C 

900  CONTINUE 
GO  TU  5000 
C 

C  J 
C 

1000  CONTINUE 
GO  TU  5000 
C 

C  K 
C 

1100  CONTINUE 


GO  IU  5000 


r 

u 

C  LIST 
C 

i  u  00  CONI  1 NOE 

ir(Ni.n.O)GO  TO  1270 
UH ITE (5,1210) 

1210  FORMAT (X,/,'  DETt'X'TYPE'2X'ANG(TF!)'2X' ANG(PF!I)'?X'R' , 
4  3X'CAL'?X'GA1N' ) 

1F(N1.GT.Q)GC  TO  1255 
CO  1250  1=1,20 

WRITE (5, 1220)  <VAR(I,J),  J=i,7) 

1223  F0RMAT(X1FS.2,X,K5,F8.2,2X,Fa.2,5X,F5.2,SX,2(H0.3,X)) 
1250  CONTINUE 
GO  TO  1263 

1255  WRITt(S,1220)  (VAK(Ni,J>,J=i,7) 

1260  GO  TO  6000 

1270  IF  (Ni.LT.-i)GO  TO  1280 
IHN2.GT.0  .AND.  N2.LE. 10)GC  TO  1276 
WRITE  <5 , 1271 > VARR ( 1 ) 

1271  FORMAT  ( '  (i>- INCIDENT  LASER  ENERGY  <J)  ',F10.3> 

URI  l'E(5, 1272)VARR(2) 

1272  FORMAT ( '  (2)-TARGEl  ANGLE  (DEG)  ' ,K10 .3) 

URI TE(5,12?3)VAR8( 3) 

1273  FORMAT ( '  <3)-SH0T  #',F1Q.3) 

DU  127?  1=4,11) 

URITLtS, 1278)1, VARR(I) 

1273  F ORMAT ( '  ('12')', FiO. 3) 

127?  CONTINUE 
GO  TO  6000 

1276  I=N2 

UR  I  f  E ( 5 , 1 278 ) N2 , V ARR ( N2 ) 

1277  FORMAT ('  VASR('12'>  =  ',F10.3> 

GO  TO  6000 

1280  CONTINUE 

DU  1231  1=1,20 
IF (IL(I) . EQ . 0 )CO  TO  1281 
WRi7E(S, 1282)1, IL(I) 

1232  FO!:MAT(X,'PART  DLT*'I2'  SAME  AS  LIGHT  DETt'12) 

1281  CONTINUE 
GO  10  6000 

C 

L  M 

c 

1303  CONTINUE 
GO  TO  5000 

r 

C  N 


c 

1400  CONTINUE 
CO  T'J  SOUQ 
C 

C  0 

c 

1S00  CONTINUE 
GO  TU  SOOG 
C 

C  P 

C 

1609  CONTINUE 
CO  TO  5000 
C 

c  q 
c 

1790  CONTINUE 
GO  TO  5000 
C 

C  R 
C 

1890  CONTINUE 
GO  TO  5300 
C 

t  S 

c 

1700  CONTINUE 
GO  TO  5000 
C 

C  TERMINATE  PROGRAM 

C 

2000  CONTINUE 
IGO  l'=-97 
GO  TO  8000 
C 

C  U 

C 

2190  CONTINUE 
GU  T'J  5003 
C 

C  V 

C 

2200  CONTINUE 
CALL  VEL 


^200  CONTINUE 
GO  TU  SO U 0 
C 

C  X 

c 

0400  CONTINUE 
GO  TU  SCOO 

r 

w 

C  Y 
C 

0530  CONTINUE 
GO  TU  SOQO 
C 

c  : 
c 

2600  CONTINUE 
GO  TO  5001) 

SO  0  C  URITE(S,50i0> 

5010  FOKMAT  < '  INVALID  CONKAND' ) 
5020  GO  TO  6000 
6000  CONTINUE 
6010  RETURN 
END 


5  URITi:(5,  iO ) 

10  f 01' MAT  ( '  WHAT  10  THE.  SHOT  *  ?  ';> 

REAB(S,*)EF:N 

1!  (ERN.GT.  OGD  TO  23 
WlTIT2(5, 11) 

11  FOSKAT ( '  THE  SHOT  *  MOST  BE  GREATER  THAN  ZERO.') 

GO  TO  S 

20  CONTINUE 

wr;ir:<  s,igo) 

130  FORMAT ( '  WHAT  IS  THE  INCIDENT  LASER  ENERGY  (J)?  '%) 
READ(S,*)EAl!S 
V!M7E(5,90) 

90  rOL’MATCX, 'WHAT  IS  THE  TARGET  ANGLE  (DEG>?'5> 

RLAD(5,*)VARR<2) 

OPEN  ( UNIT  =1 ,  F  ILE= '  T  F"MP .  TMP ' ) 

UR ITE (5, 113) 

110  FORMAT ( '  WE  WILL  DO  FORWARD  ANGLES  FIRST.') 

CALL  ECAL(EFL,EFP,i ) 

ULTTE<5,120) 

120  FORMAT ('  NOW  WE  WILL  DO  THE  BACKWARD  ANGLES. ') 

CALL  ECAL(EBL,LBP,2) 

CLGSE(UN1 !=i) 

CALL  ETOTIEFL, TEFL) 

CALL  ETOT  < EFP , IEFP ) 

CALL  STOTfEl'L,  fE5L) 

CALL  ETOT  C EBP ,  TEE'P  > 

FFEL’=TEPP/(TEBP  :T£FP) 

PEP  =  <  TEPP+TEFP  >/EABS 
WRITE (5, 120) 

130  FORMAT (X, 'ENERGY  IN  LIGHT  "FORWARD"  EVERY  S  DBG:') 
WRITE(5,20i) (CFL(I) ,1=1,91,5) 

WK'ITE(5,131) 

131  FCRMAT(X, 'ENERGY  IN  PARTICLES  "FORWARD"  EVERY  5  DEG  '  > 
U:; I TE( 5, 231) (EFP ( I)  ,1=1 ,91,5) 


WKITE(S,132) 

rORMAKX, 'ENERGY  IN  LIGHT  "SACK*  EVERY  3  DEG  ') 

URI  IE (5,201) (EEL (I) ,1=1,91,5) 

UKITE(5,133) 

fGRMAKX,  'ENERGY  IN  PARTICLES  TACK  *  EVERY  5  DEG') 

WRITE (5 ,221) (EBP (I), 1=1,91,5) 

FORMAT  ( X ,  &F  i  0 . 3  > 

UR  1TE ' 5 ,203 ) 1 LFl , TLFP , TEBL, TE3P , FPEB , PEP 

FORMAT ('  TOTAL  FORWARD  ENERGY  IN  LIGHT*'  ,F15 .5 , '  JOULES',/, 

4  '  TOTAL  FORWARD  ENERGY  IN  PARTICLES*' ,F1C .5, '  JOULES', 

4  '  TOTAL  BACKUAKW  ENERGY  IN  LIGHT*' ,F 1C .5, '  JOULES',/, 

4  '  TOTAL  BACKHAND  ENERGY  IN  PARTICLES8' ,F10.S,'  JOULES' 

4  '  FRACTION  OF  PAS  HOLE  ENERGY  BACKWARDS* ',(-10  5, 

4  /,  '  TOTAL  PARTICLE  ENERGY/ENERGY  ADSORBED8' ,f  1(1 .5) 

MIJTir.'C  IliOCfO'j 

WHITE'S, 303) 

FOR  “AT  (XI3XF6 . 3XF£}.3XR5,FS.2,2X,Ft;.2,5X,F5 .2,5X,2'F1C .  3,  IX) ) 

FORMA !  ( X , / , '  DIU'X' INPUT' X ' CONV  INF  'X ' TYPE ' 2X ' ANG(TH) ' 2X , 

4 '  ANG  ( PHI ) '  ?X '  R '  SX '  GAL '  7X '  GAIN '  / ,  OX '  V  'EX'J'iSX'  DIG '  EX '  DEG '  ?X '  CM ' 

toy/  i  /ra'i'  \ 

*j  /  A  w ,  ^nu  / 

OPEN ( UNI T=1 , FILE*' TEMP . TtfP ' ) 

DO  300  1*1,20 

RLADd  ,30l)TJ,TJ2 

IF (TJ .EG. 0 . ) GO  TO  300 

WRITEfS, 302)1, TJ,TJ2,<VAR(I,J),J=2,?> 

CONTINUE 
DO  400  1=1,20 
READ(i,30i)TJ,TJ2 
IF ( T J . £Q . 0 . ) GO  TO  400 


WRITE'S, 302)1 

400 

CONTINUE 

CLOGC(UNIT=i ) 

30 1 

ruRf1AT(2Fi2.5 

305 

T  0SMAT(/,35X' 

4/6.1,/) 

T0RMAT(/,3SX'  THE  SHOT  *  IS: 'TS.0/,30X  ,'THL  TARGET  ANGLE  II 


c 


i 


i 

c 


200 

C 

201 


\  \  »  >  \  o  '  \  v  \  \  v  \  .  \  '  '  V  \  V  ‘  ;  \  '  \  V  \  \  V  '•  \  \  \  V  V  \  \  v  \  \  •  \  \  '  \  •.  \  \  \  \  •  \ 

/*»•>! '  I  t  !!  t  :>  t  t  t  J  f  t  I  !  t  i  t  .  f  t  }•  t  I  i  .{  J  f  t  t  S  /  f  J  /  J  .  t  f  i  t  f  i 

SJLRCUTIrtE  ZEAL  '.EL,EP .  ITY  > 

DIMENSION  EL-15, Lu<i>, £2(23. 2< 

DO  i  1=1,23 
EG(I,i)=0 . 

EC(1,2)=U 

w,RiTE(5,4Cl)(ES(J,l),E5':J,2)JJ=i,20) 

lNUri-91 

I.TJrt  IS  T-E  NUMBER  01  JNT  POINTS 
1:2112(5,103; 

FORMAT  ( '  HCv)  f1A«Y  LIGHT  DETECTOR  I  SHUTS?'}) 

RZAD(5,*)iE 


If  ( 

IE.LT 

Qxrn  tq  < 

ISO 

ChL 

l  ex; 

EL, IE, ITY' 

rr 

0  ) 

MO  t 

ttw  A 

te;s, 

401; (ES< Jj 

.1), 

r c :  t  o )  t-< 

W 

,23) 

•  :dt 

<*!.  A 

;C/C 

203  ) 

cnir 

r  u-i 

MAT  ( ' 

HIV  I 

■  wW  I.Mit  1 

PAR 

T  T  r?  -  r.G  jr  r 

1  1‘^LL  UklLL 

Tr*t*  T*jpnrc'> 

1  AlO  U)  J‘ 

READ'S .*)IZ 
L!RITEiS,23i)lE 
FORMA) ;X,'IE=', IS) 

IF ( IE.  LT  .  0 )CC  TO  353 
CALL  EX(EP,IE,1TY,£S,--1  ) 


c 

WI:l(5,45i)(ES': 

j,d,es; 

2S0 

DO  403  1=1,23 

WRlTE(i,431>ES(I 

,i;,es;i 

r 

w 

lif:ITE(5,40l)ES;  1 

,1  >  ,ES(  I 

•TGI 

FORMA l(2ri2.S) 

403 

continue 

RE  I  URN 
EnD 

SUBROUTINE  EX C E .  IE , ITY >£S , FLAG ) 

smsksssssssssCCMCN  BLOCK  =  =  ========  ===  =  =  ============ 


C 

L 


INTEGER  CHAR.CMND 

COMMON/CLASER/IADE ( 72 ) , NDfAR , C!  !AK ,  CKHD ,Ni, N2 , S3 , N4 , IGOT , 
4  ’.-(?( 33,7) ,';A?R (1 C ) 

,'CENG/EDL (71 ) , EBP ( 71 ) , EFL <  71 ) ,  EFP  ( ?  i ) ,  IE  ( 20 ) 


DIMENSION  E ( l ) , E2(29 ,2) ,ETE«P ( IS ) , LET ( 1 3 i 

HOT:  SCrtE  REASON  CAN'T  IMPLICITLY  DID  EG 

THE  ABOVE  DIMENSIONS  LIMIT  ♦  Ifll'Ul  PARTICLES  TU  10 

,;RI7Z(S,10)IC 

rCHMA! (X, ' 12=' ,13) 


o  «■->  n  *">  o  r>  o 


SUBROUTINE  ETOT  (E,TQT) 

DIMENSION  Eii) 

l  THIS  PROGRAM  USES  1 KL  INTERPOLATED  VALUES  Cl  E 

C  TO  CALCULATE  TOTAL  ENERGY  ASSUMING  P'.il  CYME  TRY  AND 

L  SYML.TR Y  ABOUT  THE  LINE  NORMAL  TO  THE  TARET 
PI=3.14i5?2fcU4 
XANG=H 1/130 

C  I.L.  XANG=i  DEGREE 
RI=3.0 
C 

INTEGRATE  I'Y  TRAPIZUIDAL  RULE 

T0i=E!?i)/2.Q 
DC  100  1=2,00 
RI=RI+XANG 

00  1 CT  =  TOT +SIN(RI )*L( I ) 

NOW  MULTIPLY  TIMES  STEP  SIZE 

TQT=lOT*XANG 
TGT=2.C*PI*TOT 

USE  MAS  BEEN  MADE  Or  T ME  FACT  THAT  SIN(0)=0  AND 
UIN(90)=1 . 0 
RETURN 
END 


SUBROUTINE  VEL 


THIS  SUBROUTINE  DSLS  THE  VELOCITY  CALCULATIONS 
srarssssssarsssssIjn^^ON  8L0CK=«====«====  ====="™" 
INTEGER  CHAR,CNND 

CONMON/CLASER/ IADE  (72) ,  NCHAR ,  CHAR ,  CMND ,  Ni ,  N2 ,  N3 ,  N4 ,  I  GOT , 
i  VAR  C  3S , 7 ) , VARS ( i 0 ) 

4  /CLNG/EBL(?1) ,EBP(?i) ,EFL(91) ,EFP(?i) ,ILi23) 
i  /CVEL/ VB  <  91 > , VL ( 9i > , RNF ( ?i ) , R  Nt ( 9 1 ) 

EQUIVALENCE  ( SEAL’S ,’MRR  ( i  > ) 

DIMENSION  RPF(9i) ,KPB«9i) 


VF  IS  AN  ARRAY  CONTAINING  THE  FORWARD  VELOCITY 
VB  IS  AN  ARRAY  CONTAINING  THE  BACKWARD  VELOCIiY 
RMF  IG  AN  ARRAY  CONTAINING  THE  FORWARD  NASS 
RNB  IS  AN  ARRAY  CONTAINING  THE  BACKWARD  NASS 
I  RNF  13  THE  TOTAL  !  CRUARD  NASS 
TRNB  IS  THE  TOTAL  BACKWARD  NASS 
ALL  ARRAYS  HAVE  SLQ  THETA  VALUES  SPACED  1  DEG 
STARTING  AT  3  (93)  AND  GOING  TO  93  (273) 


IFiVAKR(3).GT.G)GO  TO  10 
URITC(S,?> 

FORMAT  ('  YOU  NUS!  CCNPUTE  ENERGY  FIRST!!') 

GO  10  S00 

WRITE(S,ii)VARR(3) 

FORNAT ( '  UE  ARE  SUING  EXPERIMENTAL  SHOT  *',FS.3) 
0R1TE(S,?C)VARR(2) 

FORNAT ( '  THE  TARGET  ANGLE  IC'FV.2'  DEGREES' ) 
URIT£(S,iiO> 

FORNAT ( ’  WE  WILL  DO  FORWARD  ANGLES  FIRST.') 

CALL  VCAL(VF,i) 

CALL  NCAL(RNF.V) ,EFL,EFF,i,RPF> 

WRITER, 120) 

FORMAT ( '  NOW  WE  WILL  DO  THE  BACKWARD  ANGLES.') 
CALL  7CAL(VB,2) 

CALL  NCAL ( RNC , VS , EBL , EBP , 2 , RP  t ) 

CALL  MTOT(PMF.TRNF,i) 

CALL  NTOT(RNB,TRMB,2) 

CALL  VP R N.F ,  V F , TR  MF , VPF , i ) 

LfiLL  VP  ( RMB ,  VB ,  l  i! NB ,  VPB ,  2 ) 

CALL  :HII(RNF,TRNF,TitHF,i) 

CALL  T!!H(RMB,TRNB,T!IHB,2) 


A 


CALL  PT0T(RPB,TRPB,TRPPB,2) 

w;;itl<s,2C1) 

Uk  1 1 L ( S ,  30 C ) '  VP  ( I  > ,  1 =i ,  9 1 , 5 ) 

U;:ITC<S,202> 

UUi 1E( S,35C) (RMF(  1)  ,1=1 ,?1 ,S> 

UNITES, 203) 

WRITE (5, 380) (RPF (I) ,1=1 ,91,5) 
yr;ITC(S,204> 

UKHE!S,350)  (VB(I)  ,1=1,91 ,5) 

W.'?iTE<5>2u5> 

WRITE<S,303>(RMBU>,I=i,9i,S> 

WPI !  2(5,206) 

WKi  !L(S,33C)(RPB(I),I=i,9i,E) 

FORMA7(X,6Fi0.3> 

FCRMAKX,'  "FORWARD *  VELOCITY  EVERY  5  DEO : ' ) 
rCRMAl (/, '  "FORWARD"  MACS  EVERY  S  DEG:') 

FOSMAT(/,'  “FORWARD"  MOMENTUM  EVERY  5  DEG:') 

FCi<KAT(/, '  "BACKWARD"  VELOCITY  EVERY  5  DEG:') 

FGKrtAI (/, '  "BACKWARD"  MASS  EVERY  S  DEC:') 

FORMAT  (//"BACKWARD"  MOMENTUM  EVERY  5  DEO:') 

TMHB=THHB*180. 

VR I TE  <  S ,  40  8 )  VPF ,  VPB ,  TRMF ,  TRM8 ,  TSIdF  ,THHB ,1  RPF , TRPPF ,TRPB , TRPPB 
FORMAT  (X,  'FORWARD  PARALLEL  VELCCi  TY*'  ,F1D  .5/E-?  CM/S',/, 

4  X , ' BACKWARD  PARALLEL  VELOCITY*', K10.S/E-7  CM/S',/, 

4  X,' TOTAL  MASS  IN  FORWARD  ANGLES*', F12.5/E-7  0',/, 

4  X/ TOTAL  MASS  IN  BACKWARD  ANGLES*', FI2.S,'E-7  G ',/, 

4  X,'M=O.S  TOTAL  MACS  (FORWARD)  A!',F12.5,'  DEG',/, 

4  X,'M=0.5  TOTAL  MASS  (BACKWARD)  A1',F12.S,'  DIG',/, 

4  X, 'TOTAL  FORWARD  MOMENTUM*' ,F12.S , '  DYNE-SEC',/, 

4  X/TCTAL  PARALLEL  FORWARD  MUMLNTUM*' ,F12.S, '  DYNE-SEC 

4  X, 'TOTAL  BACK  MUMLNTUM* ' ,P12.S, '  DYNE-SEC',/, 

4  X,' TOTAL  BACK  PARALLEL  MOMENTUM*' ,F12.5 , '  DYNE  SEC') 

RE! URN 
END 


c 


SUBROUTINE  VCAL  (V.1TY  ) 
DiiiLNSICN  V <  i  > 


iiiO 


35  0 


C 

c 


INUM  IS  THE  NUMBER  OK  INI  POINTS 
«niTE(5,180) 

FORMAT! '  HOW  MANY  VELOCITY  INPUTS?'*) 

READ(S,*)IE 
If  (IL.LT.O)GO  TO  30 U 
CALL  VEX <  V  > IE , I TY ) 

RETURN 

END 

>»>>>)>>»>>>)»>>>»»)>>)>>»»>>»>>))>)»>>>>>>)>> 

SUBROUTINE  VEX<V,ILfITY) 

===============CCi1M0N  BLOCK  =========================== 

INTEGER  CHAR  ,CMND 

UUNMUM/CL A5ER/ 1 ADE  <?2>  >  NCHAR ,  CHAR ,  LMND , Ni , N2 , N3 ,  N4 ,  IGOT , 
4  VAR (30 APR (10) 


DIMENSION  V(i),tTEHP(i8>,DET(il}> 

C  THE  ABOVE  DIMENSIONS  LIMIT  t  INPUT  PARTICLES  fO  10 
DO  120  J=i,IE 
URI  !E(5,iS0) 

ISO  FORMAT ('  ENTER  (DEI# /VELOCITY)  '*) 

READ<  5 , *» ISET , ETErtP ( J ) 

1ULT=IBET  +20 
IF( IDLT . LE . 20 > GO  TO  S00 
IFdDET.GT. 30)00  TO  SOD 
DET(J)=FLOAT (IDLT) 

100  CONTINUE 
C 

C  CONVERT  TO  REAL  UNITS  AND  GET  ANGLE  OF  DETECTOR 
C 

DO  200  J=i,IE 
I5ET=DE7(J) 

ETEMP  ( J ) =LTLMP  ( J )  *VA!(  ( IDET ,  6 )  CVAR  ( 1DET ,  ? ) 

20  0  BE f ( J ) =VAR ( IDET , i) - VARR ( 2 ) 

C 

L  INTERPOLATE 
C 

NDEG=2 

CALL  R INT ( ETEMP , IE , DLT , V , ITY , NDEG  ) 

rrft  DC*plp».| 

w  V  ^  1  4»  I  Uf  I  * 


L 


C 

(' 

u 

C 

1 


c 

c 

c 

c 


iOO 

c 

C 

C 


C 

C 


END 


>>:>>))>>»>>>>>>>/'>>>>»>>>>»»»>»>>>>>)>»)>>)) 

SUBROUTINE  MTOT(RM,TO I >ITY) 

DIMENSION  RM(i) 


THIS  CALCULATES  TKL  TOTAL  MASS  (LI  f!!£R  FORWARD  11  Y=1 
CR  BACKWARD -ITr=2) 

PI'3. 141592654 

XANG=?I/188 

I.L.  XANG=1  DEGREE 

■i  t  -  o  r. 

»v  X“  U  .  V 


INTEGRATE  BY  TRAP  1201 DAL  RULE 
TQT=RM(?i)/2 . 0 

T.n  *  fin  7-o  n?| 

ri=ri+xang 

T0T=TCT4SIN(RI)*KM(1> 


NOW  MULTIPLY  TIMLS  STEP  SIZE 


TOI=TOT*XANG 

TOT=2.0*PI*TOr 

USE  HAS  BEEN  MADE  Of-  THE  FACT  THAT  SIN<0>=0  AND 

SIN(?0)=1 . 0 

RETURN 

END 


c>  >>>>»>>»»»;>>>»>>>»»>>»>>>)>>>>»»)>>>;>>  >>)>>>>> 

Sl-BRQU  I' I  HE  THH  ( RM  ,TRM,  THH ,  ITY ) 

C 


C  THIS  FINDS  THE  ANGLE  AT  WHICH  HALF  OF  THE  TOTAL  MASS  IS 
C  ACCUUNIED  FOR 

e 


DIMENSION  RM(i) 


P 1=3 . 141592654 
XANG=Pl/10Q . 

HTRM*0 . 5CTRM 
TCT-0 . 0 
1=1 

C  5IN(D)  =  0 

RI=0 . 

10  l=l!t 

RI=RI‘XANG 

1 0 1 =T0  T  *SIN( R I ) *IIrtU ) *XANG*2 . 0  *P I 
I;  ( lOT.GE.HTRrt)GO  TO  100 
If (I  LE.90ICO  TO  10 
I=I‘l 


c>>  »>>»»>>>»»»»>>>>>>>>>>>»»>>>»>»>>»>  >»>>»»»> 

‘JUSRUUT I NE  HCAL  ( KM ,  V ,  EL ,  EP ,  I T Y ,  W ) 

DIMENSION  RM(iJ}y(i),EL(i ),EP(i),:iP{i) 

DO  ICO  1=1,91 

l!  (V(I)  EQ.O . )GC  TO  too 

TS£Mr‘=2 . 8*  ( EP  ( I ) )/ V  ( I ) 

RP(I)=fE«P 

NM(I)=TEMP/V(I) 

100  CONTINUE 

C  WkITE(5,20D)(RM(J),J=i,9i,S) 

200  FOriHA  f<X,SF12.5) 

RETURN 

END 


< 


c 


c 

c 

c 

c 


t 

c 

c 

c 
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subroutine:  PTCT(RP .tot^tqtp , ity) 

DIMENSION  PP(1) 


\  \  \  \  \  \  \  \ 
I  I  i  I  S  /  I  f 


\  \ 
i  / 


THIS  CALCULATES  THE  TOTAL  MOM  (EITHER  FORWARD  IlY=i 
OR  DACLUARD-ITY=2)  BOTH  TOT  AND  TOT  PARALLEL. 

PI=3. 141592654 
XANG-PZ/iSO 
I.L.  XANG=i  DEGREE 
RI-0 .0 


INTEGRATE  BY  TRAPEZOIDAL  RULE 


TOTP=fl. 0 
TOT=RP(?i)/2. 0 
DO  100  1=2,90 
RI=RI :XANG 
TrMP=GIN(RI)*RP(l> 

T0TP=T0TP+CQ5'(RI)*TEMP 

TOT=TOT*TEMP 

NOW  MULTIPLY  TIMES  STEP  SIZE 

TCT--TOT*XANG 

l'OTP=TCTP*XANG 

TM=2.8*PI*TCT 

T0TP=TCTP*2.Q*PI 

USE  HAS  BEEN  MADE  Of  THE  FACT  THAI  SIN(0)=0  AND 

SIN(90)=1 . 0 

RETURN 

END 


t.J  O 


SUSP  OUT  INE  R INT  •;  El  Er.P ,  IE ,  ANG ,  E .  I T Y ,  ({BEG ) 

D1  NLNG  I  CM  ETEMN  ( i > ,  ANG  ( 1 ) ,  E  ( i ) ,  A '  i  3 ) ,  £  i '  i  0 ) ,  DU i  0  5 

EiCMP  IS  Y ARRAY 

ANU  iS  T,!E  X  ARRAY 

III  IS  THE  NO  UR  TT'J 

E  IS  THE  OUTPUT  INTERPOLATED  T  ARRAY 

II Y  TELLS  IT  FORWARD  OR  SACK 

DO  1  J=i,iC 

A  ( J )  =  0 . 0 

Bi(J)=0.B 

B2(J)=0.0 

CALL  COEF ( A , C i , EE , ETEHR , ANG , NDEG+ 1 , 1 E ) 

URl!'EiS,iS)  (A(I>  ,Bi(I>,B2(I),I=i,lil) 
rCRKAKX^FtO.3) 

IF(  ITY.EGJ .  i)ANGLL=- i 
IF(ITY.EQ.2)ANGLE=17? 

DU  300  J=i ,91 
ANGLE=ANGLE+i 
FO=i.Q 

FNE’J=ANGLE-I'1(2> 

1=1 

kINlT=A(  I  )<FOt,A'I*i)*FNi:W 
If (I.G£.NDEG)U0  TO  21 0 
I  =  M 

FH0L1)=FNEU 

F NLU=ANuLE*FNEW  Bi ( 1+ 1 ) *FNCU-B2 i I : l)?FO 
FQ=I  HOLD 

RINlT=A(  I+i  )*FM£W-s  KiWTT 

03  TO  ICC 

E(J>=RINTT 

CONTINUE 

NEC  VEL  AND  ENERGIES  ARE  NON  PHYSICAL  SET=0. 

DO  SCO  1=1,71 

I;  <L< I ) . LT . 0 .  )L< I >=0 . 0 

RETURN 

END 

>>>!>>  >>>>>>)>>>>>)  >>>>>>>>>>>)  >))»>>>>»>>>>)  >»>>>: 

SL’BROU  f INE  CCER ; A , I'i  ,B2 , Y ,X , M ,H> 

M  IS  THE  NUMBER  CF  TERMS=NDECfl 
N  IS  il'E  NUMBER  CF  INPUT  POINTS 


DIMENSION  A(i),Bl(l),B2(l),Y(i) ,X(i),FCLD(iC),FNEW(iO> 

THIS  COMPUTES  THE  INTERPOLATING  LULF  A,Bi,B2 

YFFPi)D=0.0 

XI  PRUD=0 .  C 

OLBPPD=FLUAT(N) 

DO  ICG  1=1, N 
XFrROO=XFPROD:-X(I) 


YFPF:DD=YFPRQD+Y(I) 

AC i )=Y;  PROD/FLUAKN) 

JYi  <  2 )  =XFP  ROD/FLOAT  ( N ) 

Fi=i 

CALCULATE  F2 
DO  L’OO  1=1,  N 
l'!JLD(lJ=i.O 
FNElUI)=X(I)-£i(2) 

CALCULATE  ACT), Bl(l) ,1:2(1) 

K-2 

IF (K.GT.fi)GO  TO  iOOU 
X:'PROD=0 . 0 
Yf PRCB=U .  0 
FFP2Oi)=0 . 0 

CCrif'UTE  INNER  PRODUCTS 
DO  A00  1=1, N 

YFPROD=YFPROD*Y ( I >*FNEU( I ) 

PCLYSU=FNE«(I)SI  NtWCI) 

F:TT:OD=rFPROIi:'i1OLYC!,’ 

XiPKOD=XFPP.OD+Xa>*PCLYOa 
A  (K )  =YFPROD/FFPROD 
if (K.EO.H)GO  TO  iOCQ 

CALCULATE  COEF  FO  XUTH  OKTAFIG  POLY 

81  (K  >l)=XFPROD/f  FPKOD 
D2(KM)=FFP80D/QLDPRD 
DO  SOD  1=1, N 
FHOLD=FNEW(I) 

F‘Oa)=(X(I)-SliK.4))*FNEW(I)-D2(X+l)*FOLD;i> 

FCLDfl )=FHOLD 

0LDFKD=!7PRCD 


Fluid  flow  of  ablated  plasma  from  planar  targets 
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Using  a  novel  diagnostic  technique,  the  streamlines  of  the  flow  of  plasma  from  a  laser-irradiated  target  are 
observed  for  the  first  time  Implications  for  far-field  ion  measurements  and  similarity  to  a  simple  fluid  flow 
are  noted. 


Although  laser  fusion  will  use  spherical  pellets,  many 
laser-matter  interaction  experiments  are  performed 
on  flat  targets.  If  the  results  of  planar  experiments  are 
to  be  scaled  to  spherical  geometry,  the  effect  on  plas¬ 
ma  profiles  of  differences  in  material  flow  from  the 
target  must  be  taken  into  account.  The  one -dimensional 
spherically  divergent  flow  from  a  uniformly  irradiated 
pellet  affects  plasma  profiles  at  distances  from  the 
target  which  scale  with  the  pellet  radius.  While  it  has 
long  been  recognized  that  the  more  complex  multidi¬ 
mensional  flow  from  a  finite  focal  spot  on  a  planar 
target  must  diverge  at  distances  from  the  target  sur¬ 
face  which  scale  with  the  focal  spot  radius. 1  the  first 
published  results  of  multidimensional  calculations  which 
show  the  fluid  flow  pattern  are  only  to  appear  soon.2 
In  the  present  work,  a  new  diagnostic  technique  allows 
the  first  experimental  identification  of  the  fluid  flow 
lines  from  a  planar  target. 

As  shown  tn  Fig.  1,  foil  targets  are  irradiated  in  the 
near  field  of  an  /.' 6  lens  at  a  6°  angle  of  incidence  by 
100-500  J  of  Nd  laser  (*,  =  1 . 054- urn )  radiation.  With 
an  incident  pulse  length  of  3.  5  nsec,  and  an  irradiated 
spot  diameter /-/„ (90  j  energy  content)  of  1.4  mm.  aver¬ 
age  irradiances  within  J,lQ  are  1.  5-8  TW  cm!.  The 
focal-intensity  distribution,  as  determined  by  equivalent 
focal-plane  cameras,  is  relatively  flat-topped,  with  ap¬ 
proximately  2:1  intensity  modulation  across  the  central 
region.  A  camera  with  an  array  of  pinholes  is  used  to 
obtain  time-integrated  images  of  x-ray  emission  at  90'' 
to  the  target  normal.  The  15-pm  beryllium  filters 
used  on  the  camera  allow  imaging  of  photons  with  energy 
hv  1  keV. 

A  linear  array  of  aluminum  spots,  25  pm  in  diameter 
and  10  mu  thick,  are  embedded  as  tracer  material  in 
300  ,.m  thick  polystyrene  (CH)  target  foils,  which  are 
then  aligned  so  that  die  spots  fall  along  a  diameter  of 
the  laser-irradiated  region.  Since  aluminum  line  emis¬ 
sion  is  more  int.  nse  than  the  CH  continuum  emission 
in  the  spectral  band  of  the  pinhole  camera,  stream¬ 
lines  of  material  flowing  from  the  aluminum  spots  are 
identifiable  bv  the  tracks  of  strong  x-ray  emission 
[see  Fig.  2f.i)[.  The  |>erturbation  of  the  flow  pattern 
due  to  the  addition  of  a  second  material  should  be  mini¬ 
mal  because:  <  1-  the  full  /  ionized  mass  densities  for  A1 
and  CH  differ  bv  only  12  at  a  given  electron  density 
in  the  blowoff  plasma,  and  12)  the  average  velocity  of 
ablated  ions  as  determined  by  Faraday  cups  is  the  same 
for  the  two  materials. 

The  observed  flow  pattern  resembles  that  of  the  irro- 


tational  flow  of  a  nonviscous  compressible  fluid  from  a 
eircularcr’fice.’  Near  the  target  surface,  the  observed 
streamlines  curve  away  from  the  target  normal  at  a  rate 
which  is  greatest  for  streamlines  whose  source  points 
at  the  target  surface  are  nearest  the  edge  of  the  fluid 
source.  Beyond  Z  =  r0.  where  Z  is  the  distance  from 
the  target  surface  and  r0  is  the  fluid  source  radius,  the 
curvature  is  negligible  and  material  flows  outward  in  a 
straight  path.  The  dearth  ofpublished  multidimensional 
code  results  showing  the  fluid  flow  pattern  leads  us  to 
make  an  interesting  comparison  between  the  data  and 
the  aforementioned  potential  flow.  The  latter,  at  least 
for  low  Mach  numbers  (nearly  incompressible  flow)  and 
nearly  paraxial  streamlines,  is  well  approximated  by 
the  solution  to  Laplace’s  equation  subject  to  the  boundary 
conditions  at  the  source  [see  Fig.  2(b)].  To  quantify 
the  comparison,  we  consider  the  streamline  angles  6 
relative  to  the  target  normal,  which  depend  upon  Z  and 
upon  the  source  position  r  of  the  streamline  in  the  tar¬ 
get  plane.  Fixing  Z  where  >•„  is  the  orifice  ra¬ 

dius  for  the  Laplacian  flow  and  the  observed  fluid  source 
radius  for  the  x-ray  data,  the  dependences  of  8  upon  r  r0 
are  displayed  in  Fig.  2(c);  the  similarity  of  the  two  cases  is 
striking.  Of  course,  the  model  chosen  does  not  accu¬ 
rately  represent  the  case  of  the  laser-ablated  target, 
where  fluid  emerges  at  low  velocity  from  the  solid  sur¬ 
face  and  is  subsequently  accelerated  to  sonic  or  super¬ 
sonic  speeds  (where  the  flow  is  highly  compressible). 
Nevertheless,  the  resemblance  of  the  data  to  the  La¬ 
placian  does  suggest  that  some  form  of  steady  jet  flow 
is  present.  Here,  jet  flow  is  meant  to  denote  the  fluid 
expansion,  as  opposed  to  the  collisionless  expansion 
of  material  from  a  circular  orifice. 

The  fluid  nature  of  the  flow  is  not  surprising,  since 
the  plasma  radius  r0  =  850  pm  is  much  greater  than  the 
ion-ion  collision  mean-free-path  (  =  1  pm  at  0.1  critical 
density  with  T  =  400  eV  as  determined  from  x-ray  spec¬ 
tra).  One  would  also  expect  steady  flow,  since  the  laser 
pulse-length  (3.  5  nsec  FWHM)  is  greater  than  the  sound 
transit  time  (~  2  nsec)  across  )„.  Experimental  evi¬ 
dence  that  the  flow  is  steady,  at  least  for  a  large  frac¬ 
tion  of  the  duration  of  x-ray  emission,  is  obtained  from 
streak  photography  of  images  of  3  2  w, emission  from  the 
it  ■' 4  surface.4  The  position  along  the  laser  axis  of  this 
surface  appears  stationary  for  the  duration  of  the  3  2U'0 
emission,  a  period  of  up  to  3  nsec  near  the  time  of 
peak  laser  intensity. 

The  fluid  flow  affects  the  far-field  particle  diagnostics 
such  as  ion  collectors  or  plasma  calorimeters;  we  be¬ 
lieve  that  these  effects  have  nut  been  previously  noted. 
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FIG.  1.  Schematic  illustration  of  experimental  apparatus. 


First,  the  angular  distribution  of  ablated  mass  observed 
by  these  diagnostics  may  be  solely  determined  by  the 
flow.  Due  to  the  shape  of  the  jet  near  the  target  sur¬ 
face,  a  fluid  velocity  direction  is  imparted  to  each 
fluid  element.  As  the  density  decreases  away  from  the 
target,  the  flow  eventually  becomes  collisionless,  and 
the  material  drifting  in  the  direction  set  by  the  jet  ex¬ 
pands  due  to  random  thermal  motion.  Since  directed 
ablation  velocities  are  generally  significantly  greater 
than  ion  thermal  velocities,  the  thermal  expansion  will 
not  greatly  increase  the  angular  range  of  material  flow 
set  by  the  jet.  Indeed,  charge  collector  and  plasma 
calorimeter  measurements  indicate  that  half  of  the 
ablated  mass  appears  in  a  cone  about  the  target  normal 
with  a  40°  half-angle5;  this  is  consistent  with  the  angu¬ 
lar  spread  observed  in  the  x-ray  images,  as  well  as 
with  that  given  by  the  simple  model  of  Fig.  2(b). 

Ion  detectors  are  furttier  affected  because  the  jet,  ex¬ 
cept  for  any  smearing  due  to  thermal  expansion  after 
the  flow  becomes  collisionless,  maps  regions  of  target 
surface  into  space.  Since  each  point  on  the  target  sur¬ 
face  feeds  mass  into  a  specific  solid  angle,  a  single  ion 
detector  far  from  the  target  samples  only  a  portion  of 
the  irradiated  region.  To  measure  the  parameters 
averaged  over  the  focal  sput,  then,  one  should  use  an 
array  of  detectors  at  a  range  of  angles.  The  error 
that  may  result  if  too  few  detectors  are  used  depends 
upon  the  uniformity  of  the  laser  illumination  at  the  tar¬ 
get.  the  ratio  of  the  directed  ablation  velocity  to  the 
plasma  thermal  velocity,  and  the  solid  angle  subtended 
by  the  detectors  used. 


(a)  (b)  (c) 


FIG.  2.  Comparison  of  observed  flow  and  simple  fluid  model, 
(a)  Three  x-ray  images  corresponding  to  pinhole  diameters  of 
13,  2®,  and  54  pm  .  with  laser  incident  from  right,  (b)  Pattern  of 
Laplacian  flow  from  circular  orifice.  Dotted  lines  are  con¬ 
stant  velocity-potential  surfaces  and  solid  lines  are  fluid 
streamlines,  (c)  Dependence  of  streamline  angle  8  MZ  -  r0/3  up¬ 
on  source  position  r  in  the  target  plane.  Solid  curve  is  for 
Laplacian  flow  pattern  fl  am  <  bl,  and  points  ( r )  with  error 
bars  are  data  from  x-ray  images  such  as  (a). 


effects  on  the  shape  of  this  flow.  Detailed  comparisons 
should  then  be  possible  with  the  results  of  multidimen¬ 
sional  hydrodynamics  calculations,  which  are  more  ap¬ 
propriate  to  the  problem  than  the  simple  model  pre¬ 
sented  here.  In  the  present  study,  streamlines  were 
identified  by  introducing  localized  sources  at  the  target 
surface;  it  should  be  noted  that  the  same  phenomenon 
may  be  observed  if  localized  sources  arise  due  to  unin¬ 
tentional  effects,  such  as  filamentation  of  the  incident 
laser  beam.6 
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In  conclusion,  a  new  diagnostic  technique  allows  the 
first  experimental  observation  of  the  streamlines  of 
material  flow  from  the  surface  of  a  laser-irradiated 
target.  Knowledge  of  the  flow  is  required  if  the  results 
of  planar  experiments  are  to  be  properly  scaled  to  the 
spherical  laser  fusion  problem.  Implications  of  the 
fluid  nature  of  the  flow  on  far-field  ion  diagnostics  in 
planar  experiments  are  noted,  as  is  a  resemblance  of 
the  observed  flow  pattern  to  a  very  simple  fluid  flow. 
Future  studies  will  include  parameteric  variation  of 
laser  intensity  and  spot  size,  to  assess  the  possible 
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Spot  spectroscopy:  Local  spectroscopic 
measurements  within  laser-produced  plasmas 

M.  J.  Herbst,  P.  G.  Burkhaiter,  J,  Grun,"  R.  R.  Whitlock,  and  M.  Fink61 

Sava l  Research  Laboratory.  Washington.  D.  C.  203  75 
(Received  8  March  1982,  accepted  for  publication  7  May  1982) 

Use  of  a  locally  embedded  tracer  in  laser-irradiated  solid  targets  yields  a  localized  source  of 
diagnostic  x-ray  line  radiation  in  the  blowoff  plasma.  This  technique  potentially  eliminates 
problems  of  chord  integration  over  regions  of  varying  density  and  temperature  in  an 
inhomogeneous  plasma,  and  reduces  complications  due  to  plasma  opacity  effects  in  the 
interpretation  of  spectra.  Spectra  obtained  in  an  experimental  test  of  this  new  technique  are  of  a 
quality  superior  to  those  obtained  from  standard  laser-produced  plasmas,  and  should  provide  the 
best  tests  to  date  of  spectroscopic  models  for  these  plasma  conditions. 

PACS  numbers:  52.50.Jm,  52.25. Ps 


INTRODUCTION 

Density  and  temperature  profiles  are  directly  related  to 
energy  absorption  and  transport  processes  in  plasmas 
produced  by  laser  irradiation  of  solid  targets.  Therefore, 
measurements  of  these  plasma  parameters  play  a  central 
role  in  experiments  addressing  the  physics  relevant  to 
laser-driven  nuclear  fusion.1  At  the  densities  and  tem¬ 
peratures  of  interest,  x-ray  spectroscopy  provides  one 
viable  method  for  performing  these  measurements.2-3 

There  are  two  major  drawbacks  to  the  use  of  standard 
spectroscopic  techniques  in  these  plasmas.  First,  the  mea¬ 
surements  are  integrated  along  a  diagnostic  line  of  sight 
through  an  inhomogeneous  plasma.  Second,  effects  of 
plasma  opacity  on  the  radiation  exiting  the  plasma  are 
non-negligible  at  the  higher  densities  of  interest.4-5  One 
could  use  straightforward  inversion  techniques  to  obtain 
spatially  resolved  cmissivities  from  the  chord-integrated 
data  if  the  plasmas  were  optically  thin  and  axisym- 
metric,3  but  the  additional  complication  of  spatially  vary¬ 
ing  opacity  makes  inversion  much  more  difficult.  More 
typically,  one  predicts  the  plasma  profile  using  a  hydro¬ 
dynamics  computer  code  and  compares  the  experimen¬ 
tally  observed  spectrum  with  those  calculated  for  chord 
integrations  across  the  inhomogeneous  plasma.6,7  One  is 
dependent  in  that  case  upon  the  correctness  of  the  hy¬ 
drodynamics  model  as  well  as  the  atomic  physics  model; 
one  would  obviously  prefer  to  make  direct  determinations 
of  density  and  temperature,  independently  of  the  hydro¬ 
dynamics. 

In  the  present  work,  a  new  spectroscopic  method  is 
described  which  circumvents  the  problem  of  plasma  in¬ 
homogeneity  and  reduces  the  effects  of  plasma  opacity. 
Spot  spectroscopy,  the  new  technique,  differs  from  stan¬ 
dard  spectroscopic  techniques  in  that  the  source  of  di¬ 
agnostic  radiation  is  locally  embedded  into  the  solid  tar¬ 
get  before  laser  irradiation.  As  recently  demonstrated,* 
the  tracer  material  which  is  ablated  by  the  laser  can  be 
confined  collisionallv  in  the  hydrodynamic  flow  from  the 


target;  this  yields  a  source  of  diagnostic  lines  which  is 
localized  within  the  blowoff  plasma.  In  Sec.  I,  we  explain 
the  target  fabrication  procedure  and  describe  a  method 
for  in  situ  measurement  of  the  spectrograph  magnifi¬ 
cation.  Data  obtained  with  the  new  technique  are  com¬ 
pared  with  standard  spectroscopic  data  in  Sec.  II,  which 
also  contains  a  discussion  of  caveats  and  remaining  ques¬ 
tions.  Finally,  we  assess  the  new  technique  in  Sec.  III. 

I.  EXPERIMENTAL  TECHNIQUES 
A.  Target  fabrication 

For  the  proof-of-principle  experiment,  we  choose  poly¬ 
styrene  targets  with  embedded  aluminum  tracers.  The 
tracer  is  so  chosen  because  collisional  radiative  equilib¬ 
rium  (CRE)  computer  models  suggest  the  utility  of  alu¬ 
minum  lines  for  measurement  of  temperatures  and  den¬ 
sities  of  interest.’  The  polystyrene  is  chosen  because  it 
is  a  common  and  easily  fabricated  target  material.  Also, 
its  lower  atomic  number  Z  increases  the  ratio  of  the 
aluminum  tracer  emission  to  the  surrounding  plasma 
emission  and  reduces  the  opacity  of  the  surrounding 
plasma  to  the  tracer  emission. 

A  two-step  process  is  used  to  fabricate  these  targets. 
The  first  step  is  the  evaporation  of  the  tracer  material 
onto  a  soaped-glass  substrate  using  the  standard  tech¬ 
nique  of  vacuum  evaporation  from  a  hot  filament  source. 
The  thin  film  of  TEEPOL  610  detergent10  on  the  glass 
prevents  the  aluminum  from  adhering  better  to  the  glass 
than  it  will  to  the  polystyrene.  To  achieve  local  deposition 
of  the  aluminum,  bimetal  masks  into  which  the  desired 
patterns  have  been  chemically  etched  are  magnetically 
clamped  in  intimate  contact  with  the  glass  substrate;  for 
the  present  experiment,  these  patterns  are  circular  spots 
of  various  diameters.  In  the  second  step  of  the  process, 
polystyrene  dissolved  in  butyl  acetate  is  cast  over  the 
aluminum-coated  glass  substrate.  When  dry,  the  com¬ 
posite  is  lifted  and  the  aluminum  is  embedded  in  the 
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polystyrene;  the  surface  of  the  tracer  is  flush  with  that 
of  the  plastic. 

B.  The  proof-of-principle  experiment 

Except  for  the  presence  of  the  tracer  in  the  target,  the 
experimental  arrangement  for  the  new  technique  is  much 
the  same  as  for  standard  x-ray  spectroscopy.  As  shown 
in  Fig.  1,  the  targets  are  irradiated  with  50-250J  of  Nd 
laser  radiation  (1.054  wavelength)  in  a  4-5-ns  pulse. 
We  choose  to  operate  at  two  focal  conditions  for  this 
experiment.  First,  at  the  best  focus  of  our  / /6  input  lens, 
we  use  focal  diameters  (90%  laser  energy  content)  near 
100  um,  and  spatially  averaged  irradiances  up  to  5 
X  10'4  W/cm2.  Because  larger  focal  diameters  are  re¬ 
quired  to  minimize  edge  effects  in  planar  experiments, 
targets  are  also  placed  in  the  quasinear  field  of  the  fo¬ 
cusing  lens;  intensities  between  1012  and  1013  W/cm2  are 
obtained  with  focal  diameters  of  800-1000 

The  aluminum  line  radiation  in  the  5-8  A  spectral 
region  is  detected  using  an  x-ray  crystal  spectrograph 
with  a  viewing  axis  parallel  to  the  target  surface.  Spectral 
dispersion  is  achieved  through  use  of  a  PET  crystal,  and 
spatial  resolution  in  the  direction  on  the  film  which  is 
orthogonal  to  the  spectral  dispersion  is  accomplished  with 
an  entrance  slit.  As  shown  in  Fig.  1,  the  slit  is  oriented 
to  yield  spatial  resolution  in  the  direction  along  the  laser 
axis;  resolution  in  the  other  two  spatial  dimensions  results 
from  the  knowledge  of  the  position  of  the  localized  tracer 
within  the  focal  spot.  The  slit  is  located  4-5  mm  from 
the  target,  which  represents  a  tradeoff  between  light  col¬ 
lection  efficiency  and  slit  survivability;  even  at  this  dis¬ 
tance,  the  25-um-thick  nickel  slit  substrate  is  occasionally 
destroyed  by  target  debris  during  a  laser  shot.  Another 
compromise  must  often  be  made  between  spatial  reso¬ 
lution  and  detectability  of  relevant  spectral  lines.  To  im¬ 
prove  spatial  resolution  along  the  laser  axis,  one  must 
reduce  the  slit  width,  sacrificing  spectral  intensity;  to 
recover  the  spectral  intensity,  one  may  choose  to  increase 
the  implanted  spot  diameter,  reducing  spatial  resolution 
in  the  other  two  dimensions.  For  the  purposes  of  this 
study,  slit  widths  between  4  and  30  and  tracer  spot 
diameters  between  40  and  250  pm  are  used.  The  spectra 
arc  recorded  on  Kodak  No-Screen  film,  for  which  a  film 
model"  based  upon  an  absolute  calibration12  is  available. 
A  beryllium  filter  of  nominally  15  (im  thickness  protects 
the  film  from  visible  exposure. 

C.  In  situ  measurements  of  spectrograph 
magnification 

In  principle,  one  may  use  the  measured  slit-to-filin  and 
slit-to-source  distances  to  calculate  the  magnification  of 
the  spectrograph.  However,  given  the  uncertainties  with 
which  the  exact  paths  of  diffracted  rays  of  various  wave¬ 
lengths  can  be  measured,  we  prefer  a  direct  in  situ  mea¬ 
surement  of  the  spectrograph  magnification. 

To  measure  the  magnification  in  silu.  we  devised  a 
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Flo.  1.  Experimental  configuration  for  spot  spectroscopy 

special  double  target.  An  aluminum  foil  is  placed  in  the 
usual  target  position,  mounted  to  a  glass  substrate  to 
prevent  target  bowing.  A  few  hundred  microns  in  front 
of  the  A1  foil,  toward  the  laser,  a  foil  of  copper  backed 
by  glass  is  mounted.  The  foils  are  aligned  so  that  part 
of  the  incident  laser  beam  strikes  the  edge  of  the  Cu  foil, 
producing  Cu  line  emission,  while  the  rest  misses  the  Cu 
and  propagates  the  additional  distance  to  the  A1  surface. 
The  spectrograph  image  then  contains  Cu  and  A!  lines, 
but  the  A1  lines  extend  beyond  the  end  of  the  Cu  lines 
by  a  distance  L  =  ML  corresponding  to  the  distance  L 
between  the  two  target  foils  multiplied  by  the  spectro¬ 
graph  magnification  M.  As  L  may  be  accurately  mea¬ 
sured  under  a  microscope,  M  is  experimentally  deter¬ 
mined  by  the  measurement  of  L'. 

One  potential  complication  with  this  procedure  arises 
unless  the  spectrograph  axis  is  closely  aligned  with  the 
target  surface;  if  it  is  not,  the  observed  separation  be¬ 
tween  the  Cu  and  A1  may  be  lessened  due  to  geometrical 
foreshortening.  To  assure  that  this  is  not  a  problem,  we 
add  a  flag  to  the  target,  a  2-mm-wide  and  125-um-thick 
glass  strip  on  the  side  between  the  irradiated  target  and 
the  spectrograph  slit.  This  flag  is  offset  toward  the  laser 
from  the  Cu  surface  and  is  mounted  with  its  broad  face 
parallel  to  the  Cu  and  A1  target  surfaces.  Given  its  lo¬ 
cation,  it  casts  a  shadow  on  the  spectrograph;  this  shadow 
is  observed  on  the  spectrum  as  a  gap  in  the  A1  and  Cu 
spectral  lines.  By  comparing  the  dimension  of  this  gap 
with  the  thickness  of  the  flag,  the  angle  at  which  the  flag 
and,  therefore,  the  targets  are  viewed  can  be  determined; 
any  angular  misalignment  causes  an  increase  in  the 
shadow  dimension  which  depends  on  the  known  flag 
width. 

II.  RESULTS  AND  DISCUSSION  • 

A.  Results  of  experiment  and  comparison  of 
techniques 

To  allow  a  comparison  to  be  made  between  spot  spec¬ 
troscopy  and  standard  spectroscopic  techniques,  we  use 
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(a)  (b) 


(c)  (d) 

Fig.  2.  Qualitative  comparison  between  data  obtained  with  an  Al  foil  target  and  da  .  acquired  using  a  CH  plastic  target  locally  embedded  with 
Al.  The  vertical  axis  for  (a)  through  (d)  corresponds  to  the  laser  axis,  with  the  laser  incident  from  above  and  the  target  surface  al  bottom  (a) 
is  obtained  with  Kodak  3490  film,  and  (b)  through  (d)  are  on  Kodak  No-Screen  Film.  (a)  x-ray  pinhole  camera  images,  obtained  with  an  array 
of  pinholes  (5-55  >im  diameter),  for  an  Al  foil  target.  Laser  intensity  /,  =  3  X  I0':  W/cm\  due  to  100  J  focused  to  1100-pm-diam  spot  (9l)'-c 
energy  content),  (b)  Spectrum  obtained  for  Al  foil  target  at  /0  a  8  X  I0IJ  W/cm*’.  witi:  1 10  J  focused  to  600-iim-diam  spot,  (c)  Pinhole  images 
as  in  (a),  but  for  a  CH  target  embedded  with  180-nm-diam  Al  spot.  225  J  of  laser  energy  focused  to  900-/im-diam  spot  yields  /0:lx  I0,; 
W/cm:.  (d)  Spectrum  from  same  shot  as  (c). 


both  simple  Al  foil  targets  and  CH  targets  with  embed¬ 
ded  Al  spots.  A  pinhole  camera  filtered  to  detect  photons 
with  energy  hv  Z  I  kcV  looks  at  the  blowoff  plasma  along 
the  same  axis  but  from  a  position  diametrically  opposite 
the  spectrograph. 

Typical  images  of  an  aluminum  foil  plasma,  obtained 
with  a  pentagonal  array  of  pinholes  with  diameters  be¬ 
tween  5  and  55  gm,  arc  shown  in  Fig.  2(a).  An  x-ray 
spectrum  obtained  with  an  aluminum  foil  target  can  be 
seen  in  Fig.  2(b),  although  for  somewhat  different  con¬ 
ditions  than  apply  for  Fig.  2(a).  The  horizontal  edge  at 
the  bottom  of  the  spectrum  is  the  target  surface,  from 
which  regions  of  strong  line  emission  extend  over  dis¬ 
tances  of  several  hundred  microns. 

The  x-ray  pinhole  images  and  spectrum  for  a  CH  tar¬ 
get  with  embedded  Al  can  be  seen  in  Figs.  2(c)  and  2(d), 
respectively;  this  data  is  analogous  to  that  shown  in  Figs. 
2(a)  and  2(b)  for  an  aluminum  foil  target.  Note  that 
Figs.  2(c)  and  2(d)  are  obtained  on  the  same  data  shot, 
but  that  the  conditions  of  this  shot  do  not  exactly  cor¬ 
respond  with  those  of  either  Fig.  2(a)  or  2(b).  The  com¬ 
parison  between  the  techniques,  therefore,  is  qualitative 
rather  than  quantitative. 

Inspection  of  the  pinhole  camera  images  reveals  that, 
indeed,  the  source  volume  for  the  Al  x-ray  emission  can 


be  significantly  reduced  through  the  use  of  a  CH  target 
with  an  embedded  Al  spot.  As  expected  on  the  basis  of 
Ref.  8,  a  collisionally  confined  channel  of  aluminum  flow 
from  the  embedded  tracer  can  be  identified  as  the  track 
of  stronger  x-ray  emissivity  in  the  images  of  Fig.  2(c). 
This  is  a  much  smaller  source  volume  than  the  aluminum 
foil  plasma  in  Fig.  2(a)  [the  20%  difference  in  laser  focal 
spot  diameter  between  Figs.  2(a)  and  2(c)  is  not  a  major 
factor  in  this  reduction].  Even  larger  reductions  in  source 
size  are  observed  when  smaller  embedded  spots  are  used. 
To  date,  we  have  embedded  spots  as  small  as  25  jim  in 
diameter,*  a  limit  imposed  by  the  availability  of  masks 
for  tracer  evaporation.  In  principle,  the  ultimate  limit  is 
set  by  the  ability  of  the  plasma  to  collisionally  confine 
the  tracer;  as  ion-ion  mean  free  paths  are  estimated  to 
be  less  than  1  pm  even  at  electron  densities  dow  n  to  10% 
of  critical  density  for  the  incident  laser,  there  is  consid¬ 
erable  room  for  further  source  reduction. 

A  dramatic  improvement  in  spectral  resolution  results 
from  the  reduction  in  source  volume.  With  an  Al  foil 
target,  as  in  Fig.  2(b),  the  spectral  widths  of  the  lines 
near  the  target  surface  arc  quite  broad;  the  widths  exceed 
those  expected  for  other  broadening  mechanisms,  and 
reflect  source  broadening  due  to  the  large  Al  plasma  ex¬ 
tent  transverse  to  the  laser  axis.  In  contrast,  line  widths 
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near  the  target  surface  are  much  narrower,  as  in  Fig. 
2(d),  when  a  CH  target  with  embedded  Al  is  used.  Sev¬ 
eral  sets  of  lines  which  appear  to  the  naked  eye  to  be 
totally  unresolved  in  Fig.  2(b)  are  clearly  isolated  in  Fig. 
2(d).  These  data  are  obtained  with  the  larger  focal  spot; 
a  similar  effect  is  observed  at  the  smaller  focal  spot. 

Two  other  important  results  of  the  reduced  source  vol¬ 
ume  are  that  the  effects  of  plasma  inhomogeneity  can  be 
virtually  eliminated  and  that  the  effect  of  plasma  opacity 
can  be  greatly  reduced.  Without  embedded  tracers,  one 
might  try  to  minimize  inhomogeneity  and  opacity  effects 
by  reducing  the  size  of  the  plasma  to  be  studied.  For 
laser-produced  plasmas,  this  would  be  accomplished  by 
minimizing  the  laser  focal  diameter.  One  could  not  ex¬ 
pect,  however,  to  approach  the  small  plasma  sizes  which 
seem  possible  with  embedded  tracers.  Even  if  one  could 
localize  the  laser  energy  to  micron  scales,  lateral  trans¬ 
port  of  energy  by  plasma  processes  would  tend  to  broaden 
the  spectroscopic  source,  increasing  the  effect  of  plasma 
opacity.  Additionally,  this  source  would  still  have  prob¬ 
lems  due  to  its  inhomogeneity;  one  would  still  be  sampling 
regions  of  widely  different  density  and/or  temperature. 
With  the  plasma  resulting  from  an  embedded  tracer,  one 
can  hope  to  obtain  not  only  very  small  source-plasma 
diameter,  minimizing  opacity  effects,  but  also  a  nearly 
homogeneous  source,  since  the  embedded  spot  diameter 
can  be  made  much  smaller  than  the  irradiated  focal  di¬ 
ameter.  The  reduction  in  opacity,  of  course,  is  realized 
only  if  the  tracer  is  embedded  in  a  target  material  with 
lower  opacity  than  the  tracer  at  the  wavelengths  of  in¬ 
terest. 

Given  that  the  source  volume  for  emission  is  much 
smaller  with  the  embedded  Al  tracer  than  for  an  Al  foil, 
one  might  have  anticipated  greater  problems  with  de¬ 
tectability  of  the  spectrum  than  were  actually  encoun¬ 
tered.  In  Fig.  2(d),  the  spectral  lines  are  detectable  sev¬ 
eral  hundred  microns  from  the  target  surface,  just  as  far 
as  they  are  in  the  spectrum  shown  for  the  pure  Al  target 
(though  under  somewhat  different  conditions)  in  Fig. 
2(b).  In  fact,  even  continuum  emission  is  detectable  near 
the  target  surface  in  Fig.  2(d);  this  continuum  is  not 
observable  when  plastic  targets  without  the  embedded 
spots  are  irradiated,  indicating  that  the  source  of  the 
continuum  is  the  Al  tracer  material.  For  the  stronger 
lines  in  the  spectrum,  such  as  the  heliumlike  and  hydro¬ 
genlike  resonance  lines  (la2  l5,0-la2/7l/>i  and  \s-2p,  re¬ 
spectively)  and  the  helium-like  intercombination  line  ( 1  s1 
'S0-\s2p'P,),  exposures  near  the  target  surface  tend  to 
saturate  the  film  except  on  shots  where  smaller  embedded 
spots  (40  *im  diameter)  and  narrower  slits  (4  wide) 
are  used.  Weaker  spectral  lines  such  as  the  higher  Ryd¬ 
berg  members  are  lost  in  this  limit.  Spot  spectroscopy 
does  appear  to  be  practicable,  then,  at  least  with  laser 
energies  comparable  to  ours. 

B.  Caveats  and  remaining  questions 

Spot  spectroscopy  has  advantages  over  standard  spec¬ 
troscopic  techniques  in  laser  produced  plasmas;  however, 
there  are  a  few  caveats  to  be  mentioned.  First,  source 


localization  requires  a  coll isional  plasma.  This  restriction 
may  prevent  use  of  the  technique  at  much  higher  laser 
intensities  or  at  much  longer  laser  wavelengths,  where 
higher  plasma  temperatures  and  the  generation  of  su- 
prathermal  ions  produce  longer  ion-ion  mean  free  paths, 
destroying  the  confinement  of  the  embedded  tracer.  Of 
course,  spot  spectroscopy  is  not  limited  to  laser-produced 
plasmas;  any  method  such  as  ion-beam  irradiation  or 
electrical  discharge  which  creates  collisional  plasma  from 
an  initially  solid  target  may  be  used  to  generate  the  spec¬ 
troscopic  source. 

A  second  caveat  concerns  the  ability  of  the  technique 
to  serve  as  a  plasma  diagnostic.  If  we  are  to  regard  the 
plasma  flowing  from  the  embedded  tracer  as  a  local  spec¬ 
troscopic  probe  for  density  and  temperature  measure¬ 
ments,  then  we  want  to  be  certain  that  conditions  in  the 
tracer  plasma  correspond  to  those  in  the  surrounding 
plasma.  Certainly,  the  visualizations  of  the  hydrody¬ 
namic  flow  pattern,  as  shown  in  Ref.  8,  strongly  suggest 
that  there  is  pressure  equilibrium  between  Al  tracers  and 
CH  plasmas.  In  the  future  we  will  determine  by  com¬ 
parison  with  other  diagnostics  when  there  is  also  density 
and  temperature  equilibration. 

Even  if  equilibration  is  shown  for  the  present  data, 
spectroscopic  models  show  that  each  tracer  element  may 
be  used  only  within  limited  ranges  of  density  and  tem¬ 
perature9;  to  diagnose  a  wide  range  of  parameters,  one 
must  use  a  variety  of  tracer  materials.19  The  required 
equilibration  must  be  demonstrated  for  each  of  these 
diagnostic  elements. 

Another  issue  is  the  effect  of  time  integration  upon  the 
data.  Much  as  with  chord  integration,  time  integration 
can  complicate  interpretation  of  the  spectra  if  plasma 
parameters  are  varying  during  the  x-ray  emission.  In  the 
present  experiments,  we  believe  that  the  laser  pulse  is 
varying  sufficiently  slowly  so  that  the  plasma  profiles  are 
nearly  in  a  steady  stale  during  the  x-ray  emission.  Ex¬ 
perimental  evidence'4  for  this  is  provided  by  streak  pho¬ 
tography  of  imaged  %  w0  emission  (where  u)0  is  the  in¬ 
cident  laser  frequency),  which  indicates  that  the  location 
of  the  nj 4  surface,  and,  therefore,  all  higher  density 
surfaces,  is  nearly  steady  for  a  few  nanoseconds  around 
the  peak  of  the  laser  pulse  (when  the  x-ray  emission  is 
observed  to  occur15).  Even  though  we  believe  time-inte¬ 
gration  effects  to  be  minimal  for  our  case,  one  can  hope 
to  obtain  time  resolution  with  spot  spectroscopy  in  one 
of  two  ways.  First,  one  may  replace  the  film  with  an  x- 
ray  streak  camera  to  obtain  a  time  history  on  a  single 
shot,  though  one  cannot  simultaneously  obtain  the  com¬ 
plete  spectrum  of  wavelengths  A  for  all  positions  r  relative 
to  the  target  surface.  An  alternative  method  is  to  deposit 
a  tracer  of  lesser  thickness  at  a  prescribed  depth  beneath 
the  surface  of  the  target  (this  involves  a  more  difficult, 
though  realizable,  target  fabrication  procedure).  By 
varying  the  depth  d  and  thickness  t  of  the  tracer,  one 
controls  the  time  interval  during  which  tracer  emission 
occurs.  On  a  single  shot,  then,  one  obtains  information 
for  all  A  and  but  at  only  a  single  time;  f  and  d  are 
varied  on  a  shot-to-shot  basis  to  generate  a  time  history. 
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III.  CONCLUSIONS 

The  advantages  of  the  new  technique  are  clear  for 
measurements  of  density  and  temperature  profiles  in 
plasmas  created  by  Nd-laser  interaction  with  solid  tar¬ 
gets.  Localization  of  the  source  of  line  radiation  allows 
spectroscopic  measurements  to  be  made  with  three  di¬ 
mensions  of  spatial  resolution.  Since  this  alleviates  the 
problem  of  chord  integration  over  regions  of  widely  vary¬ 
ing  density  and  temperature,  one  generally  will  not  need 
to  account  for  the  profile  along  which  one  is  integrating 
in  order  to  obtain  agreement  with  spectroscopic  models.5 
Therefore,  the  new  technique  also  provides  improved 
checks  of  spectroscopic  models.  As  added  bonuses:  (1) 
the  reduced  source  broadening  improves  spectral  reso¬ 
lution  in  the  data  obtained,  and  (2)  the  interpretational 
problems  raised  by  opacity  effects  in  the  plasma  are  much 
reduced,  especially  if  the  surrounding  target  material  has 
much  lower  Z  than  the  tracer. 

The  present  work  represents  a  significant  step  toward 
obtaining  improved  measurements  of  density  and  tem¬ 
perature  in  laser-produced  plasmas.  Detailed  analysis  of 
the  present  data,  including  extensive  comparisons  with 
spectroscopic  models,  is  the  subject  of  a  future  work. 
Further  experiments  are  also  required  to  fully  answer  all 
of  the  questions  raised  in  Sec.  I  IB. 
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EXPERIMENTAL  METHODS  FOR  STUDYING  THE  RAYLEICH-TAYLOR 
INSTABILITY  OF  ABLATIVELY  ACCELERATED  TARGETS* 
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S.P.  Obenschain,  B.H.  Ripin,  J.A.  Stamper  and  R.R.  Whitlock 

We  describe  new  diagnostic  methods  for  the  study  of  hydrodynamic 
Instabilities  in  ablatively  accelerated  targets.  These  methods  include  face- 
on  x-ray  backlighting  that  does  not  require  a  backlighting  laser  beam  (for 
growth  rate  measurement) ,  and  a  tracer  dot  technique  (for  tracking  ablation 
plasma  flow).  The  targets  in  our  experiments  are  periodically  perturbed  to 
provide  initial  conditions  for  the  growth  of  the  Rayleigh-Taylor  instability. 
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INTRODUCTION 


Designs  of  hollow,  directly  irradiated  inertial  fusion  pellets  are 

constrained  by  the  growth  rates  of  the  Rayleigh-Taylor  (R-T)  instability*  in 

the  ablation  region.  During  the  past  few  years  we  have  used  nonperturbed 

planar  targets  to  model  the  dynamics  of  such  fusion  pellets  early  in  the 
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implosion  phase,  before  spherical  convergence  effects  become  dominant.  We 

have  successfully  accelerated  these  targets  up  to  160  km/sec  while  keeping  the 

target  relatively  cool  (<  10  eV)  and  uniform  (15%).^’^  The  target  velocities 

were  consistent  with  ablation  pressures  and  ablation  rates  measured  on  the 
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thrust  (laser)  side  of  the  target.  ’  Moreover,  target  nonuniformities 
observed  in  these  experiments  were  shown  to  be  related  to  nonuniforaities  in 

<j  i 

the  irradiance  profile:  ’  We  have  not  as  yet  observed  any  nonuniformities  of 
accelerated  nonperturbed  targets  that  could  be  definitely  ascribed  to  a 
hydrodynamic  instability.^ 

To  study  the  development  of  hydrodynamic  instabilities  such  as  the  R-T 
during  the  ablation  phase  we  have  initiated  experiments  using  fcC**1eratc« , 
periodically-perturbed  targets.^  ***  Using  perturbed  targets** '**  for  R-T 
studies  has  these  advantages: 

-  By  seeding  the  instability  with  well-defined  initial  conditions 
experiment  and  theory  may  be  compared  more  directly  than  is  otherwise 
possible . 

-  The  known  periodicity  makes  it  easier  to  differentiate  between  R-T 
effects  which  result  from  the  initial  periodic  perturbation  and  other 
effects  which  are  not  related  to  the  perturbation. 
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-  By  varying  the  wavelength  and  amplitude  of  the  initial  perturbation  we 
control  the  linear  growth  rate  and  the  number  of  e-foldings  necessary 
to  reach  the  nonlinear  regime.  Consequently,  the  linear  and  nonlinear 
evolution  of  the  instability  is  easier  to  study. 

Our  targets  fall  into  two  classes  (Fig.  1),  targets  that  are  corrugated 
(wiggly)  with  a  constant  areal  mass  density  and  targets  that  have  a 
periodically  varying  mass  per  unit  area.  The  first  type  of  target  has  an 
initial  amplitude  perturbation  only.  The  second  type  of  target  develops  both 
an  amplitude  and  velocity  perturbation  early  in  the  laser  pulse  due  to  the 
differential  acceleration  across  the  target  surface.  Also,  in  the  second 
target  type  changes  in  the  areal  mass  density  across  the  target  surface  -  and 
hence  the  instability  growth  rate  -  may  be  more  easily  related  to  the  initial 
target  perturbations. 

To  get  measurable  growth  rates  we  accelerate  the  targets  to  10^  cm/sec 
with  a  5  nsec  FWUM,  1.06  pm  laser  pulse.  Under  these  conditions  about  23 
classical  e-foldings  (/  /kg  dt  )  are  expected  for  perturbation  periods  of 
10  pm:  Perturbations  with  periods  as  large  as  300  pm  are  expected  to  grow 
with  about  4  classical  e-foldings  (Fig.  1).  Thus,  even  if  actual  R-T  growth 
rates  are  reduced  from  /kg  by  a  factor  of  2  or  3  as  some  predict,*'*-1-’ 

measurable  growth  rates  can  still  be  achieved  in  our  parameter  regime. 

The  following  sections  describe  proof-of-principle  tests^’111  of  two 
experimental  methods  designed  to  study  the  linear  and  nonlinear  evolution  of 
the  R-T  instability.  The  first  method  is  a  face-on  x-ray  backlighting 

technique  with  which  we  hope  to  measure  areal  mass  nonuniformity  development 

and,  hence,  the  linear  growth  rates  of  the  instability.  The  unique  feature  of 
this  technique  is  that  it  does  not  require  a  separate  x-ray  source  or  a 
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separate  backlighting  laser  beam.  The  second  method  uses  x-ray  emission  from 

tracer  dots ^  to  map  out  the  ablation  plasma  flow  from  accelerated  perturbed 

targets.  This  method  nay  be  capable  of  detecting  the  vortex  shedding^  and 

Kelvin-Helnholtz  like  rollup  that  are  predicted  to  occur  on  the  ablation  side 

18 

of  unstable  targets. 

FACE-OX  X-RAY  BACKLIGHTING 

The  R-T  instability  in  a  periodically  perturbed  target  causes  mass  to 
flow  laterally  near  the  ablation  surface  from  the  thinner  to  the  thicker  parts 
of  the  target.  Consequently,  measurements  of  areal  mass  across  the  target 
surface  at  some  known  time  into  the  laser  pulse,  or  measurements  of  the  rate- 
of-change  of  areal  mass  density  across  the  target  surface,  can  be  related  to 
the  instability  growth  rates. 

We  deduce  areal  mass  variations  from  the  attenuation  of  an  x-ray  pulse 
that  has  passed  through  the  spike  and  bubble  structure  in  a  target  predicted 
to  be  R-T  unstable.  These  x  rays  originate  from  a  thin  layer  of  magnesium 
that  has  been  burled  within  a  periodically  perturbed  target  made  of  carbon. 
In  these  initial  experiments  the  depth  at  which  magnesium  is  buried  is  chosen 
so  that  the  5  nsec  FWHM  laser  pulse  burns  through  to  the  magnesium  at  2.5  nsec 
past  the  peak  laser  intensity.  The  sequence  of  events  is  as  shown  in  Fig. 
2:  First,  the  outer  carbon  layer  Is  heated  and  ablates,  accelerating  the 

target  but  emitting  relatively  few  x  rays  above  1  KeV.  Once  the  outer  carbon 
layer  is  depicted  (and  hopefully  the  R-T  is  well  developed)  the  thin  magnesium 
layer  rapidly  heats  up  emitting  the  short  burst  of  >  1  KeV  x  rays  that  are 
used  to  radiograph  the  target. 

Figure  3  shows  the  experimental  setup.  The  C/Mg/C  perturbed  target  is 
viewed  by  three  pinhole  cameras  and  a  PIN-diode  all  filtered  with  1/2  mil 


U 


beryllium*  (An  x-ray  streak  camera  may  also  be  used  in  the  future.)  One  of 
the  two  pinhole  cameras  on  the  laser  side  of  the  target  provides  a  two- 
dimensional  image  of  the  x-ray  flash  source;  the  second,  with  a  carbon  step 
filter  in  front,  provides  a  calibration  of  the  film  density  as  a  function  of 
the  carbon  areal  mass  that  the  x  rays  pass  through.  The  third  camera,  at  the 
target  rear,  records  the  radiograph  of  the  target.  The  diode  is  used  to 
measure  the  duration  of  the  Mg  x-ray  pulse  and  its  timing. 

This  face-on  x-ray  backlighting  method  has  some  very  nice  features:  No 
separate  backlighting  laser  beam  is  needed;  alignment  of  diagnostics  is  easier 
than  in  a  comparable  experiment  with  a  separate  backlighting  source;  the 
diagnostics  are  simple,  and  the  system  is  self-calibrating  on  each  shot*  It 
is  also  possible  to  account  for  x-ray  source  nonuniformities  during  data 
reduction  by  comparing  the  x-ray  source  photograph  with  the  radiograph.  On 
the  other  hand,  the  target  Is  somewhat  complicated  and  there  may  be  some  space 
versus  time  distortion  due  to  possible  non-simultaneity  of  the  Mg  emission 
across  the  target. 

Our  first  experiments  tested  whether  the  magnesium  pulse  was  sufficiently 
brighter  than  the  less  intense  but  longer  duration  carbon  x  rays,  whether  the 
x-ray  film  density  is  a  sufficiently  sensitive  function  of  the  carbon  target 
areal  mass  density,  and  whether  the  magnesium  x-ray  pulse  is  sufficiently 
short  to  make  a  meaningful  growth  rate  measurement  with  a  pinhole  camera.  To 
check  the  ratio  of  magnesium  to  carbon  x-ray  emission  we  shot  a  perturbed 
target  in  whLch  the  buried  magnesium  layer  and  the  overlaying  carbon  had  the 
shape  of  a  disk  of  about  the  same  diameter  as  the  laser  focal  spot.  The  laser 
focal  spot  was  offset  with  respect  to  this  disk  so  that  a  part  of  the  beam  was 
incident  on  pure  carbon,  exciting  carbon  x  rays  only,  and  a  part  was  incident 
on  the  C/Mg/C  layer.  The  laser  irradiance  and  thicknesses  of  C/Mg  layers  were 
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chosen  so  that  the  magnesium  would  light  up  at  about  2  nsec  past  the  peak  of 
the  pulse  and  deplete  1/2  nsec  later.  Figure  4  describes  the  experimental 
parameters  and  shows  photographs  of  the  x-ray  source  as  well  as  the  resulting 
radiograph.  The  radiograph  shows  a  clear  picture  of  the  perturbed  target. 
The  exposure  levels  caused  by  x  rays  from  the  layered  C/Mg/C  region  were  about 
6  times  greater  than  exposures  from  the  pure  C  region  for  a  pinhole  diameter 
of  6-8  um  and  a  pinhole  camera  magnification  of  2.5.  The  signal/noise  ratio 
in  this  experiment  is  therefore  adequate. 

It  is  interesting  to  note  that  the  radiograph  in  Fig.  4  shows  some  large- 
scale  nonuniforraities  that  are  not  correlated  with  the  initial  target 
perturbation.  Some  of  these  nonuniforraities  appear  in  the  x-ray  source  photo 
as  well  and  are  probably  related  to  imperfections  in  the  carbon/magnesiura 
coatings.  These  types  of  nonuniformities  can  be  accounted  for  in  data 
reduction;  also,  improvements  in  the  fabrication  of  the  layered  targets  can  be 
made.  Other  small-scale  nonunif ormities  appear  in  the  radiograph  but  not  in 
the  x-ray  source.  These  nay  be  caused  by  the  instability  mechanism  itself  or, 

perhaps,  by  variations  in  target  thickness  caused  by  a  nonunif orm  laser 

beam.  In  either  case,  two-dimensional  pinhole  pictures  will  be  necessary  to 
monitor  and  understand  the  phenomenon. 

The  functional  relationship  between  the  film  density  and  the  areal 
density  of  carbon  that  the  x  rays  passed  through  depends  on  the  location  of 
the  film's  H-D  curve  -  or  equivalently  on  the  source  intensity,  spectrum, 

duration,  the  amount  of  intervening  carbon,  pinhole  diameter,  camera 

magnification  and  so  on.  Because  so  many  variables  affect  the  result,  in  the 
future  we  intend  to  take  in  situ  calibration  measurements  on  each  shot  through 
an  array  of  pinholes. 

For  this  experiment  calibration  was  done  on  separate  shots.  An  example 
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of  a  film  density  versus  carbon  areal  density  curve  for  an  8  pm  diameter 
pinhole,  camera  magnification  of  2  and  other  conditions  similar  to  Fig.  4  is 
shown  in  Fig.  5.  Under  these  conditions  the  useful  range  of  film  density 
corresponded  to  areal  mass  density  less  than  1.5  mg/cm^.  The  ratio  of  C/Mg 
layer  to  C  x-ray  induced  density  (signal/noise)  was  high.  A  larger,  13  pm 
diameter,  pinhole  showed  reduced  signal/noise  but  carbon  thicknesses  up  to  2 
mg/cra  were  easily  distinguishable. 

We  measured  the  duration  of  the  magnesium  x-ray  pulse  using  a  filtered 

MRD-500  photodiode,  which  has  a  1/2  nsec  rise  time.  However,  since  this 

measurement  is  spatially  integrating,  nonuniformities  in  the  laser  pulse  which 
ablate  the  target  at  different  rates  as  well  as  imperfections  in  the 

carbon/magnesium  layers  broaden  the  diode  pulse;  the  result  is  an  upper  limit 
on  the  x-ray  pulse  duration  for  a  given  location  on  the  target.  This  upper 
limit  was  measured  to  be  2  nsec  FWHM.  In  future  experiments  we  will  deposit 
the  magnesium  layer  locally  within  the  focal  spot,  improve  the  quality  of  our 
targets,  and  improve  the  temporal  resolution  so  as  to  measure  the  actual  x-ray 
pulse  duration. 

What  if  the  magnesium  x-ray  pulse  is  not  significantly  shorter  than  2 
nsec  -  can  meaningful  measurements  still  be  made?  The  answer  is  yes  - 

although  the  flexibility  of  this  backlighting  method  would  be  reduced.  For 
example,  during  the  latter  part  of  the  laser  pulse  the  number  of  classical 
integrated  R-T  growth  periods  (  J  Vkg  dt}  changes  by  only  20%  since  the 
laser  power,  and  therefore,  target  acceleration  are  falling  rapidly  (Fig. 
2).  A  measurement  of  such  a  slowly  varying  quantity  is  possible  with  a  2  nsec 
x-ray  pulse.  But  a  2  nsec  FWHM  x-ray  pulse  is  too  long  to  allow  multiple 
flash  experiments  where  laterally  offset  magnesium  layers  are  buried  at 
varying  depths  in  the  target  to  temporally  and  spatially  resolve  the  evolution 
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of  areal  mass  perturbations  with  pinhole  camera  diagnostics.  An  x-ray  streak 
camera  would  have  to  be  used  under  such  circumstances. 


TRACER  METHOD  FOR  FLOW  VISUALIZATION 

The  ablation  plasma  flow  from  a  R-T  stable  and  R-T  unstable  target  is 

predicted  to  be  very  different.  If  a  R-T  stable  target  is  irradiated  by  a 

relatively  uniform  laser  beam  the  ablation  plasma  flows  in  a  fan  pattern  which 

resembles  that  of  an  incompressible  fluid  squirting  from  a  circular 

nozzle.1®  xn  a  linearly  unstable  R-T  target,  however,  non-colinear  pressure 

and  density  profiles  in  the  ablation  plasma  are  thought  to  generate  vorticity 

near  the  ablation  surface  causing  the  ablating  plasma  to  be  shed  from  the 

'  '  17 

target  in  a  circulating  pattern  resembling  that  of  a  Karman  vortex  street. 

In  nonlinearly  unstable  R-T  targets  the  vorticity  is  thought  to  be  advected  up 

the  sides  of  R-T  spikes  producing  a  Kelvin-Helmholtz-like  rollup  of  the  spike 

tips.  *  The  degree  of  rollup  may  depend  on  the  de'ails  of  initial 

perturbation  (square  vs  sinusoidal),1^  and  the  wavelength  of  the 

1 8 

perturbation.  These  phenomena  may  be  responsible  for  the  reduction  of  the 
linear  R-T  growth  rates  below  classical  (/kg)  and  for  the  eventual  saturation 
of  the  instability.1^’1® 

Figure  6  shows  the  nonlinear  Kelvin-Helmholtz-like  rollup  of  the  spike 

tips  of  a  strongly  perturbed  R-T  unstable  target  as  predicted  by  the  NRL  FAST 
1 8 

2D  simulation.  The  corresponding  flow  pattern  is  shown  in  Fig.  7. 

We  are  attempting  to  apply  tracer  methods  used  successfully  in  other 
experiments ^ to  detect  the  above  phenomena  in  periodically  perturbed 
targets  predicted  to  be  R-T  unstable.  The  experimental  arrangement  is  shown 
in  Fig.  8.  A  row  of  aluminum  or  silicon  tracer  dots  is  embedded  into  or 
deposited  on  the  surface  of  a  periodically-perturbed  polystyrene  (CH)  or 
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carbon  target.  (For  the  experiments  below  the  tracer  material  was  aluminum 
deposited  on  a  carbon  surface.)  When  the  target  is  irradiated  the  tracers  and 
the  target  surface  heat-up  and  ablate  away.  Since  the  tracers  are  much 
stronger  emitters  of  x  rays  in  the  spectral  band  of  the  pinhole  camera  (>^  1 
keV)  than  either  carbon  or  polystyrene,  the  location  of  tracer  ions  in  the 
ablation  plasma  flow  is  identified  through  their  emission.  We  assume  that  the 
tracers  do  not  significantly  perturb  the  dynamics  of  the  ablation  plasma 
flow.  This  assumption  is  reasonable  since  the  speed  of  Al,  Si,  and  CH  or  C 
target  ions,  measured  with  time-of-f light  ion  collectors,  is  the  same.  Also 
for  a  uniform  electron  density  in  the  ablation  plasma  the  fully-ionized  ion 
mass  density  for  Al  and  CH  (C)  targets  differ  by  only  12%  (4%):  Si  and  C 
fully-ionized  ion  mass  densities  differ  only  by  0.2%.  Embedded  Si  tracers 
should  not  perturb  the  target  dynamics  significantly  either  since  Si  and  C 
densities  are  about  equal  (2.32  vs  2.25  gm/cm  ). 

Examples  of  the  targets  we  use  are  shown  in  Fig.  8.  The  target  types 
are:  wiggly  CH  targets  with  Al  dot  tracers,  areal  mass  perturbed  carbon 
targets  with  Al  dot  tracers,  and  areal  mass  perturbed  carbon  targets  with 
rectangular  aluminum  tracers  that  lie  only  on  the  thicker  parts  of  the  target 
and  have,  therefore,  the  same  periodicity  as  the  target  perturbations. 

For  the  experiments  described  below  the  targets  were  viewed  edge-on  with 
a  pinhole  camera.  Since  this  gives  no  time  resolution  the  expected  circular 
plasma  motion  will  show  up  as  a  smearing  of  the  tracer  image.  In  the  future 
we  Intend  to  also  use  an  x-ray  streak  camera  to  get  the  desired  temporal 
resolution  and  spatially  two-dimensional  information. 

Our  experiments  show  that  the  plasma  flow  from  perturbed  and  nonpcrturbed 
targets  can  indeed  be  very  different.  In  the  nonperturbed  case  each  tracer 
dot  naps  into  a  distinct  flux  tube  as  the  ions  flow  away  from  the  target  (Fig. 
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8,9).  In  contrast,  flux  tubes  from  perturbed  targets  merge  together 

Indicating  that  lateral  plasma  flow  occurs.  This  is  seen  most  clearly  in  Fig. 
9  which  resolves  the  flow  from  a  200  pm  period,  2:1  areal  mass  perturbed 
target  using  six  tracer  dots  per  period.  A  similar  target  with  a  50  pm  period 
(Fig.  10)  shows  such  a  merging  of  flux  tubes  also  but  with  poorer  resolution 
(2  tracer  dots/period).  The  lateral  plasma  flow  deduced  from  these  figures 
resembles  that  expected  from  R-T  unstable  targets,  but  more  mundane  mechanisms 
associated  with  dynamic  breakup  of  these  heavily  perturbed  targets  cannot  be 
excluded.  Future  experiments  with  improved  time  resolution  and  more  extensive 
parameter  variation  will  differentiate  between  the  possible  flow  mechanisms. 

It  is  interesting  to  note  that  the  x-ray  emission  from  the  thicker 
portions  of  the  heavily  perturbed  target  in  the  edge-on  photos  of  Fig.  9  and 
Fig.  10  is  greater  than  the  emission  from  the  thinner  portions.  The  same 
phenomenon  is  seen  in  Fig.  11  which  shows  a  frontal  view  of  the  emission  from 
a  heavily  perturbed  carbon  target.  Since  the  height  of  the  initial 
perturbation  step  in  all  these  targets  is  much  smaller  than  the  critical  to 
ablation  distance,  these  brightness  variations  appear  to  be  due  to  phenomena 
occurring  during  the  target  acceleration.  For  example,  similar  phenomena 
caused  by  a  R-T  spike  shielding  a  R-T  bubble  from  laser  radiation  have  been 
observed  in  numerical  simulations.^  Effects  not  related  to  the  R-T  may  also 
explain  our  observations:  For  instance,  if  the  thinner  part  of  the  target 
accelerates  ahead  of  the  thicker  part,  the  critical  to  ablation  distance  may 
be  smaller  near  the  thicker  part  than  near  the  thinner  part.  If  this  were  so, 
laser  energy  would  transport  preferentially  to  the  thicker  part  causing  it  to 
emit  more  brightly.  In  any  case  care  must  be  taken  experimentally  that  the 
target  self  emission,  which  can  acquire  the  periodicity  of  the  imposed 
perturbations,  does  not  distort  areal  mass  measurements.  Cases  in  which  the 


10 


i 


ratio  of  backlighter  to  target  x-ray  emission  is  not  very  large  or  the 
contrast  of  the  backlighter  x  rays  as  they  pass  through  the  thicker  and 
thinner  parts  of  the  target  is  not  very  large  require  special  care. 

Very  interesting  and  unexpected  plasma  flow  patterns  were  seen  in  a  few 
cases.  The  most  interesting  q f  these  is  shown  in  Fig.  12.  Here  the  target 
was  an  areal  mass  perturbed  type  with  tracers  that  were  in  phase  with  the 
150  pm  perturbation  period.  The  plasma  flow  from  most  of  this  target  exhibits 
the  lateral  expansion  that  has  been  observed  in  Fig.  9  and  Fig.  10.  But  the 
flow  from  one  location  behaves  very  differently:  The  plasma  near  the  target 
surface  in  this  location  is  compressed  and  is  dimmer  than  plasma  from  other 
comparable  regions  of  the  target.  Further  from  the  target  surface  a  dramatic 
decompression  occurs  and  coincidently  the  plasma  emission  increases.  Still 
further  from  the  target  surface  the  plasma  becomes  compressed  and  dimmer  but 
this  time  the  compression  ratio  is  much  smaller.  (The  second  compression 
cycle  is  seen  on  the  negative  though  it  may  not  be  clearly  visible  on  the 
figure.)  This  phenomena  is  all  the  more  amazing  since  it  is  seen  in  a  time- 
integrated  pinhole  picture.  The  flow  pattern,  therefore,  resembles  a 
standi ng-wa ve -like  phe nome no n . 


CONCLUSION 

We  are  developing  methods  to  measure  the  growth  rate  and  ablation 
characteristics  of  hydrodynamically  unstable  targets.  In  our  experiments 
perturbations  which  may  seed  the  instability  are  purposely  introduced  in  the 
targets.  The  first  experiments  using  a  face-on  x-ray  backlighting  method  and 
a  tracer-dot  method  have  been  described.  Both  methods  show  pronise. 
Parametric  studies  will  be  attempted  in  the  near  future. 
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FIGURE  CAPTIONS 


Fig.  I 


Fig.  2 


Fig.  3 


Fig.  A 


Hydrodynamic  instabilities  are  studied  using  accelerated, 

periodically  perturbed  targets.  Amplitude  modulated  targets  with 

constant  areal  mass  and  areal  mass  modulated  targets  are 

illustrated.  The  number  of  classical  Rayleigh-Tay lor  e-foidings 

during  the  laser  pulse  is  shown  as  a  function  of  final  target  speed 

to  illustrate  where  measurable  growth  rates  may  be  achieved  in  the 

experiment.  X  is  the  wavelength  of  the  Rayleigh-Tay lor  mode. 

Schematic  of  the  temporal  behavior  of  the  laser  pulse,  magnesium  x- 
t _ 

ray  flash,  and  /  /kg  dt  integrated  over  the  laser  pulse.  The 

—oo 

exact  duration  of  the  x-ray  flash  is  uncertain  as  indicated  by  the 
dashed  lines.  If  the  x-ray  flash  is  bright  and  of  short  duration 
the  instability  growth  may  be  temporally  resolved  using  pinhole 
cameras. 

Setup  of  the  face-on  x-ray  backlighting  experiment.  The  function 
of  each  diagnostic  is  described  in  the  text.  The  arrow  indicates 
the  incident  laser  direction. 

Radiograph  of  a  perturbed  target  and  an  image  of  the  x-ray 

source.  On  this  shot  the  buried  magnesium  layer,  which  has  the 

shape  of  a  disk,  is  offset  with  respect  to  the  focal  spot.  The 

oval  brighter  images  are  caused  by  the  x  rays  emitted  by  the 

magnesium  layer.  The  dimmer  part  of  the  image  is  caused  by  the 
emission  of  carbon  x  rays  only.  The  parameters  of  this  shot  are: 
target  perturbation  period  =  50  pm  with  peak  (valley)  thickness  = 
8.7  pm  (5. A  pm),  upper  carbon  over  layer  thickness  =  1  pm,  and  I'.g- 
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Fig.  5 

Fig.  6 

Fig.  7 

Fig.  8 

Fig.  9 


1*1  O 

layer  thickness  =  0.1  pm;  average  irradiance  =  1  x  10  J  W/cm  . 
Optical  density  of  x-ray  film  as  a  function  of  the  areal  density  of 
a  cold  carbon  filter  for  an  8  pm  pinhole  and  a  camera  magnification 
of  2.  The  irradiance  and  C/Mg  layer  thicknesses  were  similar  to 
those  of  Fig.  4.  The  Mg/C  x-ray  source  curve  and  the  C  x-ray 

source  curve  were  obtained  on  separate  shots. 

Density  contours  of  a  nonlinearly  unstable  Rayleigh-Taylor  target 
predicted  by  the  NRL  Fast  2D  code.  The  target  parameters  and  tine 
of  the  "snapshot"  are  shown  in  the  figure. 

Predicted  ablation  plasma  flow  patterns  for  the  same  conditions  as 
in  in  Fig.  6. 

Experimental  arrangement  and  targets  used  to  visualize  the  ablation 
plasma  flow  from  periodically  perturbed  targets.  Upper  left:  the 
experimental  arrangement.  Lower  left:  fan-like  plasma  flow  from 
an  accelerated  planar  target.  Upper  right:  corrugated  plastic 
target  with  aluminum  tracer  dots.  Lower  right:  carbon  targets 

with  perturbed  areal  mass  and  tracer  dots.  Target  dimensions  are 
exaggarated  for  clarity  in  these  sketches. 

Plasma  flow  from  a  flat  (left)  and  a  perturbed  target  (right).  The 
perturbation  period  is  200  urn  with  an  areal  mass  ratio  (thick/thin) 
of  2.  There  are  approximately  6  tracer  dots  per  perturbation 

period.  Unlike  the  flat  target  case  where  flux  tubes  of  individual 

tracer  dots  are  visible,  the  individual  flux  tubes  in  the  perturbed 
case  are  not  distinguishable.  The  perturbed  target  photos  were 
taken  on  the  same  shot  through  three  different  diameter  pinholes. 
The  smallest  pinhole  photo  is  on  the  bottom;  the  largest  pinhole 
photo  on  the  top. 


Fig.  10  Plasma  flow  from  a  target  perturbed  with  a  50  pm  period.  The 
photographs  of  the  perturbed  target  were  taken  on  the  same  shot 
through  three  different  diameter  pinholes.  The  smallest  pinhole 
photo  is  on  the  bottom;  the  largest  pinhole  photo  on  the  top. 
There  are  about  two  tracer  dots  per  perturbation  period. 

Fig.  11  Pinhole  photograph  of  >_  1  KeV  x-ray  emission  taken  from  the  laser 
side  of  a  heavily  perturbed  carbon  target. 

Fig.  12  Unusual  plasma  flow  from  a  perturbed  target  diagnosed  with  aluminum 
tracers  matched  to  the  perturbation  period.  Note  the  compression 
and  expansion  of  ablation  plasma  on  the  right  laser  side  of  the 
photograph. 
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MEASUREMENTS  OF  LASER  COUPLING  AND 
PLASMA  PROFILES  IN  LONGER-SCALELENGTH  PLASMAS 

I.  INTRODUCTION 

To  maximize  applicability  of  current  laser-plasma  coupling  experiments  to  high-gain  laser  fusion, 
these  experiments  should  include  three  very  important  characteristics.  First,  laser-target  coupling  stu¬ 
dies  should  be  performed  in  plasmas  with  longer  density,  temperature,  and  velocity  scalelengths  to 
minimize  the  extrapolation  which  must  be  made  from  present  experiments  to  the  multimillimeter- 
scalelength  plasmas  anticipated  for  high-gain  pellets.  Second,  measurements  should  be  made  over  a 
wide  range  of  scalelengths  to  provide  a  number  of  benchmarks  for  the  plasma  and  hydrodynamics  codes 
which  are  used  to  make  the  extrapolation.  This  is  important  because  the  absorption  physics  in  long- 
scalelength  plasmas  may  be  qualitatively  (as  well  as  quantitatively)  different  from  that  in  short- 
scalelength  plasmas;  the  various  processes  which  determine  absorption  efficiency  depend  differently  on 
plasma  scalelengths,  and  the  processes  which  determine  the  absorption  efficiency  in  short-scalelength 
plasmas  might  not  be  the  same  processes  which  are  most  important  in  long-scalelength  plasmas. 
Finally,  it  is  important  that  a  complete  set  of  plasma  scalelength  measurements  be  made  in  these  exper¬ 
iments.  If  one  measures  absorption  or  scattered  light  and  does  not  specify  the  scalelengths.  free  param¬ 
eters  are  available  which  may  be  used  in  a  hydro-code  calculation  to  fit  the  measurements;  for  example, 
the  flux  limiter  may  be  varied  (which  alters  the  plasma  profiles)  to  obtain  agreement  with  the  experi¬ 
ment.  In  this  case,  the  physics  assumptions  in  the  calculation  are  not  being  properly  tested,  and  no 
confidence  is  developed  in  the  ability  of  the  hydro-code  to  predict  laser-target  coupling  in  reactor-sized 
plasmas. 

Standard  measurements  of  laser  coupling  with  solid  targets  exhibit  weaknesses  on  all  three  points. 
As  illustrated  in  Fig.  1,  these  generally  involve  focusing  the  laser  to  a  small  spot  in  order  to  obtain  the 
higher  intensities  of  interest  for  high-gain  scenarios.  Planar  geometry,  rather  than  spherical,  is  com¬ 
monly  used  because  somewhat  longer  plasma  scalelengths  can  be  produced  with  fixed  laser  power  and 
intensity.  Even  so,  scalelengths  produced  by  a  small  focal  spot  are  still  much  shorter  than  reactor 
scalelengths  because  of  the  divergent  How  of  plasma  from  the  finite  planar  spot.  Moreover,  variations 
of  laser  intensity  and  plasma  scalelength  cannot  be  made  independently.  If  the  laser  power  is  varied  to 
change  the  irradiance  at  constant  spot  size,  the  plasma  scalelength  varies  also;  if  laser  spot  size  is  varied 
to  control  plasma  scalelength,  laser  irradiance  and  focal  quality  change  as  well.  Measurements  of  den¬ 
sity  and  temperature  profiles  have  not  been  made  routinely  in  these  experiments,  and  techniques  for 
velocity  profile  measurement  have  not  been  developed. 

Significant  physics  results  have  been  obtained  in  previous  experiments  designed  to  address  laser 
coupling  in  longer  scalelength  plasmas;  these  experiments  have  also  had  limitations.  Using  underdense 
gas  targets  to  eliminate  complicating  elfects  due  to  the  presence  of  a  critical  surface,  investigators  have 
made  positive  identifications  of  stimulated  Brillouin  backscatter1 4  and  sidcscatter'  and  of  stimulated 
Raman  scattering. ”  They  have  also  observed  the  saturation  of  Brillouin  backscatter ' v  and  sidcscatter.' 
Measurements'1  and  calculations11  suggest  that  coupling  to  these  plasmas  is  dependent  upon  the 
dvnamics  of  ionization  and  healing  waves,  as  well  as  upon  the  hydrodynamic  expansion  lateral  to  the 
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laser  axis;  thus,  the  dynamics  of  the  gas  target  plasmas  are  as  complicated  as  the  dynamics  of  solid  tar¬ 
get  plasmas  but  are  less  similar  to  the  dynamics  of  reactor  pellet  plasmas.  Furthermore,  the 
scalelengihs  in  these  underdense  gas  targets  are  often  so  long  that  they  can  be  treated  by  the  theory  of 
parametric  instabilities  in  homogeneous  plasmas;  therefore,  one  does  not  check  the  inhomogeneous 
theory  of  relevance  to  reactor  pellets.  More  recent  experiments  with  inhomogeneous,  overdense  gas 
targets12'14  overcome  this  latter  problem.  However,  the  important  efTect  of  the  flow  velocity  gradient 
along  the  laser  axis  is  still  not  present. 

At  NRL,  longer-scalelength  solid-target  experiments  in  progress  are  designed  to  include  effects  of 
flow  velocity  gradients  while  overcoming  the  abovementioned  deficiencies  of  standard  solid-target 
experiments.  To  produce  longer  scalelength  plasmas,  two  separate  beams  from  our  Nd  laser  (X  =  1.054 
/am)  may  be  focused  in  different  planes,  as  shown  in  Fig.  2.  One  beam,  operating  with  up  to  250  J  in  a 
3-5  ns  pulse,  is  brought  to  a  large  diameter  (=1  mm)  spot  at  the  target  to  produce  plasmas  with  density 
scalelengihs  L„  =  /i/V n  of  250-600  /am  at  rt,/ 10.  The  second  beam  is  tightly  focused  inside  the  first  to 
give  the  high  intensities  of  interest  for  laser-plasma  coupling  experiments.  By  reducing  the  pulse  dura¬ 
tion  of  the  second  beam  to  tl  <  0.5  ns  (with  laser  energies  as  high  as  25  J),  we  can  minimize  the  per¬ 
turbation  to  the  long-scalelength  hydrodynamic  flow  of  the  first  beam  that  is  caused  by  this  higher 
intensity  region.  Both  high  intensity  and  long  scalelength  can  be  achieved  in  these  experiments,  and 
the  two  parameters  may  be  varied  independently  since  the  scalelength  is  determined  by  the  larger  spot 
size  and  the  intensity  by  the  smaller  spot  size.  This  new  method  for  performing  longer-scalelength  cou¬ 
pling  experiments  has  evolved  from  prepuise  techniques  pioneered  by  NRL  in  the  mid-1970's.15  As  will 
be  discussed  later  in  this  paper,  we  are  developing  new  plasma  diagnostic  techniques  to  make  improved 
measurements  of  density  and  temperature  profiles,  as  well  as  the  first  measurements  of  velocity 
profiles. 

II.  PRELIMINARY  LONGER-SCALELKNC  fH  EXPERIMENTS 

Preliminary  observations  of  laser-plasma  coupling  have  been  made  using  this  two-beam 
configuration  to  produce  longer  scalelengihs.  As  indicated  in  Fig.  3,  a  significant  increase  in  backscatter 
through  the  lens  is  observed  when  the  large-focus  beam  is  combined  with  the  tightly  focused  beam;  this 
occurs  both  for  long  pulselcngth  and  short  pulselength  operation  of  the  higher  intensity  beam,  although 
the  maximum  intensity  achievable  in  the  former  case  is  iewer  due  to  laser  protection  considerations. 
Despite  the  large  increase  observed  in  backscatter  from  the  carbon  targets  used  in  these  experiments, 
the  maximum  backscatter  is  10  percent;  thus,  the  backscatter  accounts  for  only  a  small  fraction  of  the 
incident  energy.  Measurements  of  scattered  light  at  other  angles,  which  have  not  yet  been  reduced, 
may  alter  this  picture. 


While  the  increase  in  backscatter  in  these  preliminary  longer-scalelength  experiments  was  entirely 
predictable  from  parametric  instability  theory,16  two  other  observations  are  more  surprising.  Time- 
integrated  images  of  second-harmonic  emission,  taken  at  90°  to  the  incident  laser  axis,  arc  shown  in 
Fig.  4.  When  a  single  high-intensity  beam  is  focused  on  a  carbon  target,  as  in  the  left  hand  image,  a 
localized  region  or  copious  2w„  emission  is  observed.  When  this  same  high-intensity  beam  is  incident 
upon  the  longer  scalelength  plasma  set  up  by  the  large-focus  beam,  the  intensity  of  the  2w0  emission  is 
not  detectable,  as  in  the  right-hand  image. 

The  second  surprising  observation  is  in  the  x-ray  emission.  Time  integrated  spectra  of  1-50  keV 
x-rays  are  obtained  with  arrays  of  PIN  diodes  and  photomultiplier-scintillator  detectors.  With  the  single 
high-intensity  beam  incident  upon  a  carbon  target,  x  rays  in  the  10-50  keV  spectral  range  are  easily 
detected,  as  indicated  in  Fig.  5.  However,  with  the  same  high-intensity  beam  incident  upon  the  longer 
scalelength  plasma,  the  intensity  of  x  rays  in  this  spectral  region  falls  below  detection  thresholds  (see 
f  ig  5),  The  upper  bounds  on  the  x-ray  flux  at  50  and  20  keV  represent  reductions  from  the  single¬ 
beam  case  by  factors  of  at  least  30  and  1000.  respectively!  These  reductions  are  consistent  with  obser¬ 
vations  made  m  longer-scalelength  gas-jet  targets  at  KMS.,J 
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At  least  two  explanations  might  be  offered  to  account  for  the  observations  of  2 to0  and  energetic  x 
ray  emission.  First,  a  longer  scalelength  near  critical  density  may  directly  reduce  the  importance  of 
resonance  absorption;  this  could  explain  the  disappearance  of  both  2io0  emission  and  high-energy  \- 
rays,  and  would  indicate  that  the  experimental  conditions  are  such  that  the  relative  importance  of  vari¬ 
ous  absorption  processes  is  changing.  The  observations  may  also  be  explained  by  a  reduced  laser  inten¬ 
sity  at  the  critical  surface,  due  to  increased  inverse  Bremsstrahlung  absorption  or  Brillouin  scatter  in  the 
underdensc  plasma.  The  increased  reflectivity  through  the  lens  in  these  preliminary  experiments  is 
insufficient  to  account  for  this,  but  measurements  of  scattered  light  at  other  angles,  which  have  not  yet 
been  reduced,  may  alter  the  picture. 

III.  PLASMA  PROFILE  MEASUREMENTS 

As  mentioned  above,  measurements  of  plasma  density,  temperature,  and  velocity  profiles  are  at 
least  as  important  in  these  experiments  as  measurements  of  absorbed  energy,  scattered  light,  and  x 
rays.  Without  these  measurements,  scalelengths  in  hydrocodes  may  be  varied  to  fit  the  experimental 
observations,  and  our  confidence  in  the  ability  of  these  codes  to  extrapolate  the  physics  correctly  to 
reactor  conditions  is  not  increased. 

To  overcome  the  deficiency  of  previous  coupling  experiments.  NRL  is  developing  new  methods 
for  measurement  of  plasma  profiles.  All  of  these  new  methods  have  evolved  from  a  technique 
presented  at  the  1981  Anomalous  Absorption  Meeting,  which  allowed  the  first  visualizations  of  material 
flowlines  front  laser-irradiated  targets.1  We  use  a  special  target,  as  shown  in  Fig.  6.  which  has  Al  tracer 
dots  embedded  within  a  polystyrene  (CH)  substrate.  Here,  the  dots  are  25  in  diameter  and  are 
spaced  by  50  jim  center-to-cenler.  Regions  of  Al  flow  from  each  of  the  dois  are  easily  identified  as 
tracks  of  stronger  cntissivity  in  pinhole-camera  images  of  x-ray  emission  above  I  keV,  because  the  Al 
line  emission  in  this  spectral  region  is  much  stronger  than  either  ihe  C  or  H  coniouum  emission. 

Examples  of  flow  visualizations  obtained  in  this  wav  are  shown  in  Fie  7  These  three  images, 
obtained  with  different-sized  pinholes  on  a  single  shot,  show  tracks  of  stronger  volume  emissivity 
resulting  from  the  Al  dots.  Images  obtained  with  smaller  pinholes  provide  better  spatial  resolution, 
while  those  obtained  with  larger  pinholes  give  greater  sensitivity  and  .How  visualization  of  the  flow 
further  from  the  target.  A  more  thorough  treatment  of  the  flow  pattern >  may  be  found  in  Ref  17.  the 
one  property  relevant  to  the  present  discussion  is  the  distinctness  of  the  streamlines.  Since  these  x-ray 
images  arc  time-integrated  over  the  4-nsec  x-ray  pulsewidth  (FWIIM),  the  absence  of  smearing  of  the 
streamlines  suggests  that  we  are  in  a  nearly  steady  flow  regime  Further  experimental  evidence  for  this 
flow  condition  is  obtained  from  streaked  images  of  3u>,/2  emission,  which  imply  a  nearly  stationary  rr,/4 
surface.  Is 

If  we  arc.  in  fact,  in  a  steady  flow  regime,  we  may  use  the  flow  visualizations  in  conjunction  with 
a  measured  density  profile  to  infer  the  fluid  velocity  profile  in  the  plasma  As  indicated  in  Fig.  8.  the 
method  relies  upon  the  steady-state  mass  conservation  equation,  which  tells  us  that  the  product  of  the 
density  n  in  a  streamlubc,  the  cross-sectional  area  A  of  the  streamtube,  and  the  fluid  velocity  v  in  the 
sireamtube  is  a  constant,  C.  If  we  measure  the  density  profile  by  interferometry  or  some  other  tech¬ 
nique,  and  we  obtain  the  area  A  from  the  flow  visualizations,  we  need  only  to  evaluate  the  constant  C 
in  order  to  obtain  absolute  values  for  the  velocity.  (Relative  values  may  be  obtained  without  ('.) 
Experimentally,  we  may  obtain  C  by  imposing  a  boundary  condition  far  from  the  target  where  the  velo¬ 
city  has  reached  a  constant  value;  in  this  region,  we  equate  the  velocity  with  that  measured  far  from  the 
target  by  time-of-llight  ion  detectors. 

Dala  from  a  proof-of -principle  experiment  have  been  reduced  using  this  model,  as  shown  in  lie 
9  Ihe  density  profile  obtained  interferometrically  shows  a  l/r:  dependence  of  the  density  with  dis¬ 
tance  :  far  from  the  target  ^bottom  of  Fig  9)  (n  this  same  region,  the  flow  visualizations  show  a 
spherically  divergent  flow,  giving  I  proportional  to  .  Since  the  product  of  density  and  area  is  constant 
in  this  region,  the  model  of  Fig  8  implies  that  the  flow  velocity  is  constant  flop  of  lag  9),  it  is  in  this 


region  that  we  equate  the  velocity  to  that  measured  in  the  far  field  by  ion  detectors.  In  closer  to  the 
target,  the  density  variation  changes  and  we  begin  to  infer  a  change  of  the  flow  velocity  from  its  far- 
field  value.  Interestingly,  the  region  of  nearly  constant  velocity  extends  inward  approximately  to  the 
nj  10  surface.  This  could  explain  the  fact  that  Brillouin  scatter  from  solid  targets  often  seems  to 
emanate  from  this  lower  density  region13;  velocity  gradients  may  stabilize  the  Brillouin  process  at  higher 
density.  To  infer  the  velocity  profile  for  higher  densities  than  0.3  the  limit  of  the  3u>0  inter¬ 
ferometric  probe  under  these  conditions,  a  technique  for  measuring  the  higher  densities  closer  to  the 
target  is  needed. 

To  access  higher  densities,  we  are  developing  an  improved  x-ray  spectroscopic  method,  as  illus¬ 
trated  in  Fig.  10. 19  Mere,  a  polystyrene  target  with  a  single  embedded  A1  spot  is  used,  but  now  the  Al 
is  serving  as  a  localized  source  of  x-ray  line  emission  for  spectroscopic  diagnosis.  The  Al  line  radiation 
is  collected  by  an  imaging  spectrograph  with  the  entrance  slit  oriented  to  give  spatial  resolution  along 
the  laser  axis.  This  technique  has  three  major  advantages  over  standard  laser-plasma  spectroscopy 
(where  the  target  is  entirely  Al).  First,  emission  is  obtained  only  from  a  known  volume  within  the 
plasma,  rather  than  by  integrating  along  the  diagnostic  line-of-sight  through  regions  of  differing  density 
and  temperature;  thus,  local  measurements  are  made  directly,  without  the  need  for  deconvolution. 
Second,  effects  of  plasma  reabsorption  of  the  emitted  radiation  are  greatly  reduced,  particularly  if  a 
higher  2  spot  is  embedded  in  a  lower  Z  target;  this  simplifies  interpretation  of  measured  spectra. 
Third,  the  reduced  size  of  the  emitter  leads  to  better  spectral  resolution  because  source  broadening  of 
the  detected  radiation  is  reduced. 

To  demonstrate  this  third  point,  data  obtained  with  this  new  technique  are  compared  with  stan¬ 
dard  spectroscopic  data  in  Fig.  11.  In  Fig.  11  A,  a  set  of  pinhole  images  from  an  Al  target  show  that  a 
large  volume  of  plasma  participates  in  the  x-ray  emission,  as  one  would  expect.  In  contrast,  pinhole 
images  from  a  single-spot  target  show  a  well-confined  channel  of  Al  line  emission  within  a  larger- 
diameter  region  of  plastic  continuum  emission,  as  seen  in  Fig.  11C.  The  result  of  this  reduced  emitter 
diameter  is  that  better  spectral  resolution  is  seen  in  the  spectrum  obtained  with  a  spotted  target,  in  Fig 
1 1 D,  than  is  seen  in  the  standard  spectrum  of  Fig.  I  IB.  More  details  on  this  new  technique  may  be 
found  in  Ref.  19. 

With  this  improved  spectral  resolution,  one  expects  that  belter  comparisons  between  data  and 
spectroscopic  models  will  be  made  and  that  these  will  lead  to  improved  density  and  temperature  meas¬ 
urements.2"  A  standard  technique  for  plasma  density  measurement  relies  upon  the  ratio  between  the 
lle-like  intcrcombination  line  (Me  22P)  and  the  Me-like  resonance  line  (Me-a)21;  the  measured  varia¬ 
tion  of  this  ratio  with  distance  ;  from  the  target  is  shown  in  the  left-hand  graph  of  Fig.  12.  Collisional- 
radialivc  equilibrium  mode's  indicate  that  this  ratio  should  fall  monotonieally  with  increasing  density 
except  for  plasma  opacity  effects  very  near  to  the  target.22  For  e  ^  200  /am,  this  expected  behavior  is 
observed  and  leads  to  the  solid  portion  of  the  density  profile  curve  in  the  right-hand  graph  of  Fig.  12; 
this  yields  densities  up  to  0.5  nt.  For  z  <  200  pm,  we  observe  an  apparent  levelling  off  or  even  an 
increase  in  the  Me  2’P/lle-a  ratio.  At  this  lime,  we  believe  that  this  is  due  to  interfering  spectral  lines 
which  nearly  underlie  the  Me  2’P  line  and  increase  relative  to  it  as  the  density  increases.  In  previous 
experiments,  the  spectral  resolution  was  not  sufficient  to  distinguish  these  from  the  intcrcombination 
line;  with  our  greater  spectral  resolution,  we  are  on  the  verge  of  resolving  these  lines.  Were  we  to 
naively  apply  previous  models  to  the  apparently  saturated  line  ratio,  we  would  obtain  a  density  shelf  at 
5  x  lO  "  clcctrons/cm’  for  r  <  200  pm.  In  fact,  we  believe  the  Me  2’P/Me-a  ratio  continues  to 
decrease,  as  indicated  by  the  speculative  dashed  curve.  This  would  imply  higher  densities,  as  shown  by 
the  dashed  curve  in  the  right-hand  graph  of  Fig.  12;  densities  up  to  2/;(  should  be  accessible  with  this 
technique  We  are  currently  working  on  eliminating  the  problem  by  using  a  more  detailed  model, 
which  should  lesult  in  reliable  measurement  of  these  higher  densities. 

For  temperature  determinations,  we  rely  upon  intensity  ratios  between  lle-like  and  M-like  Ryd¬ 
berg  lines  Intensity  variations  of  three  of  the  stronger  lines  wiih  distance  from  the  target  are  shown  in 
the  left-hand  side  ot  Fig.  12,  and  temperatures  inferred  from  the  ratios  of  I  hose  lines  are  shown  in  the 


4 


right-hand  graph  of  the  same  figure.  For  z  <  200 /am,  the  temperatures  are  speculative  because  the 
inferred  temperatures  are  dependent  upon  the  inferred  densities  which  are  speculative  at  this  region,  as 
discussed  above. 

One  final  note  concerns  a  new,  untested  idea  to  improve  the  velocity  profile  measurement.  The 
method  described  earlier  in  Figs.  8  and  9  suffered  from  three  drawbacks:  (1)  it  required  the  assump¬ 
tion  of  steads  slate,  (2)  it  required  simultaneous  measurement  of  the  density  profile,  and  (3)  it  gave 
only  a  time-averaged  velocity  profile.  To  overcome  all  three  of  these  problems,  we  propose  the  use  of 
layered  tracer  dots,  as  shown  in  Fig.  13.  By  making  the  tracer  dot  from  alternating  layers  of  two 
materials  which  have  distinguishable  x-ray  emissivity,  we  should  obtain  a  direct  measure  of  the  velocity 
profile  with  time  resolution.  We  may  follow  the  front  of  each  of  these  layers  as  it  transits  the  blowoff 
plasma  by  imaging  the  x  rays  onto  the  slii  of  a  streak  camera;  the  resultant  z.  I  curve  of  the  front 
directly  measures  the  velocity  profile.  The  time  evolution  of  this  profile  is  obtained  from  the  multiple 
layers;  each  interface  gives  us  the  velocity  profile  at  one  time. 

IV.  STATUS  REPORT 

To  summarize,  studies  of  laser  coupling  in  longer-scalelength  plasmas  which  better  simulate  reac¬ 
tor  plasmas  are  underway  at  NRL,  as  indicated  in  Fig.  14.  Preliminary  observations  do  show  increased 
backscatter.  as  expected.  The  concomitant  decrease  in  second-harmonic  and  high-energy  x-ray  emis¬ 
sions  may  indicate  that  we  are  entering  a  new  physics  regime  where  longer  scalelengths  are  reducing  the 
relative  importance  of  resonance  absorption. 

Because  we  realize  that  these  experiments  are  best  used  to  benchmark  the  codes  that  are  used  for 
scaling  to  reactor-sized  plasmas,  we  believe  that  measurements  of  plasma  profiles  are  at  least  as  impor¬ 
tant  as  measurements  of  laser-target  coupling.  We  have  demonstrated  one  technique  and  proposed  a 
second  improved  technique  for  velocity  profile  measurement;  these  are  the  first  techniques  available  for 
measurement  of  the  velocity  profile,  which  is  important  in  evaluating  parametric  instability  mechan¬ 
isms.  In  addition,  we  have  demonstrated  an  improved  method  for  x-ray  spectroscopic  measurements, 
and  we  are  developing  the  necessary  improved  models  for  extraction  of  plasma  density  and  tempera¬ 
ture.  We  hope  to  nuke  more  definitive  statements  on  laser-coupling  and  profile  measurements  in 
longer-scalelength  plasmas  in  the  near  future. 
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Fig.  2  —  Novel  configuration  of  two-beam  laser  allows 
production  of  longer-scaleiength  plasmas  for  coupling 
experiments.  One  beam  with  3-5  nsec  pulselength  is 
brought  to  large  focai  diameter  (=1  mm)  to  produce 
longer  scalelengths  (t„  250-600  /am  at  nr/IOi. 

while  a  second  beam  is  tightly  focused  within  the  first 
to  produce  higher  intensities  of  interest  The  pul- 
selenglh  of  the  high-intensit>  beam  may  be  shortened 
to  minimize  its  perturbation  to  the  long-scalelength 
How, 


3-7  FOLD  BACKSCATTER  INCREASES 
AT  MID-1013  (5  ns)  AND  MID-1014  {<0.5  ns) 
WHEN  LARGE-FOCUS  BEAM  PRESENT 


Fig  3  —  Summary  of  preliminary  backscaticr  observations 
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At  =  5.4ns 

Fig  4  —  Comparison  of  time-integrated  images  of  second  harmonic  emission  for  two  modes  of 
operanon  The  laser  is  inculeni  from  (he  led  and  (he  images  are  negatives  With  a  single  high- 
imcnsily  beam,  the  image  on  (he  lefi  shows  copious  2oi0  emission  When  (he  same  high-iniensity 
beam  is  focused  into  a  longcr-scalelengih  plasma,  a-  shown  in  Fig  2.  ihe  2<u0  emission  disappears 
The  delay  of  ihe  short-pulse.  high  inicnsiiv  beam  is  5  4  nsec  afier  the  peak  of  the  longer  pulse 
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X-RAY  FLUX  (Ergs/sterad/keV) 


HIGH-ENERGY  X-RAYS: 
STANDARD  vs.  LONGER  SCALELENGTH 


I  5  —  Comparison  of  high-energy  x-ray  spectra  for  two  modes  of  operation  With  a  smele.  short-pulse,  hich-intensuv 
beam,  an  x-rav  lul  corresponding  to  Tt,  7  kcV  is  observed  This  tail  falls  below  the  detection  threshold  for  a  long- 
stale  length  target  plasma  The  delay  of  the  short-puKe.  high-intensity  beam  is  5  4  nsec  after  the  peak  of  the  longer  pulse 


F  ig  7  —  Time  integrated  images  of  x-ray  emission  obtained  lor  single  laser  shot  with  13.  28,  and  54  ^ m  diameter 
pinholes  La>cr  is  mudent  from  left,  and  images  are  negatives  Camera  filter  is  nominal!)  15  gm  thick  Be  230  J  of 
laser  energy  is  focused  to  940  gin  diameter  spot  <90  energy  content)  to  give  /,,  =  8  x  It)1-'  W/crrr  on  50  gni  thkk 
CH  target 


STEADY-STATE  MASS  CONSERVATION 

n  A  v  =  constant,  C 

C 

v  =  - 

nA 


BOUNDARY  CONDITION: 

Vqq  from  time-of-flight  detectors 

Fig.  8  —  Model  for  inferring  fluid  velocity  profile. 
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ELECTRON  DENSITY  FLOW  VELOCITY 

(cm-3)  (cm/sec) 


I*ig  10  —  lixpcrimcni.il  arrangement  for  spot  spectroscopy 


(A) 


(B) 


Fig  II  —  Qualitative  comparison  between  data  obtained  with  At  foil  target  and  data  acquired  using  CM  plastic  target  locally 
embedded  with  Al.  the  vertical  axis  for  (A)  through  (I)>  corresponds  to  the  laser  axis,  with  the  laser  incident  from  above  and 
the  target  surface  at  bottom  (A)  is  obtained  with  Kodak  3490  film,  and  (B)  through  (D)  are  on  Kodak  No-Screen  Film.  (A  I 
X-ray  pinhole  camera  images,  obtained  with  an  array  of  pinholes  (5-55  jam  diameter).  Tor  an  Al  foil  target  Laser  intensity 
/0  =  )  «  I0IJ  W/cm1,  due  to  100  J  focused  to  1 100  aim  diameter  spot  OO'm  energy  content).  (B)  Spectrum  obtained  for  Al  foil 
target  at  /„  e*  8  x  10l}  W/cnv.  with  HO  J  focused  to  600  rim  diameter  spot.  (C)  Pinhole  images  as  in  (A),  but  for  CM  target 
embedded  with  180  nm  diameter  Al  spot.  225  J  of  laser  energy  focused  to  900  rim  diameter  spot  yields  /0=8x  I2IJ  W/cm:. 
(O)  Spectrum  from  same  shot  as  (C). 
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log  I  (arb) 


ne  AND  Te  FROM  AI-SPOT  SPECTROSCOPY 


z  (^m)  2  (|im) 

I90  ~  7.3  x  1012W/cm2  dgQ~1mm  E  =  273 J  tl  =  4ns  I11SO6) 

big  12  —  Preliminary  reduction  of  spot  spectroscopy  data.  Left-hand  graph  shows  spatial  variation  of  intensities  of  He-like  a 
and  p  lines  and  of  1 1 -like  a  line,  as  well  as  of  ratio  of  He-like  intcrcombination  to  resonance  lines  Right-hand  graph  shows 
density  and  temperature  information  extracted  front  line  ratios  within  spectrum  Note  that  dotted  portions  of  density  and  tem¬ 
perature  profiles  are  based  upon  speculative  extrapolation  of  He-’P/Hea  ratio  shown  in  graph  at  left 
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STATUS  REPORT 


I.  LONGER-SCALELENGTH  EXPTS. 
UNDERWAY 

•  INCREASED  BACKSCATTER 

•  DECREASED  2  co0  &  HIGH-ENERGY 
X-RAYS 

II.  DIAGNOSTIC  DEVELOPMENT 
CONTINUING 

•  FIRST  TECHNIQUES  FOR  V  PROFILE 

•  SPOT  SPECTROSCOPY  FOR  n,  T 

Hg  14  —  Summary  of  work-lo-dalc 
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NEW  MEASUREMENT  TECHNIQUES  USING  TRACERS 
WITHIN  LASER-PRODUCED  PLASMAS 


i.  lNTRODUCT  Iuii 

Many  of  the  physical  processes  wnich  are  important  to  tne  ablative- 
coiiipression  approach  to  laser-pellet  fusion  are  intimately  related  to 
plasma  density,  temperature,  and  velocity  profiles.  Inverse  bremsstranl 
absorption  and  parametric  instabilities  are  closely  coupled  witn  underderse 
plasma  profiles.  Similarly,  energy  transport  Detween  the  cri  t  ica-!- de’nsVtT 
surrace  ana  the  ablation  surface  is  interrelated  with  overdense  plasma 
profiles.  Tnerefore,  I  aser-pl  asrna  interaction  experiments  can  better 
address  these  physics  issues  if  they  include  measurements  of  the  plasma 

profiles.  If  experiments  are  performed  in  planar  rather  than  spherical 
geometry,  tne  problem  becomes  two-dimensional  since  the  plasma  profiles  anc 
the  above-mentioned  physical  processes  are  affected  by  lateral  flow  of 
energy  out  of  the  finite  focal  spot.  Or.e  tnen  desires  measurements  of  the 
area  over  whicn  the  absorbed  energy  is  distributed  ana,  if  possible,  of 
both  axial  and  lateral  profiles. 

A  new  genre  of  diagnostic  methods  is  being  developed  at  the  Naval 
Research  Laboratory  to  address  these  and  otner  problems.  The  unique 

feature  ot  tnese  new  methods  is  the  use  of  tracer  materials  which  are 

locally  embedded  into  the  solid  target  to  be  irradiated.  Tracers  which  are 
localized  within  tne  focal  region  remain  localized  within  the  blowoff 
plasma,  due  to  the  collisional  or  fluid  nature  of  the  flow  of  material 
ablated  from  the  target.  By  imaging  characteristic  radiation  from  a  linear 
array  of  such  tracer  dots,  we  have  demonstrated  the  first  visualizations  of 
the  hydrodynamic  flow  of  material  witnin  laser  produced  plasmas,  this 
technique  and  its  applications  are  discussed  in  Sec.  II.  As  described  in 
Sec.  Ill,  improved  measurements  of  plasma  density  and  temperature  profiles 
have  also  been  demonstrated,  by  using  the  tracer  dots  as  local  sources  of 
spectroscopic  diagnostic  lines.  Comparisons  witn  these  measurements 
provide  the  most  direct  calibration  of  hydrodynamics  computer  codes, 

because  plasma  profiles  are  the  experimental  observables  most  closely 
coupled  to  mechanisms  of  energy  absorption  and  transport  in  1 aser- produced 
plasmas.  Tnerefore,  comparisons  between  experiments  and  the  Naval  Researcn 
Laboratory  hydro-codes  are  found  throughout  Secs.  II  and  III.  In  Sec.  IV, 
the  use  of  layered  tracer  spots  is  proposed  for  time-resolved  measurements 
of  fluid  velocity  profiles^  and  for  improved  measurements  of  mass  ablation 
races . 

II.  FLUID  FLOW  VISUALIZATIONS  AND  THEIR  USES 

The  simplest  use  of  tracer  dots  is  for  visualization  of  the 
hydrodynamic  flow  of  material  from  the  laser  irradiated  target,  for  which 
the  experimental  arrangement  is  as  shown  in  Fig.  1.  The  distinguishing 
feature  of  this  conf i gurat i on  is  the  target:  thick  polystyrene  [  ( CH )  x ] 

witn  a  linear  array  of  26  um  diameter,  aluminum  tracer  dots.  These  dots 
are  embedded  so  that  their  initial  surfaces  are  flush  with  that  of  tne 
surrounding  target,  and  their  thicknesses  are  chosen  to  exceed  tne  expected 
ablation  deptn  on  any  yiven  shot.  Nitn  the  tracer  spots  aligned  to  tall 
along  a  diameter  of  the  laser-i  rradiated  region,  between  10  and  60U  J  of 
Nd-laser  energy  =  1.054  M,n)  is  brought  onto  the  target  in  a  3-6  nsec 

pulse.  A  camera  witn  an  array  of  pinholes  is  used  to  obtain  time- 

M.inuv  ripi  jpprowd  [)*T*»mb**r  9.  19**‘J. 
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X-RAY  PINHOLE  CAMERA 


Fig.  1  Experimental  configuration  for  visualization  of  hydrodynamic  flow 
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integrated  images  of  x-ray  emission  at  yQ°  to  the  target  normal .  The  15  pm 
beryllium  filter  used  on  the  camera  allows  imaging  of  photons  with  energy 
hu  >  1  keV.  Since  aluminum  line  emission  is  more  intense  than  carbon  or 
hydrogen  continuum  emission  in  the  spectral  band  of  the  pinhole  camera, 
streamlines  of  material  flowing  from  the  aluminum  spots  are  identifiable  as 
tracks  of  strong  x-ray  emission  in  the  images  (see  Fig.  2). 

Flow  vi  sua 1 i zat i on  experiments  have  been  performed  at  a  variety  of 

laser  focal  conditions  with  incident  intensities  l^y  (averaged ^wi  thi n  the 
90%  energy-content  diameter  dgy)  between  lO1^  and  2  x  1014  W/cm  .  On  each 
shot,  a  series  of  laser  isointensity  patterns  is  recorded  in  an  equivalent 
focal  plane  to  that  of  the  target,  with  a  factor  of  2.2  in  intensity 
between  successive  patterns.  From  these  records,  we  may  plot  focal 
intensity  distributions,  such  as  those  shown  in  Fig.  3.  (These  one¬ 

dimensional  profiles  represent  azimuthal  averages  of  the  real  two- 
dimensional  profiles  of  the  laser;  these  are  used  as  input  to  the  hydro- 

codes.)  Note  that  the  focal  distribution  near  best  focus  of  our  f/6  lens 
is  qualitatively  different  from  that  in  the  quasi-near  field.  The  most 

obvious  difference  is  near  the  center  of  the  distribution,  where  large 
focal  spots  tend  to  exhibit  local  minima  in  intensity.  Equally  important, 
however,  are  differences  in  the  wings  of  the  distributions,  as  will  be 
discussed  in  Sec.  II. B.  It  should  also  be  noted  that  intensities  I^q 
(averaged  over  50%  energy  content  diameter  d^)  and  I  ^  (peak  intensity) 
tend  to  be  higher  than  IgQ  by  factors  of  2  and  4,  respectively. 

While  we  have  not  conclusively  demonstrated  that  the  presence  of  the 
tracers  does  not  perturb  the  visualized  flow,  there  is  some  evidence  that 
the  perturbation  should  be  minimal.  First,  the  fully  ionized  mass 
densities  for  A1  and  ( CH )  differ  by  only  12%  at  a  given  electron  density 
in  the  blowoff  plasma.  Moreover,  the  average  velocity  of  ablated  ions  as 
determined  by  Faraday  cups  is  the  same  for  the  two  materials. 

In  the  following  sections,  the  various  physics  issues  that  may  be 
addressed  by  flow  visual i zation  methods  are  discussed.  The  shape  of  the 
visualized  flows,  as  shown  in  Sec.  II. A,  is  related  to  the  outward 
hydrodynamic  flow  of  plasma,  and  provides  a  valuable  check  for  the 
hydrodynamics  codes.  In  Sec.  II. B,  we  present  the  observed  sizes  of  the 
fluid  sources  under  various  irradiation  conditions;  this  should  indicate 
the  degree  to  which  absorbed  laser  energy  spreads  laterally  as  it  flows 
inward  from  the  absorption  region  to  the  ablation  region  near  the  edges  of 
the  focal  spot.  Lateral  energy  transport  can  also  be  important  near  the 
center  of  the  focal  spot  if  the  irradiation  there  is  nonuni  form.  ,b  Flow 
visualizations  obtained  with  intentionally  perturbed  laser  beams 
qua  1 i tati vely  show  this  effect,  as  described  in  Sec.  II. C.  Another 
proposed  application  for  these  Mow  visualizations,  outlined  in  the  final 
section,  is  for  diagnosis  of  the  Kaylei gh-Taylor  instability/ 

A.  Shape  of  the  Visualized  Flows 

While  the  si ze  of  the  observed  fluid  source  varies  with  laser 
intensity  ano  spot  size,  the  shape  of  the  flow  pattern  seems  to  be 
remarxably  constant.  As  shown  in  Tfg.  2,  streamlines  curve  away  from  the 
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target  normal  at  a  rate  which  depends  upon  r/r  ,  where  r  is  the 
displacement  of  a  given  streamline  from  the  center  of  the  fluid  source  and 
rQ  is  the  radius  of  the  fluid  source.  Each  streamline  attains  a  constant 
direction  within  a  distance  z  -  r  from  the  target  surface;  beyond  z  -  c  , 
all  streamlines  are  following  straight-line  paths  corresponding  to  a  far- 
field ,  spherically  diverging  flow. 

The  shape  of  the  visualized  flows  is  in  good  agreement  with  NRL's  two- 
dimensional  hydrodynamics  codes.  As  shown  in  Fig.  4,  the  measured 
streamline  angles  e  relative  to  the  target  normal  are  well  predicted  by  a 
cylindrical  (r  and  z)  hydrodynamics  simulation  of  the  experiment.  The 
computer  model,  described  elsewhere  in  these  proceedings,  uses  a  sliding- 
grid  Eulerian  mesh  with  fl ux-corrected  transport,  classical  T^'^  thermal 
conduction,  and  inverse  bremsstrahlung  absorption.  The  calculation  uses 
the  measured  laser  pulse  duration  and  an  azimuthal  average  of  the  measured 
laser  focal  distribution  (such  as  those  in  Fig.  3).  The  slight  asymmetry 
of  the  experimental  data  about  r  =  0  is  a  gentle  reminder  that  real  laser 
profiles  are  not  perfectly  symmetric,  but  prohibitively  expensive,  three- 
dimensional  codes  would  be  required  to  accurately  simulate  such  focal 
nonuniformities.  While  Fig.  4  shows  agreement  at  one  distance  z  from  the 
target,  sets  of  such  plots  at  varying  z  convince  us  that  the  code  is 
correctly  calculating  the  shape  of  the  flow. 

The  fluid  nature  of  the  flow  is  not  surprising,  since  the  plasma  radii 
are  much  greater  than  the  ion-ion  collision  mean  free  path.  The  fact  that 
time-integrated  flow  visualizations  compare  well  with  snapshots  of  the 
code,  as  in  Fig.  4,  suggests  that  we  may  also  have  a  nearly  steady  flow; 
this  is  not  entirely  unreasonable,  since  the  laser  pulse  length  exceeds  the 
acoustic  transit  time  across  a  plasma  radius.  Experimental  evidence  that 
the  flow  is  nearly  steady,  at  least  during  the  period  of  maximal  x-ray 
emission,  is  obtained  from  streak  photography  of  images  of  3/2u>0  emission 
from  the  n  /4  surface.  The  position  along  the  laser  axis  of  this  surface 
appears  relatively  stationary  for  the  duration  of  the  3/Zw0  emission,  a 
period  of  up  to  3  nsec  near  the  time  of  peak  laser  intensity. 

The  shape  of  the  flow  affects  far-field  particle  diagnostics  such  as 
ion  collectors  or  plasma  cal orimeters.  First,  the  angular  distribution  of 
ablated  mass  observed  by  these  diagnostics  may  be  determined  by  the 
col  1 isional  flow  near  the  target,  which  imparts  a  fluid  velocity  direction 
to  each  fluid  element.  As  the  material  drifts  through  the  collisionless 
region  further  from  the  target,  some  lateral  spreading  of  material  results 
from  random  thermal  motions;  since  directed  ablation  velocities  are 
generally  much  greater  than  ion  thermal  velocities,  this  should  not  greatly 
increase  the  angular  range  set  by  the  fluid  flow.  Indeed,  charge  collector 
and  plasma  calorimeter  measurements  indicate  that  half  of  the  ablated  mass 
appears  in  a  cone  about  the  target  normal  with  a  40°  half -angle,  this  is 
consistent  with  the  angular  spread  inferred  from  flow  visual  izations. 

Ion  detectors  are  further  affected  because  the  flow,  except  for  any 
smearing  due  to  thermal  expansion  in  the  collisionless  region,  maps 
different  regions  of  target  surface  into  different  diagnostic  solid 
angles.  A  single  ion  detector  at  a  specific  angle  far  from  the  target  then 
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50° 


VECTOR  FLUID  FLOW  AT  Z  =  85  /xm 
2-D  CODE  vs  EXPERIMENT 


RADIUS  (/xm) 

g.  4  Comparison  between  measured  and  calculated  streamline  angles  6 

relative  to  the  target  normal  as  functions  of  the  radial  displacement 
from  the  center  of  the  flow.  Values  are  taken  at  one  distance  z  = 

85  um  from  the  target  surface.  The  calculation  uses  measured 
experimental  conditions:  In0  =  2  x  1012  W/cm  ,  dgU  =  650  urn,  laser 

energy  is  35  J. 


samples  only  a  portion  of  the  irradiated  region.  To  measure  parameters 
averaged  over  the  focal  spot,  one  should  use  an  array  of  detectors  at  a 
range  of  angles.  The  error  that  may  result  if  too  few  detectors  are  used 
depends  upon  the  uniformity  of  the  laser  illumination  at  the  target,  the 
ratio  of  the  directed  ablation  velocity  to  the  ion  thermal  velocities,  and 
the  solid  angle  subtended  by  the  detectors  used. 

B.  Sizes  of  the  Visualized  Flows 

The  sizes  of  the  visualized  flows  are  well-defined  and  useful 

observables.  From  the  flow  visualizations,  we  define  the  fluid  source  size 
to  be  the  diameter  of  the  region  at  the  target  surface  from  which 
streamlines  are  seen  to  emanate.  This  quantity  does  have  physical 

signi f icance ;  as  discussed  below,  it  most  likely  reflects  the  diameter 

within  which  the  plasma  temperature  is  high  enough  for  excitation  of  A1  K- 
lines.  It  is  also  a  well-defined  quantity,  in  the  sense  that  it  is  not 

extremely  sensitive  to  the  size  of  pinhole  used  to  detect  the  x-rays.  As 
seen  in  Fig.  2,  the  same  number  of  streamlines  is  observed  in  the  three 

largest  pinhole  images,  despite  a  factor  of  17  variation  in  pinhole 

collection  efficiency.  This  may  indicate  the  rapidity  with  which  the 
temperature  is  falling  at  the  edge  of  the  x-ray  images.  The  advantage  of 
this  observable  over  the  size  of  standard  pinhole  images  lies  mainly  in  the 
ease  of  observation.  The  size  of  standard  images  is  usually  sensitive  to 
pinhole  size.  In  order  to  define  a  plasma  size,  one  must  densitometer  the 
images,  deconvolve  the  pinhole  size,  and  invert  the  data  to  obtain 

isointensity  contours  of  x-ray  emission. 

The  observed  fluid  source  diameter  dQ  varies  with  laser  intensity  and 
spot  size.  As  indicated  in  Fig.  5,  dQ  increases  with  laser  energy  at  fixed 
spot  size  and  increases  with  spot  size  at  fixed  laser  energy.  One 
interesting  observation  is  that  the  sensitivity  of  d0  to  laser  energy  is 

greater  for  the  smaller  focal  spot  sizes.  Several  possible  explanations 
exist  for  such  an  observation:  there  is  the  potential  effect  of  higher 
lateral  thermal  conductivity,  which  could  result  from  higher  plasma 

temperatures  at  the  higher  intensities,  as  well  as  the  possibility  of 

increased  lateral  flow  due  to  fast  electrons.  A  closer  look  at  the  shapes 
of  the  wings  of  the  laser  focal  distributions  in  Fig.  3  suggests  another 
explanation.  By  locating  the  position  within  the  focal  distributions  at 
which  a  given  intensity  occurs  and  plotting  the  variation  of  that  position 
with  increasing  laser  energy,  we  generate  curves  such  as  shown  in  Fig.  5. 
Within  experimental  uncertainty  (approximately  a  factor  of  2),  the 

intensity  for  both  curves  shown  in  Fig.  5  is  2  x  1011  W/cm  .  The 

correlation  with  the  data  points  suggests  that  the  edge  of  the  fluid  source 
consistently  occurs  at  that  particular  intensity  in  the  wings  of  the  laser 
focal  distribution.  This  could  be  explained  by  the  existence  of  a 

threshold  intensity  for  heating  to  a  temperature  at  which  the  A1  K-lines 
are  excited;  this  heating  could  be  directly  due  to  laser  energy  deposition 

or  could  be  due  to  laser  heating  to  a  temperature  at  which  the  thermal 

conductivity  allows  lateral  transport  of  energy  from  higher  intensity 

regions.  It  is  not  likely  that  intensity  thresholds  for  ionization  of  the 
target  play  any  role  in  determining  the  fluid  source  size,  since  the  two- 
dimensional  hydrodynamics  codes  calculate  the  flow  pattern  correctly  (see 


VISUALIZED  FLUID  SOURCE  DIAMETER  (fj. m) 


Fig.  5  Variation  of  fluid  source  diameter  (data  points)  with  laser  energy  for 
the  two  laser  focal  conditions  shown  in  Fig.  3.  Curves  show  the 
location  of  2  x  1011  W/cmZ  intensity  in  the  wings  of  the  focal  profile 
as  function  of  laser  energy  for  these  same  two  focal  conditions. 
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Fig.  4)  while  assuming  initial  ionization  everywhere  on  the  target. 

If  the  fluid  source  size  is  entirely  determined  by  the  shape  of  fne 
wings  of  the  laser  focal  distribution,  one  might  infer  that  the  effects  of 
lateral  energy  transport  at  the  edges  of  the  focal  spot  are  minimal  or  that 
the  lateral  energy  flow  does  not  result  in  a  high  plasma  temperature  at  the 
focal  periphery.  This  empirical  observation  may  be  purely  coincidental, 
however.  furtner  comparisons  should  be  made  between  the  experiment  and 
code  at  a  variety  of  intensities  and  focal  spot  sizes  before  such  a 
conclusion  can  be  drawn. 

C.  Flow  Visualizations  with  Intentionally  Perturbed  Laser  Intensity 

Distributions 

Lateral  eneryy  transport  between  absorption  and  ablation  surfaces  can 
also  be  addressed  by  flow  visualizations  in  the  presence  of  intentionally 
imposed  laser  nonuniformities.  This  diagnostic  complements  target 

acceleration  experiments,  where  velocity  modulations  Av/v  have  been 
measured  as  a  function  of  imposed  intensity  modulations  Al/I;^’3’^  with 
the  tracer  technique  one  hopes  to  access  the  modulation  in  ablation 
pressure  aP/P  which  couples  Av/v  to  Al/I. 

The  idea  is  qualitatively  illustrated  in  Fig.  6,  where  two  flow 
visualizations  are  compared.  The  upper  image  is  obtained  with  a  deeper  and 
wider  intensity  minimum  imposed  into  the  center  of  the  focal 
distribution.  Under  these  conditions,  target  acceleration  results3 
indicate  little  smoothing  of  the  intensity  modulation  by  lateral  thermal 
transport  between  absorption  and  ablation  surfaces.  Indeed,  the  flow 
visualization  shows  two  seemingly  independent  fluid  sources,  with  no 
detectable  ablation  from  the  central  region  where  the  intensity  is  low.  In 
contrast,  the  lower  image  is  obtained  at  higher  average  irradiance,  with  a 
shallower  and  narrower  intensity  minimum  in  the  center  of  the  beam;  under 
these  conditions,  smoothing  of  the  intensity  modulation  is  observed  in  the 
target  acceleration  studies.  Flow  is  now  seen  from  the  center  of  the 
focal  region,  although  compression  of  this  reyion  due  to  inward  flow  of 
surrounding  streamlines  indicates  that  there  is  residual  pressure 
modulation.  More  quantitative  analysis  of  pressure  modulations  requires 
comparisons  of  flow  visualizations  with  the  code  or  development  of  new 

techniques  discussed  in  Sec.  III.C. 

D.  Detection  of  the  Rayleigh-Taylor  Instability 

Another  potential  use  of  the  flow  visualization  method  is  for 
observation  of  hydrodynamic  instability  of  accelerating  targets.'  Two- 
dimensional  hydrodynamic  simulations  of  foils  unstable  to  Rayleigh-Taylor 
growth  have  shown  that  nonuniform  ablation  from  the  target  surface 
accompanies  the  nonlinear  state  of  the  instability.  In  the  extreme  limit 
of  bubble-and-spike  formation,  these  computer  calculations  show  that  all  of 
tne  plasma  flowiny  from  the  target  can  be  emanating  from  the  tips  of  the 

spikes,  while  none  flows  from  the  target  surface  in  the  regions  of  the 

bubbles.  Preliminary  flow  visualization  experiments  have  been  performed  to 
1oo<  for  this  nonuniform  ablation,  using  targets  with  periodically 
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perturbed  areal  densities  to  seed  instability  growth  at  given 
wavelengths.  J  The  use  of  these  perturbed  targets  also  allows  one  to 
selectively  deposit  the  tracers  at  locations  which  are  expected  to  develop 
into  the  spikes  of  the  unstable  target;  one  can  then  look  for  the  tracer 
material  to  appear  everywhere  in  the  blowoff,  signifying  the  expected 
nonuniformity  in  ablation.  These  preliminary  experiments  have  produced 
flow  patterns  which  show  qualitative^  deviations  from  the  usual  patterns, 
but  they  are  not  yet  fully  explained. 


III.  SPOT  SPECTROSCOPY:  IMPROVED  DENSITY  AND  TEMPERATURE  MEASUREMENTS 

A  natural  extension  of  the  tracer  technique  is  into  the  area  of 
spectroscopic  plasma  diagnosis.  As  pointed  out  in  Sec.  II,  it  is  the 

strong  line  emission  from  the  ablated  tracer  material  which  allows 
vi sual i zation  of  25  pm  diameter  streamlines  within  1  mm  diameter  plasmas. 
These  emitted  lines  can  also  serve  for  spectroscopic  measurements  of  plasma 
density  and  temperature.  The  new  technique,  spot  spectroscopy,  has  three 
major  advantages  over  standard  spectroscopic  methods.  Standard 
spectroscopic  measurements  are  integrated  along  a  diagnostic  line  of  sight 
through  an  inhomogeneous  plasma.  One  could  use  strai ghtforward  inversion 
techniques  to  obtain  spatially  resolved  emissivities  from  the  chord- 
integrated  data  if  the  plasmas  were  optically  thin,  4  but  effects  of  plasma 
opacity  on  the  radiation  exiting  the  plasma  are  non-negl i gi ble  at  the 
higher  densities  of  interest.  A  more  complicated  inversion  technique 

might  also  account  for  the  spatially  varying  opacity;  however,  the  more 
typical  .  procedure  has  been  to  predict  the  plasma  profile  with  a 
hydrodynamics  computer  code  and  to  compare  the  experimentally  observed 
spectrum  with  those  calculated  for  chord  integrations  across  the 
inhomogeneous  plasma.  In  that  case,  one  is  dependent  upon  the 

correctness  of  the  hydrodynamics  model  as  well  as  the  atomic  physics  model, 
and  one  would  obviously  prefer  to  make  more  direct  determinations  of 
density  and  temperature.  Two  of  the  three  advantages  of  spot  spectroscopy 
alleviate  these  problems  of  standard  spectroscopy:  1)  tne  localized  source 
of  spectroscopic  lines  should  eliminate  chord  integration  across  an 
inhomogeneous  plasma,  and  2)  the  smaller  source  size  also  reduces 
complicating  effects  due  to  plasma  opacity.  The  third  advantage  of  spot 
spectroscopy  is  that  better  spectral  resolution  can  be  realized,  since 
spectral  resolution  with  standard  x-ray  spectroscopy  is  usually  limited  by 
source  broadening  due  to  finite  plasma  size. 

A.  The  Experiment 

Except  for  the  presence  of  the  tracer  in  the  target,  the  experimental 
arrangement  for  the  new  technique  is  much  the  same  as  for  standard  x-ray 
spectroscopy.  As  shown  in  the  illustration  of  Fig.  7,  the  target  has  a 
single  aluminum  tracer  spot  embedded  in  the  center  of  the  focal  region;  the 
variety  of  laser  irradiation  conditions  used  are  the  same  as  those  used  for 
the  flew  visual ization  technique  (see  Sec.  1 1 ) .  Aluminum  line  radiation  in 
the  5-3  A  spectral  region  is  collected  by  an  x-ray  crystal  spectrograph 
wi tn  a  viewing  axis  parallel  to  the  target  surface.  The  spectrograph  slit 
is  oriented  to  yield  spatial  resolution  in  •  ne  direction  along  the  laser 


12 


axis;  resolution  in  the  other  two  spatial  dimensions  results  frc  :r~ 
knowledge  of  tne  position  of  tr.e  localized  tracer  witnir  tne  focal  set:. 
The  spectra  are  recorded  on  Kocax  No-Screen  film,  for  whicti  a  film  m;cel -- 
based  upon  an  aDsolute  cal ibratici/u  is  available.  A  beryllium  filter  :: 
nominally  IS  um  tnickness  protects  tne  film  from  visible  exposure. 

More  specific  conditions  depend  upon  tne  goals  of  tne  particular 
experiment.  As  with  standard  x-ray  spectroscopy,  one  may  vary  tne  crystal 
type  and  geometry,  the  spectrograph  slit  wictn  and  magni  fi  cat  ion ,  arc  tre 
si it-to-target  distance  in  order  to  acnieve  tne  optimal  tradeoff  between 

spectral  resolution,  spatial  resolution  along  tne  laser  axis,  arc 
sensitivity.  An  additional  variable  available  with  spot  spectroscopy  is 
the  diameter  of  tne  embedded  tracer;  this  allows  one  to  specify  spatial 
resolution  in  tne  direction  lateral  to  the  laser  axis,  but  it  also  impacts 
achievable  spectral  resolution  ana  sensitivity.  Tne  specific  choice  of 
parameters  for  tne  present  work  are  the  same  as  those  listed  in  Ref.  2. 

Data  obtained  witn  tne  new  tecnniMue  do  show  the  expected  improve:"ert 
in  spectral  resolution.  'witn  an  aluminum  foil  target,  as  shown  in  rig.  7, 
the  spectral  widths  of  tne  lines  near  tne  target  surface  are  quite  orcac; 

the  widths  exceed  tnose  expected  for  other  broadening  mechanisms,  arc 
reflect  source  broadening  due  to  the  large  A1  plasma  extent  transverse  to 
the  laser  axis.  Also  shown  in  Fig.  7  is  a  spectrum  obtained  using  a  (C-.) 
target  with  an  embedded  A1  act.  Line  widths  near  the  target  surface  a^e 
much  narrower,  and  sets  of  lines  are  clearly  resolved  that  are  not  resolved 
i n  the  A]  foi 1  spectrum. 

B.  Density  and  Temperature  Profiles 

Several  features  of  tne  experimental  spectra  may  be  used  in 

determinations  of  plasma  density  and  temperature.  In  the  present  worx,  we 
extract  temperatures  from  intensity  ratios  of  various  heliurn-lixe  (A1  XIII 
ionic  lines  to  nydrogen- 1 i ke  (A1  XIII)  ionic  lines.  These  ratios  are 

compared  with  calculated  ratios  from  a  cylindrical  plasma/atomic  phys’ts 
model  of  a  finite,  homogeneous  plasma  with  diameter  equal  to  the  measured 
tracer  plasma  diameter/1  At  least  four  methods  exist  for  density 

determination:  1)  tne  often-used  technique  based  uppp  the  ratio  between 

the  helium-like  resonance  and  intercombi nation  lines,  ^  2)  a  method  based 
upon  the  ratio  between  two  helium-like  satellite  lines  (ls2s  °S  -*■  2s2g  °P 
and  ls2p  r  2p  °P ) ,  3)  Stark  profile  analysis  of  higher  Rydberg 
members,  and  4)  profile  analysis  of  opaci ty-broadened  lines  near  the  target 
surface.  To  date,  we  have  exploited  the  first  two  methods,  though  we  have 
discovered  a  need  to  modify  the  first  method  when  comparing  witn  the 

experiment;  preliminary  work  is  underway  on  the  two  remaining  techniques. 

A  problem  with  previous  models  for  density  measurement  is  identified 
during  data  reduction.  Variations  of  intensity  of  various  spectral 
features  with  distance  from  the  taryet  are  plotted  in  Fig.  8.  Plots  are 

shown  of  the  three  strongest  helium-like  and  hydrogen- 1 i ke  lines  used  in 
determining  temperature,  but  tne  more  important  curve  in  the  figure  snows 
tne  variation  of  the  ratio  of  tne  nelium-like  intercoinoinaticn  (He0?)  and 

resonance  (He  a)  lines.  As  mentioned  above,  this  ratio  is  routinely  used 
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for  density  measurement.  Col lisional -radiative  equilibrium  (CKE)  models 
indicate  that  this  ratio  should  fall  monotonical ly  with  increasing  density 
except  for  plasma  opacity  effects  very  near  to  the  target.  1  For  z  _> 
200  ym,  this  expected  behavior  is  observed  and  leads  directly  to  measured 
densities  up  to  0.5  nc.  For  z  <  200  pm,  however,  we  observe  an  apparent 

levelling  off  or  even  an  increase  in  the  line  ratio;  were  we  to  naively 

apply  previous  models  to  this  data,  we  would  infer  a  density  shelf  of  n  = 
0.5  nc  for  z  <_  200  pm.  With  our  improved  spectral  resolution,  we  can 
almost  resolve  satellite  lines  which  nearly  underlie  the  long-wavelength 

side  of  the  intercombination  line,  and  interfere  with  the  intercombination 
intensity  measurement.  These  lines,  identified  as  the  lithium-like  s,t 

1  ri<-2o,2e  .  1  /  'ir  On  In  \  2d  2c  „  1 o~lr 


satellites  [Is 
are  seen  in  the 
intercombination 
calculated  intercombination 


2s^S,/2  -  ls(2s2piP)‘:Pa/2  and  ls‘2s  cSl/z  ♦  ls(2s2pAP)*P1/2] 
CHE  calculations  to  increase  in  intensity  relative  to  tne 
line  as  the  density  increases.  When  the  sum  of  the 

intensities  is  compared  with  the 


and  satellite  intensities  is  compared  with 
data,  injured  densities  continue  to  increase  for  z  <  200  pm,  as  shown  in 
Fig.  9.  3  The  densities  measured  in  this  way  are  then  in  reasonable 
agreement  with  those  measured  for  z  £  100  pm  using  technique  (2)  of  the 
previous  paragraph;  for  z  >  100  pm,  no  such  comparison  is  possible,  because 
of  insufficient  intensity  in  the  helium-like  satellite  lines. 


The  density  and  temperature  profiles  measured  in  this  way  are  in  fair 
agreement  with  those  calculated  using  the  two-dimensional  hydrodynamics 
code.  Error  bars  on  measured  densities  are  shown  only  for  z  <  100  pm, 
where  they  reflect  the  range  of  values  from  the  two  methods  avaTlable  in 
that  region;  error  bars  for  z  >  100  pm  are  not  necessarily  any  smaller,  but 
are  more  difficult  to  estimate.  Even  without  these  experimental 
uncertainties  at  the  lower  densities,  agreement  in  Fig.  9  is  seen  to  be 
good.  Agreement  is  not  quite  as  good  between  the  calculated  and  measured 
temperature  profiles  in  Fig.  9.  Temperatures  near  the  target  (z  <_  100  pm) 
and  beyond  the  point  of  peak  temperature  (z  >  400  pm)  are  in  good 
agreement,  but  the  details  of  the  intermediate  region  reveal  discrepancies 
as  large  as  35%;  further  comparisons  between  experiments  and  calculations 
will  be  required  to  resolve  these  differences. 


To  demonstrate  that  this  agreement  is  sensitive  to  details  of  the 
physics  contained  in  the  code,  we  also  show  profiles  in  Fig.  9  from 
calculations  without  inverse  bremsstrahlung  absorption  included.  Using  a 
simple  dump  of  laser  energy  at  the  critical  surface,  the  underdense  density 
and  temperature  profiles  are  greatly  changed  and  they  no  longer  agree  with 
the  measured  profiles.  It  is  interesting  to  note  that  the  overden se 
profiles  obtained  with  inverse  bremsstrahlung  are  very  similar  to  those 
obtained  with  the  critical  dump,  indicating  that  overdense  profiles  are  not 
sensitive  to  details  of  the  underdense  absorption  physics. 

Comparisons  such  as  these  provide  the  most  direct  calibration  of 
hydrodynamics  codes.  Density,  temperature,  and  velocity  profiles  are  the 
experimental  observables  most  closely  coupled  to  mechanisms  of  energy 
absorption  and  transport  in  laser-produced  plasmas.  The  ability  to 
reproduce  measured  profiles  is,  therefore,  a  severe  test  of  the  ability  of 
the  code  to  treat  these  physical  processes. 
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DETAILS  OF  ABSORPTION  PHYSICS  GREATLY  AFFECT  DENSITY  AND 

TEMPERATURE  PROFILES 


O  65Mm  A I  DOT  (11506) 
□  115pm  Af  DOT  (112441 


lN  »  7  *  IQ”  W/em* 
dN  3  1  mm 
rt  “  4  n*«c 


2-0  CODE.  INVERSE  BREMSSTRAHLUNG 
2-0  CODE.  CRITICAL  DUMP 


Fig.  9  Data  points  show  density  and  temperature  profiles  from  spot 
spectroscopy  with  65  pm  (o)  and  115  pm  (□)  diameter  A1  dots. 
Experimental  conditions  are  IgQ  *  7  x  10*Z  W/cmz  within  d^0  *  1  mm. 

2-D  code  calculations  for  the  same  conditions  are  shown  by  the  solid 
curve.  The  sensitivity  of  the  agreement  between  code  and  experiment 
to  details  of  the  code  physics  is  shown  by  the  dotted  curve,  generated 
with  the  2-D  code  altered  to  neglect  inverse  bremsstrahlung  and  to 
dump  all  energy  at  critical  density. 
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C.  Caveats  and  Future  Work 

At  least  two  sets  of  caveats  regardiny  this  work  should  be  pointed 
out.  First,  there  are  three  assumptions  which  we  make  when  reducing  the 
data  regarding  the  spatial  qualities  of  the  spectroscopic  information.  As 
pointed  out  in  Sec.  II. 8,  the  plasma/atomic  physics  model  with  which  the 
data  are  compared  simulates  a  homogeneous  cylinder  of  plasma;  the  implicit 
assumption  here  is  that  the  A1  tracer  plasma  is  nearly  homogeneous  in  the 
radial  direction.  We  also  assume  that  the  A1  tracer  plasma  equilibrates 
with  the  surrounding  CH  plasma  if  tracer  spot  diameters  are  kept 
sufficiently  small.  This  allows  us  to  compare  the  spectroscopic 
measurements  of  density  and  temperature  (which  reflect  conditions  within 
the  A1  tracer  plasma)  with  the  results  of  the  hydro-code  (which  simulates  a 
CH  target  with  no  A1  tracer  )  in  Fig.  y.  These  first  two  assumptions  seem 
to  be  supported  by  the  measured  profiles  in  Fig.  9,  where  reduced  data  from 
two  shots  with  different  tracer  diameters  are  shown  to  agree.  A  third 
assumption  is  that  the  spectrum  emitted  at  a  given  distance  z  from  the 
target  is  determined  by  the  density  and  temperature  at  that  position,  with 
negligible  effects  due  to  photon  pumping  from  higher  density  tracer  plasma 
closer  to  the  target.  This  assumption  will  be  tested  using  the  above- 
mentioned  plasma/atomic  physics  model. 

A  second  set  of  caveats  concerns  the  temporal  qualities  of  the 
spectroscopic  information.  The  present  data  is  time-integrated;  in  making 
the  comparison  in  Fig.  9  with  an  instantaneous  profile  from  the  hydro-code 
(at  a  time  near  the  peak  of  the  laser  pulse),  we  are  implicitly  assuming 
that  the  hydrodynamics  is  nearly  in  steady  state.  It  is  appropriate  to 
recall  the  experimental  inferences  discussed  in  Sec.  II. A,  which  support 
this  assumption.  By  comparing  the  measurements  with  a  CHE  model,  we  are 
also  making  assumptions  regarding  the  steadiness  of  the  atomic  physics.  As 
the  plasma  flows  through  a  given  volume  at  a  distance  z  from  the  target,  we 
assume  that  its  residence  time  in  that  volume  is  sufficiently  long  that 
states  of  ionization  and  excitation  approach  the  equilibrium  conditions 
appropriate  for  n(z)  and  T(z).  This  assumption  is  weakest  far  from  the 
target  surface,  where  relaxation  times  increase  with  decreasing  plasma 
density.  In  this  case,  the  state  of  plasma  ionization  and  excitation 
becomes  dependent  upon  conditions  that  exist  upstream  in  the  hydrodynamic 
flow,  and  the  temperature  determined  by  spot  spectroscopy  may  reflect  an 
ionization  temperature  rather  than  the  local  electron  temperature. 

In  future  work,  time-resolved  measurements  of  the  profiles  should 
provide  a  more  direct  answer  to  this  second  set  of  questions.  These 
measurements,  as  pointed  out  in  Ref.  2,  might  be  obtained  by  replacing  the 
film  with  an  x-ray  streak  camera.  An  alternative  approach  is  to  bury  a 
tracer  spot  of  thickness  t  at  a  depth  d  within  the  target.  By  varying  the 
depth  d  and  thickness  t  of  the  tracer  on  a  shot -by-shot  basis,  one  may  vary 
the  time  interval  during  which  the  tracer  emission  occurs. 

Future  work  may  also  allow  us  to  satisfy  the  need  for  measurements  of 
lateral  profiles.  As  mentioned  in  the  introduction,  the  use  of  finite 
focal  spots  on  planar  targets  causes  the  problem  to  become  two-dimensional, 
and  simultaneous  measurement  of  axial  and  lateral  profiles  is  desirable. 
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By  changing  the  orientation  of  the  spectrograph,  one  can  take  full 
advantage  of  its  two-dimensional  imaging  capability.  Then,  using  a  target 
identical  to  those  used  in  flow  visual izations  (Fig.  1),  one  obtains  an 
image  of  the  visualized  flow  in  each  spectral  line  detected  by  the 
spectrograph.  In  principle,  this  technique  should  allow  determination  of 
the  desired  profiles.  Quantitative  studies  of  lateral  variations  in  plasma 
pressure,  as  discussed  in  Sec.  II. C,  may  then  be  feasible. 


IV.  LAYERED  TRACERS:  VELOCITY  PROFILE  AND  ABLATION  RATE  MEASUREMENTS 

Further  applications  of  tracers  become  available  with  the  use  of 
layered  tracer  spots.  Layered -target  methods,  suggested  originally  by  G. 
Qahlbacka  and  developed  by  NRL,  facilitated  studies  of  axial  transport 
through  measurements  of  mass  ablation.  Here,  we  propose  the  marriage  of 
layered-target  techniques  with  our  new  tracer  methods  for  measurements  of 
fluid  velocity  profiles  and  improved  measurements  of  mass  ablation  rates. 

Use  of  a  layered  tracer  spot  should  allow  direct  and  time-resolved 
measurement  of  fluid  velocity  profiles.  Before  tracers,  velocity  profiles 
of  the  fluid  were  not  measured,  despite  their  importance  for  interpreting 
stimulated  Brillouin  scattering  observations.  In  Ref.  3,  we  demonstrated 
an  indirect  method  for  inferring  time-averaged  profiles,  but  tracers  also 
allow  a  direct  and  time-resolved  measurement.  As  shown  in  Fig.  10,  one  can 
fabricate  a  tracer  spot  from  alternating  layers  of  two  materials  with 
distinguishable  x-ray  emissi vities.  As  each  successive  layer  ablates  from 
the  target,  the  progress  through  the  blowoff  plasma  of  the  interface 
between  layers  may  be  observed  with  an  x-ray  streak  camera;  the  camera  slit 
is  oriented  to  give  spatial  resolution  along  the  laser  axis.  The  z(t) 
curves  generated  by  the  multiple  interfaces  may  be  inverted  to  give  the 
velocity  profiles  u(z)  with  time  resolution. 

The  above  data  also  yields  the  ablation  rate  through  each  layer  of 
the  tracer  spot.  This  is  an  improvement  over  the  original  layered-target 
measurement  of  ablation  rate,  5  because  localization  of  the  layers  within 
the  tracer  spot  alleviates  complications  due  to  variations  in  laser 
intensity  across  the  irradiated  region. 


V.  SUMMARY 

With  tracer  methods  still  in  their  infancy,  we  see  several  new 
techniques  which  are  capable  of  addressing  a  variety  of  interesting  and 
important  physics  issues.  The  measurements  of  density  and  temperature 
profiles  shown  in  Sec.  II  and  the  proposed  measurement  of  velocity  profile 
discussed  in  Sec.  IV  are  an  important  adjunct  to  experiments  related  to 
laser  absorption  in  the  underdense  plasma.  If  overdense  profiles  can  be 
measured  with  sufficient  precision,  issues  of  axial  transport  of  energy 
inward  from  the  absorption  region  may  be  treated;  the  applicability  of  the 
time-resolved  ablation  rate  measurement  to  this  problem  is  also  clear.  We 
have  already  begun  to  apply  the  flow  visual ization  technique  to  questions 
of  lateral  energy  spread  due  to  finite  focal  spot  size  or  due  to  laser 


nonuniformities,  and  proposed  spectroscopic  methods  can  make  this  more 
quantitative.  Experiments  are  also  underway  to  test  the  use  of  tracers  for 
detection  of  Rayleigh-Taylor  instability.  Many  of  these  measurements 
provide  valuable  data  for  calibration  of  the  hydrodynamics  codes  upon  which 
we  rely  for  predictions  of  high-gain  target  performance.  As  many  of  the 
current  ideas  for  use  of  tracers  are  proven  and  as  new  ideas  continue  to 
surface,  we  anticipate  an  expanding  role  for  these  methods  in  experiments 
related  to  the  physics  of  laser-matter  interaction  and  to  inertial 
confinement  fusion. 
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